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1  Introduction 


Background 

As  commercial  towboats  navigate  inland  waterways,  physical  forces  in  the  form  of 
waves  and  currents  are  generated.  In  confined  or  restricted  waterways,  these  vessel- 
induced  forces  can  have  a  significant  effect  on  the  surroimding  environment.  Evaluating 
the  impacts  of  navigation-induced  forces  on  the  riverine  environment  has  become  a  topic 
of  serious  discussion  and  intense  debate  in  recent  years. 

Technical  studies  related  to  navigation  effects  investigate  the  characteristics  of  the 
physical  forces  produced  by  moving  vessels  and  determine  the  effects  of  these  forces 
from  both  an  engineering  and  a  biological  perspective.  Bank  stability,  bed  stability, 
sediment  transport,  maneuverability  of  the  vessel  through  the  waterway,  sill  heights  at 
locks,  and  collisions  with  structures  are  a  few  of  the  reasons  evaluation  of  forces  is 
important  from  an  engineering  design  perspective.  Environmental  concerns  range  from 
propeller  jet  turbulence  impacts  on  fish  larvae,  wave  energy  on  aquatic  vegetation,  to 
sediment  resuspension  and  transport  in  critical  shoreline  habitat  areas  and  adjacent 
backwaters. 

Minimizing  shoreline  erosion  can  be  critical  for  the  protection  of  valuable  resources 
whether  they  are  man-made  or  natural.  Shoreline  erosion,  among  other  things,  can  lead 
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to  loss  of  habitat,  can  expose  important  archaeological  sites,  and  can  threaten  the 
structural  integrity  of  flood  control  and  navigation  projects.  This  research  focuses  on  the 
protection  of  banks  from  forces  produced  by  commercial  vessels  underway  on  inland 
navigable  waterways  using  stone  slope  protection. 

Riprap  is  a  classical  solution  to  bank  stability  problems  whether  related  to  wave- 
induced  forces  or  flow-induced  velocities.  Minimization  of  rock  size  can  result  in 
substantial  savings  in  a  bank  protection  project.  Design  guidance  related  to  riprap  design 
for  waves  was  originally  developed  for  protection  of  coastal  shorelines.  Beginning  with  a 
formulation  developed  by  Iribarren  (1948),  Hudson  (1957)  reformulated  an  equation,  still 
widely  accepted  and  used,  that  relates  stone  weight  to  wave  height.  Criteria  for  riprap 
design  in  open  channel  flows  has  been  studied  and  developed  by  many,  including 
Maynord  (1988). 

The  Delft  Hydraulics  Laboratory,  the  Netherlands,  has  conducted  research  and 
produced  design  methodologies  for  riprap  design  in  constricted  waterways  due  to  vessel¬ 
generated  forces.  Site-specific  studies  have  been  conducted  at  the  U.S.  Army  Engineer 
Waterways  Experiment  Station  (WES),  Vicksburg,  MS,  for  the  Termessee-Tombigbee 
Waterway  (Maynord  and  Oswalt  1986)  and  for  locks  and  dams  on  the  Ohio  River.  The 
Ohio  River  studies  for  the  U.S.  Army  Corps  of  Engineers,  Huntington  District  (Martin 
1989),  suggested  that  rather  small  stone  sizes  exhibit  stability  under  reasonable  vessel 
operating  conditions.  Verification  of  this  premise  could  result  in  considerable  savings  to 
inland  navigation  projects  by  reducing  the  size,  and  consequently,  blanket  thickness,  of 
stone  protection  in  heavily  navigated  reaches  of  commercial  waterways.  Savings  have 
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already  been  realized  on  the  Point  Marion  Lock  approach  on  the  Monongahela  River 
based  on  site-specific  studies  recommending  a  reduction  in  the  riprap  size  over  the 
proposed  stone  (Maynord  1989). 

The  U.S.  Army  Corps  of  Engineers  have  recently  been,  and  are  currently,  authorized 
to  construct  major  rehabilitation  projects  on  several  inland  locks  and  dams  with  studies 
underway  on  many  more.  In  the  approach  channels  to  these  locks,  where  the  major 
potential  causal  mechanism  for  bank  erosion  can  be  attributed  to  commercial  tow  traffic, 
a  generalized  approach  is  needed  to  minimize  the  need  for  site-specific  studies  and  to 
insure  safe  and  economical  stone  slope  protection.  Sorensen  (1986),  in  his  investigations 
of  existing  criteria  for  the  design  of  bank  protection  from  vessel-generated  waves, 
concludes,  "significant  improvements  can  be  made  in  the  quality  of  this  design  informa¬ 
tion  so  less  conservative  and  more  economical  designs  can  be  made  with  more 
confidence." 


Objective 

The  main  objective  of  this  research  was  to  develop  a  riprap  sizing  procedure  for 
protecting  bank  slopes  from  tow-induced  forces  generated  by  commercial  traffic  on 
United  States  navigable  waterways.  A  commercial  tow  consists  of  multiple  barge  units 
pushed  by  a  towboat.  The  riprap  design  applications  were  in  essentially  slack  water 
canals  or  natural  channels  with  small  ratios  of  wetted  channel  cross-sectional  areas  to 
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submerged  vessel  cross-sectional  areas.  A  secondary  emphasis  was  placed  on 
presentation  of  existing  methods  for  quantifying  tow-induced  forces. 

Although  not  recognized  as  an  objective  prior  to  this  study,  considerable  knowledge 
was  obtained  regarding  scale  effects  in  physical  model  studies.  Low  Reynolds  values 
during  testing  and  the  test  materials’  natural  angle  of  repose  were  found  to  contribute  to 
stone  instability. 

Approach 

This  research  scrutinized  existing  methods  for  sizing  riprap  by  (a)  the  review  of 
literature,  (b)  analyzing  physical  model  data  from  previous  site-specific  studies,  and 
(c)  supplemental  data  through  additional  physical  model  testing.  Applicability  of  existing 
techniques  and  guidance  were  compared  to  results  obtained  through  the  model  tests. 
Recommendations  for  riprap  sizing  were  made  to  satisfy  research  objectives. 

The  existing  and  supplemental  riprap  stability  tests  were  conducted  at  two  model 
scales  by  exposing  different  gradations  of  model  stone  on  two  bank  slopes  with  repetitive 
passes  of  a  model  tow  (towboat  and  variable  barge  configurations).  The  tow’s  draft, 
speed,  and  sailing  line  were  varied.  Several  water  depths  were  used  during  testing.  All 
tests  were  conducted  in  generally  confined  channels  with  no  ambient  currents. 
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2  Literature  Review 


Existing  Methods  for  Predicting  Physical  Forces 

Definitions 

Vessel-induced  water  motion  or  navigation  effects  can  be  classified  in  three  general 
categories  according  to  the  mechanisms  that  caused  the  effects  (a)  the  displacement 
motion  of  the  vessel  in  the  channel  or  the  primary  wave  system,  (b)  waves  produced  by 
dynamic  ship  theory  or  secondary  waves,  and  (c)  propeller  jets.  Methods  coupling  all  the 
dynamics  of  vessel  motion  that  quantify  the  three-dimensional  (3-D)  forces  produced  do 
not  exist.  Therefore,  the  three  basic  areas  of  vessel-induced  water  motion  are  studied  in 
pieces  based  on  conservation  of  energy  or  momentum,  slender  body  theory,  and  empirical 
approaches  (Blaauw  et  al.  1984). 

The  vessel’s  displacement  motion  produces  long  period  changes  in  water-surface 
elevations  and  some  rather  characteristic  current  patterns  around  and  beneath  the  vessel. 
This  phenomena  can  be  explained  in  much  the  same  manner  as  a  piston  moving  through  a 
cylinder.  As  the  vessel  moves  faster  through  the  channel  the  fluid  must  change  directions 
and  water  levels  to  satisfy  the  principles  of  energy  and  continuity.  Figure  1  shows  the 
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Figure  1.  Schematic  of  primary  navigation  effects 


primary  components  of  this  motion  in  plan  view.  Major  water  level  changes  are 
attributed  to  the  front  or  bow  wave,  drawdown,  and  transversal  stem  wave.  Major 
velocity  components  are  the  return  current  and  slope  supply  flow.  The  strength  of  these 
forces  is  primarily  a  function  of  the  vessel  speed  and  its  size  relative  to  the  channel.  The 
current  beneath  the  vessel  is  a  force  not  shown.  The  shear  forces  and  velocities  generated 
in  this  region  can  be  critical  to  the  channel  bed,  but  do  not  influence  the  stability  of  the 
channel  banks. 
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Secondary  waves  include  the  entire  wave  spectrum  formed  by  the  pressure 
disturbance  on  the  hull  of  the  moving  vessel.  The  vessel  creates  transverse  and  diverging 
waves  intersecting  at  peaks  which  form  a  distinct  pattern  in  deep  water  (Verhey  and 
Bogaerts  1989).  The  shape  of  the  vessel  hull  and  the  speed  of  the  vessel  dictate  the 
magnitude  of  these  waves.  Figure  2  is  a  diagram  showing  this  wave  pattern  in  plan  view. 


Figure  2.  Ship  wave  angles  and  definition  sketch 

A  commercial  tow  is  generally  powered  by  twin  propellers  each  having  diameters  of 
5  to  8  ft  (1 .5  to  2.4  m).  Propeller  jet  forces  can  impact  the  channel  bed  and  the  lower 
portion  of  the  bank  slope  if  the  tow  is  navigating  within  close  proximity  of  the  bank. 
Bank  failure  due  to  propeller  jets  could  occur  in  areas  where  direct  impingement  is 
inevitable  as  in  a  tight  bend,  maneuvering  into  a  lock,  or  mooring  near  a  bank. 
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Since  components  of  the  displacement  forces  and  the  secondary  waves  have  the  most 
impact  on  the  stability  of  channel  banks,  no  further  discussion  is  given  to  propeller  jet 
forces  or  velocities  beneath  the  vessel.  Terms  related  to  navigation  effects  are  found  in 
Appendix  A. 

Displacement  motion 

The  general  physics  related  to  displacement  of  a  tow  moving  through  a  chaimel  can 
be  represented  by  applying  the  principles  of  energy  or  momentum  and  continuity.  These 
theories  are  represented  in  existing  analytical  equations  regarding  prediction  of  these 
forces  particularly  in  the  determination  of  return  current  and  drawdown.  Equations  for 
return  current  and  drawdown  can  be  derived  by  applying  the  Bernoulli  principle  for 
energy  and  conservation  of  mass,  an  approach  used  by  Jansen  and  Schijf  (1953). 
Bouwmeester  et  al.  (1977)  used  the  momentum  equation  rather  than  the  energy  approach. 
Equating  the  energy  along  a  streamline  using  a  point  mid-length  of  the  tow  and  one  in  the 
undisturbed  channel,  and  assuming  superposition,  the  following  results: 


z 


ave 


(1) 


where 

Zave  ~  average  drawdown,  ft  (m) 
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Vj  =  speed  of  tow  relative  to  earth,  :^s  (mps) 

Ur  =  return  current,  fps  (mps) 
g  =  gravitational  constant  =  32.2  ft/s^  (9.81  m/s^) 

Based  on  continuity,  the  equation  follows: 

V.A,  =  (V^  .  u^)A^  (2) 


where 

A(.  =  chaimel  cross-sectional  area  before  drawdown,  (m^) 

=  wetted  cross-sectional  area  of  channel  at  mid-ship  of  the  barges,  ft^  (m^) 
and 

A  (h  -  z\  -  A.  h 

A.  =  P ^  (3) 

h 

Ab  =  cross-sectional  area  of  the  submerged  barges,  fl^  (m^) 
h  =  depth  of  xmdisturbed  water  level,  ft  (m) 

This  standard  method  for  estimating  drawdown  and  return  current  is  called  the  Schijf 
approach  or  method.  It  gives  the  "average"  drop  in  water  surface  and  return  current  near 
mid-ship.  Development  of  these  equations  was  based  on  the  following  assumptions: 

a.  uniform  cross  section 

b.  uniform  return  current 

c.  uniform  water  level  drawdown 
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d.  no  friction 


e.  negligible  ambient  current 

f.  centerline  placement  of  the  vessel 

The  drawdown  and  return  current  are  primarily  a  function  of  vessel  speed  and 
blockage  ratio,  n ,  defined  as  the  ratio  of  A^.  over  At .  The  limiting  speed,  Vl  ,  defined  as 
the  maximum  or  critical  vessel  speed,  was  determined  by  differentiating  Equations  1 
and  2  (Jansen  and  Schijf  1953).  This  speed  is  the  maximum  speed  a  vessel  can  reach  in  a 
confined  channel  of  specified  blockage  ratio  regardless  of  vessel  horsepower. 

Numerous  researchers  have  used  variations  of  this  approach  that  introduce 
coefficients  and/or  modifications  to  the  equations  that  account  for  parameters  such  as  off- 
center  sailing,  irregularity  of  channel  section,  fiiction,  etc.  These  methods  were  based  on 
model  studies  and/or  field  data.  Maynord  (1990)  presented  a  thorough  review  of  the 
many  methods  for  determining  displacement  and  return  velocities.  The  Permanent 
International  Association  of  Navigation  Congresses  (PIANC)  published  guidelines  for 
predicting  navigation-induced  forces  (PIANC  1987).  Equations  in  these  guidelines  are 
composites  of  previous  researchers'  efforts,  particularly  those  from  the  Delft  Hydraulics 
Laboratory,  the  Netherlands.  Figure  3  defines  some  channel  and  vessel  parameters  used 
in  the  following  methods  for  calculating  physical  forces. 
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Figure  3.  Definition  of  geometric  terms  used  in  predictive  equations 

PIANC  (1987)  recommended  the  average  return  current,  u^ ,  and  dra'wdown,  » 
be  estimated  using  the  Schijf  equations.  The  method  is  valid  for  waterway  topwidths, 
bt ,  to  vessel  beam  widths,  Bj ,  between  2  and  12.  They  also  recommended  adding  an 
empirical  coefficient,  aj ,  that  is  a  function  of  the  ratio  of  to  Vl  ,  to  the  energy 
equation. 

Rewriting  Equation  1  to  include  aj  gives 
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Bouwmeester  et  al.  (1977)  describes  the  maximum  depression  of  the  water  level  as 
"the  combined  effect  of  the  primary  wave  and  the  secondary  wave  system."  The  deepest 
trough  occurs  between  the  bow  and  stem.  They  define  the  transverse  stem  wave  as  "the 
slope  of  the  water  surface  due  to  deceleration  of  the  return  flow  at  the  stem."  They 
observed  that  as  the  water  level  depression  increased,  the  transverse  stem  wave  had  the 
appearance  of  a  moving  hydraulic  jump.  Furthermore,  "under  certain  conditions  the  stem 
transverse  wave  was  intensified  by  secondary  waves,"  particularly  when  the  cusp  line  of 
the  secondary  waves  intersected  with  the  crest  of  the  transverse  stem  wave. 

The  transversal  stem  wave,  ,  was  empirically  determined  as  a  function  of  the 
drawdown  in  PIANC  (1987)  and  is  defined  by  Figure  4  as  the  maximum  drawdown.  It 
can  be  determined  for  a  push-tow  unit  (towboat  and  barge  unit)  by 


z 


max 


C  Z 


ave 


(5) 


where 


c  =  1.2  +  SxlO'^F^ 


(6) 
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Figure  4.  Definition  sketch  of  water  level  terms  for  displacement  motion 


The  vessel  Froude  number,  F,, ,  is  defined  as 


F 


h 


V 

_ s 

Vgii 


(7) 


and  Yt  is  defined  as 


=  0.5bj-y-  O.Shcota 


(8) 


bt  =  top  width  of  channel,  ft  (m) 

y  ==  distance  from  centerline  channel  to  centerline  tow,  ft  (m) 
L3  =  length  of  ship  at  waterline,  ft  (m) 
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a  =  angle  slope  makes  with  horizontal 


The  depth,  h ,  in  a  trapezoidal  or  irregular  channel  can  be  approximated  using  Aj./bt . 
From  a  graphical  solution  in  Blaauw  et  al.  (1984)  the  following  linear  relationship  was 
obtained  for  z^ax  as  a  function  of  y  (the  eccentricity),  topwidth,  and  average  drawdown. 


z 


max 


\ 

+  1.2 
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(9) 


Based  on  a  multivariate  regression  of  prototype  data  collected  on  the  Illinois  River, 
Bhowmik,  Demissie,  and  Guo  (1982),  developed  two  empirical  formulas  for  predicting 
wave  heights  as  a  function  of  tow  speed.  The  equation  for  drawdown,  z ,  was  based  on 
the  following: 


\0.81 
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(10) 


where 
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V 

F  =  ^ 

The  variable,  S  ,  is  the  distance  from  the  wave  gage  to  the  sailing  line  and  d  is  the  draft 
of  the  barges  (Figure  3).  The  tow  speeds  measured  for  this  analysis  varied  from  4  to 
20  ^s  (1.2  to  6.1  mps),  and  blockage  ratios  were  quite  variable  as  well,  from  15  to  230. 
The  maximum  measured  drawdown  was  0.7  ft  (0.2  m). 

The  front  wave,  hf ,  is  the  total  peak  to  following  trough  of  the  wave  crest 
preceding  the  bow  to  the  trough  or  drawdown  (Figure  4).  The  relationship  given  in 
PIANC  (1987)  for  computing  this  value  is 


h,  =  0.1  z  +  z 

I  ave  max 


(13) 


From  a  graphical  solution  in  Blaauw  et  al.  (1984),  the  following  linear  relationship  was 
obtained 
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+  1.2 
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Fuehrer,  Romisch,  and  Engelke  (1981)  gave  equations  for  calculating  the  maximum 
change  in  water  level  for  vessels  traveling  less  than  critical  speed  and  for  those  traveling 
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greater  than  critical  speed.  They  stated  that  the  primary  wave  system  plays  the  dominant 
role  in  the  subcritical  vessel  speed  range.  As  the  vessel  increases  its  speed  beyond  this 
critical  level,  the  secondary  wave  motion  becomes  important.  Their  wave  height  equation 
for  vessels  in  the  subcritical  range  actually  represented  the  front  wave  and  drawdown, 
such  that 


(15) 


The  constant,  C,  is  a  function  of  channel  parameters  including  blockage  ratio,  channel 
depth,  and  chaimel  width. 

Based  on  model  data  collected  at  WES  and  the  Delft,  and  field  studies  conducted  by 
the  Louisville  District,  Maynord  and  Siemsen  (1991)  recommended  a  modified  Schijf 
approach  that  introduces  not  only  an  asymmetric  sailing  line,  but  also  provides  a 
mathematical  decay  function  for  return  currents  in  channels  with  higher  blockage  ratios. 
It  was  concluded  that  in  confined  channels  where  blockage  ratios  are  small,  a  uniform 
distribution  of  return  currents  is  a  reasonable  assumption.  Even  when  the  tow  is  sailing 
off  the  channel’s  center  line,  the  lateral  distributions  are  miiform  even  if  the  relative 
magnitudes  are  different  on  one  side  of  the  chaimel  versus  the  other.  In  larger  rivers  the 
strength  of  the  return  current  decays  with  distance  from  the  tow  (Maynord  and  Siemsen 
1990).  Based  on  his  recent  analysis  of  model  and  prototype  data,  Maynord  (1996) 
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presented  a  method  for  calculating  maximum  return  current,  u,^ ,  and  maximuTn 
drawdown,  ,  as  a  function  of  distance  from  the  vessel  for  blockage  ratios  up  to  52. 

The  composite  field  of  currents  formed  from  displacement  motion  include  bow 
currents,  lateral  currents,  return  currents,  flow  beneath  the  vessel,  and  slope  supply 
currents.  The  return  current  and  slope  supply  flow  have  the  most  impact  on  the  bank. 

The  average  maximum  return  current  in  a  channel,  Un„,  is  defined  by  Maynord 
(1996)  as  a  function  of  the  Schijf  return  current,  u^,  and  a  ratio  of  the  vessel  speed  to 
limiting  speed  such  that 


=  u, 
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PIANC  (1987)  defined  the  maximum  return  current  velocity,  u^n, ,  as 
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The  coefficient,  c ,  according  to  PIANC  (1987)  is  defined  by  Equation  6.  From  a 
graphical  solution  in  Blaauw  et  al.  (1984),  a  similar  linear  relationship  existed  for  return 
current.  Equation  9,  and  drawdown,  such  that 
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According  to  Bouwmeester  et  al.  (1977),  when  ships  are  traveling  at  high  speeds 
near  the  bank,  a  translatory  wave  or  moving  bore  is  formed  and  propagates  at  a  speed 
relative  to  the  water  equal  to  the  sum  of  the  return  current  and  ship's  speed.  They  called 
the  high  velocities  near  the  banks,  which  resulted  from  water  being  restored,  the  "slope- 
supply  flow." 

According  to  Blaauw  et  al.  (1984),  the  maximum  currents  in  the  stem  wave,  or  slope 
supply  current  can  be  estimated  by 
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Yr  =  specific  weight  of  stone,  Ib/ft^  (kg/m^) 
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Yw  =  specific  weight  of  water,  Ib/ft^  (kg/m^) 

D50  =  diameter  of  stone  protection  for  which  50  percent  is  lighter  by  weight,  ft  (m) 


Gradients 

Gradients  in  the  front  and  transversal  waves  generate  external  load  factors  causing 
internal  loadings.  According  to  Blaauw  et  al.  (1984),  the  gradient  in  the  longitudinal  (x) 
direction  was  of  minor  significance  to  the  relationship  between  internal  load  and  strength, 
In  the  case  of  revetment  blocks,  the  gradients  in  the  transverse  (y)  direction  can  be 
important  in  designing  filter  layers  and  the  gradients  in  the  z  direction  can  cause  uplift. 

Internal  hydraulic  gradients  in  the  y  direction  can  be  represented  by  the  product  of 
the  external  load  from  the  height  of  the  front  wave,  gradient  of  the  front  wave,  if,  and 
ship  speed  as  follows 


h  -  h 
^  dt  ^  dx  dt 
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The  gradient,  i^^g^  5  of  the  transversal  stem  wave  is  defined  as 
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The  value  of  has  a  limiting  value  of  0.1  to  0.15.  The  factor  Zq  is  a  function  of 
channel  top  -width,  ,  and  eccentricity,  y ,  such  that 


—  =  0.04  -  0.158 


(24) 


Secondary  waves 

Boat  waves  can  be  classified  according  to  criteria  in  the  Shore  Protection  Manual 
(1984).  Drawdown  will  almost  certainly  be  classified  as  a  shallow  water  wave;  whereas 
secondary  waves,  with  periods  of  one  to  five  seconds,  can  approach  a  deep  water 
classification  and  may  be  either.  Based  on  the  same  criteria,  the  waves  formed  by  tows 
in  channels,  whether  short  period  or  long  period,  are  generally  nonbreaking,  except  in 
shallow  near-shore  areas. 

Figure  5  contains  time-history  schematics  of  the  waves  generated  by  a  moving  vessel 
in  a  confined  reach.  The  first  response  contains  the  unfiltered  long  and  short  period 
waves;  the  second,  the  filtered  long  period  wave  or  drawdo-wn;  and  the  third,  the  filtered 
short  period  or  secondary  waves. 

It  has  been  observed  by  the  author  and  others  (Fuehrer,  Romisch,  and  Engelke  1981) 
that  secondary  waves  are  negligible  at  low  vessel  speeds  and  high  blockage  ratios.  At 
faster  vessel  speeds,  secondary  waves  can  become  more  significant  and  their  magnitude 
approach  or  exceed  the  magnitude  of  the  maximum  drawdown.  Generally,  secondary 
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Figure  5.  Typical  time-history  wave  response  for  a  tow  in  a  confined  channel 
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waves  have  higher  frequency,  lower  amplitude,  and  last  longer  than  the  shallow  water 
drawdown  effects.  For  vessels  in  larger  channels  and  recreational  craft,  the  diverging 
coincidental  peaks  can  dominate  the  wave  spectrum. 

Verhey  and  Bogaerts  (1989)  presented  a  predictive  equation  for  secondary  waves 
that  related  the  wave  height  to  the  sailing  line,  vessel  speed,  and  water  depth.  The 
coefficient,  a, ,  regards  the  type  of  vessel  and  its  draft.  The  equation  for  estimating 
maximum  secondary  wave  height,  H; ,  for  the  incidental  wave  peak  is 
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where: 

=  coefficient  regarding  hull  type  and  draft 
aj  =  coefficient  experimentally  determined  to  be  4.0 
s  =  distance  between  ship's  side  and  water's  edge,  ft  (m) 

Values  of  were  determined  experimentally  as  follows: 
a  I  =  1 .0  for  loaded  conventional  inland  motor  vessels 
a,  =  0.35  for  empty  conventional  motor  vessels 
=  0.5  empty  European  barges 

Verhey  and  Bogaerts  (1989)  did  not  provide  a  value  of  ttj  for  loaded  barges  stating  that 
the  wave  heights  for  this  type  were  negligible.  Blaauw  et  al.  (1984)  used  the  same 
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equation,  experimentally  determining  values  of  and  ocj  for  loaded  push-tows  of  0.80 
and  2.67,  respectively.  In  PIANC  (1987)  and  Boeters,  van  der  Knaap,  and  Verheij  (1995) 
the  coefficients  for  ttj  and  a2  were  1.0  and  4.0,  respectively. 

Equation  25  is  valid  for  nonbreaking  waves  with  Hj/h  less  than  approximately  0.6 
and  vessel  Froude  numbers  less  than  0.7  (Verhey  and  Bogaerts  1989).  For  vessels  in  this 
Froude  number  range,  the  angle  of  wave  propagation  normal  to  the  bank,  p  ,  is  given  by 
PIANC  (1987),  Verhey  and  Bogaerts  (1989),  Sorensen  (1986),  and  others  as  54.7  deg. 
The  companion  angle,  0  ,  between  the  direction  of  propagation  and  the  vessel  sailing  line 
is  35.3  deg.  Figure  2  shows  these  angles. 

From  deep  water  wave  theory,  Verhey  and  Bogaerts  (1989),  Sorensen  (1986),  and 
others  suggested  that  the  following  average  wave  period,  Tj  ,  and  wave  length,  L^;  , 
correspond  to  Hj  . 


T. 


2  7i:(cos  0 ) 


g 


(26) 


T.^ 

^wi  "  =  2ir(cos^0)— 


2  7C 


g 


(27) 


As  vessel  Froude  numbers  approach  1 .0  the  transverse  wave  dominates  the  diverging 
waves  so  that  the  wave  energy  is  concentrated  in  a  single  wave  near  the  stem.  As  the 
Froude  number  increases  above  approximately  0.7,  0  decreases  from  35  deg  to  0  deg 
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(Sorensen  1986).  This  phenomenon  is  demonstrated  when  the  transversal  stem  wave 
moves  as  a  hydraulic  jump  traveling  perpendicular  to  the  bank. 

Weggle  and  Sorensen  (1986)  studied  secondary  waves  and  developed  a  wave  height 
prediction  model  as  a  function  of  the  vessel  speed,  distance  from  the  sailing  line,  the 
water  depth,  and  the  vessel  displacement.  There  were  two  general  equations  (a)  one  for 
vessel  Froude  numbers  less  than  0.55  and  (b)  one  for  Froude  numbers  greater  than  0.55. 
The  equations  required  empirical  coefficients  for  each  specific  hull  form. 

Based  on  the  same  prototype  data  collected  on  the  Illinois  River  for  determination  of 
an  equation  for  drawdown,  Bhowmik,  Demissie,  and  Guo  (1982)  also  developed  an 
empirical  formula  for  predicting  the  maximum  wave  height,  H; ,  of  the  secondary  waves. 
It  is  given  as 


Hi 

—  =  0.133  F. 
d  ^ 


(28) 


The  draft  Froude  number,  F^  ,  is  defined  as 


F 


d 


The  maximum  measured  wave  height  was  1.1ft  (0.3  m). 


(29) 
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Current  Riprap  Design  Criteria 


Numerous  equations  have  been  developed  that  relate  velocity,  shear  stress,  or  wave 
height  to  a  stable  riprap  design.  This  study  regards  specifically  the  design  of  stone  slope 
protection  for  tow-induced  forces  which,  as  described  above,  is  a  complex  combination  of 
waves  and  currents.  A  detailed  literature  review  follows  highlighting  the  most  pertinent 
formulations,  including  general  wave-based  equations,  vessel-related  equations,  and 
formulations  based  on  shear  and/or  velocity. 

Forces  on  stone  slope 

Van  Gent  (1995)  and  Kobayashi  and  Greenwald  (1988)  both  described  three  basic 
forces  that  are  applied  to  the  slope  under  wave  attack:  the  drag  force,  ,  acting  parallel 
to  the  bank  in  the  direction  of  the  velocity;  the  inertia  force,  Fj ,  acting  parallel  to  the 
bank;  and  the  lift  force,  Fl  ,  acting  normal  to  the  slope.  The  forces  were  balanced  against 
the  buoyant  weight  of  the  stone,  Wg .  They  are  defined  as 

Fd  =  (30) 

Fi=Y,C,k,D=^  (31) 


25 


(32) 


F,  =  0.5  Y. 


W3  =  kj(Y,  -  yJd'  (33) 

k]  =  volume  shape  factor,  7c/6  for  spheres 

k2  =  area  shape  factor,  7t/4  for  spheres 

Cq  =  drag  coefficient 

Cj  =  inertia  coefficient 

Cl  =  lift  coefficient 

D  =  diameter  of  stone,  ft  (m) 

u  =  directional  velocity,  ^s  (mps) 

^  =  acceleration 

fit 


For  stability,  the  sum  of  the  forces  parallel,  Fp  ,  to  the  bank  over  the  sum  of  the  forces 
normal,  F^ ,  to  the  bank  are  balanced  against  the  internal  fnction  of  the  slope,  such  that 


— ^  tan  ({) 

N 

where 

4)  =  natural  angle  of  repose 


(34) 
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and 


X^Fp  =  Fd  +  Fi  -  Wg  since 


(35) 


=  Wgcosa  -  Fg 


(36) 


Figure  6  is  a  force  diagram  showing  the  major  force  components  and  angles. 


Figure  6.  Forces  on  surface  layer  of  riprap  embankment 
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Kobayashi  and  Greenwald  (1988)  evaluated  graded  stone  on  an  impermeable  slope 
and  Van  Gent  (1995)  on  a  permeable  breakwater.  Van  Gent  (1995)  described  an 
additional  force  at  the  wave  jfront  resulting  from  the  pressure  gradient  related  to  the  free 
surface  slope.  This  force  determines  the  direction  of  stone  motion. 

Iribarren  (1948)  balanced  the  drag  and  lift  forces  against  the  internal  friction  such 
that 


R 

(N  -  P,) 


<  tan({) 


(37) 


where 

R  =  parallel  force  resisting  motion 
N  =  normal  force  resisting  motion 
Pi  =  normal  force  lifting  stone 
The  resistant  and  normal  forces,  respectively,  are 


R  =  Wg  sina 


(38) 


N  =  Wg  cosa 


(39) 
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Pi 


(40) 


where 

H  =  height  of  attacking  wave 
k,  =  coefficient 


Substituting  into  Equation  37  yields 


Wg  sina 

-  ^  tan(j)  (41) 

Wgcosa  -  k,Y„D^H 


Due  to  the  complexity  of  the  acceleration  term  in  the  inertia  force,  Hudson  (1957) 
combined  it  with  the  drag  force  resulting  in  a  total  force,  .  Hudson  dismissed  the 
internal  friction  stating: 


For  breakwaters  constructed  by  dumping  or  placing  armor  units 
essentially  pell-mell,  the  forces  resisting  displacement  are  the  buoyant 
weight  of  the  individual  units  and  the  fiiction  between  units.  Except  for 
isolated  instances  where  wedging  action  is  involved,  fiiction  between 
armor  units  can  be  neglected,  and  the  principal  resisting  force  for  pell- 
mell-constructed  cover  layers  can  be  assumed  to  be  Wg  . 


Then  for  stability  of  a  rubble-mound  breakwater  subjected  to  breaking  waves  Hudson 
(1957)  used. 
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(42) 


F  <  W 

q  ~  B 


(43) 


where 

Cq  =  total  force  coefficient 


Though  Hudson's  equation  was  based  on  breaking  waves,  he  states, 


The  forces  that  tend  to  displace  armor  units  from  breakwater  slopes  when  the 
waves  do  not  break,  or  break  only  partially,  are  not  the  same  as  Aose  forces  that 
result  from  breaking  waves,  nor  do  they  act  in  the  same  directions.  However,  the 
order  of  magnitude  of  the  non-breaking  wave  forces,  and  the  effects  of  these 
forces  on  the  stability  of  rubble-mound  breakwaters,  should  be  approximately  the 
same  as  those  caused  by  breaking  waves. 


Hedar  (1960)  added  a  moment  to  his  force  diagram  applying  the  resultant 
hydrodynamic  force,  Fj ,  at  some  distance  from  the  stone’s  center  of  gravity.  The  angle, 
P  ,  determined  the  direction  of  the  force.  Hedar  distinguished  between  uprushing  waves 
and  downrushing  waves  in  his  theoretical  development  of  the  forces,  such  that  for 
uprushing 


cosP  +  tancf)  sinP  +  ^ 


W  ^  F 

VV  B  ^  q. 


coscj) 


tanji)  cos  a  +  sina 


(44) 
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The  coefficient,  c ,  relates  to  the  distance  from  the  center  of  gravity  the  force,  F^- ,  is 
applied;  if  c  is  assumed  to  be  1,  the  force  is  applied  to  the  outer  edge  of  the  stone.  The 
total  force  terms  can  be  replaced  with  components  along  the  slope  and  normal  to  it  such 
that 


Fp  =  F.j,cosP 


(45) 


(46) 

Substituting  these  equalities  into  Equation  44,  assuming  c  =  1,  and  rearranging  terms 
according  to  parallel  and  normal  forces,  yields 


/ 

k 


cos(j)^ 


Wg  sina 


Wgcosa  -  Fj^ 


<  tan(l) 


(47) 


From  this  formulation,  it  can  be  shown  that  Hedar's  equation  for  uprushing  waves  is 
comparable  to  Equation  34,  except  that  the  parallel  force  component  includes  a 
coefficient  that  is  a  function  of  (j) . 

A  similar  equation  for  downrushing  waves  can  be  formulated  from  Hedar's  work 
such  that 
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(48) 


Wg  cosa  - 


<  tancj) 


For  open  channel  environments,  the  critical  shear  stress  for  noncohesive  sediments 
was  obtained  by  equating  moments  of  the  gravity  and  drag  forces  about  the  support  point. 
Though  considered  important,  the  lift  force  is  often  neglected  and  assumed  to  be 
accoimted  for  implicitly  in  the  formulation  (American  Society  of  Civil  Engineers 
(ASCE)  1977). 

The  critical  shear  stress,  on  a  sloping  bed  was  given  as 

~  Yw)d  cosa(tan(i)  -  tana)  (49) 


For  a  horizontal  bed,  this  becomes 


-  yjDtancj) 


(50) 


where 

Cl  =  coefficient 
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Use  of  this  in  a  dimensionless  shear  stress  form,  plotted  against  a  critical  Reynolds  value 
results  in  the  well-known  Shields  diagram  (ASCE  1977). 

Wind  wave  equations 

Many  riprap  sizing  methods  for  protection  from  waves  are  based  on  analytical 
methods  supported  by  laboratory  and  field  studies  in  which  wave  height  is  related  to 
stone  size  (Hudson  1957;  Iribarren  1948;  Ahrens  1989;  and  others).  Van  Gent  (1995)  and 
Kobayashi  and  Greenwald  (1988)  both  used  numerical  approaches  to  calculate  wave 
forces  and  slope  stability,  and  compared  their  results  to  laboratory  scale  tests. 

Many  empirical  wave  based  equations  use  a  dimensionless  variable  first  defined  by 
Hudson  (1958)  as  the  stability  number,  .  Depending  on  the  assvunptions  used  by  the 
researcher  for  a  particular  application,  the  stability  number  can  be  a  function  of  a  number 
of  variables.  The  stability  number  is 


N 


S 


H 

Ad 


(51) 


The  weight  of  the  stone,  ,  is  related  to  the  diameter  as  follows: 


W  =  k,  D  ^  Y 

r  1  •  r 


(52) 
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Where,  kj ,  defined  above  as  the  volume  shape  factor  is  7u/6  for  spherical  stones  and  1  for 


cubic  stones.  The  stability  number  as  a  function  of  is  defined  as 


N. 


Aw 

r 


(53) 


Iribarren's  formula  (1948)  for  incipient  motion  of  stone,  including  the  friction  angle 
and  based  on  the  balance  of  forces  in  Equation  41,  yields 


K  Y,tan(})^H^ 


W  = 


(tan(j)cos  a  -  sin  a)^ 


1l 
I  Yw 


-  1 


(54) 


Then,  the  stability  coefficient,  ,  for  Iribarren's  formula  is 


tancbcosa  -  sin  a 
N  =  - i - 

®  — xl/3 


(K)'"^  tan(j) 


(55) 


According  to  a  graphical  representation  given  in  Hudson  (1957),  values  of  Iribarren's 
coefficient  K  are  approximately  0.013  for  a  2H:1V  slope  and  0.030  for  a  3H:1V  slope. 

For  a  tan  (})  of  1,  this  results  in  values  of  Nj  of  1.9  and  2.0,  respectively. 

Hudson  felt  that  the  formulation  of  Iribarren's  equation  placed  too  much  emphasis  on 
the  firiction  force.  This  meant  that  for  steeper  slopes,  small  variations  in  tan  (j)  caused 
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For  a  2H;  1 V  slope,  Nj  for  nonbreaking  waves  on  a  rubble-mound  breakwater  is 
approximately  1.86  and  for  a  3H:1V  slope,  it  is  2.13. 


Hedar  (1960)  expanded  on  Iribarren's  and  Hudson's  work  concluding  that  the 
stability  formula  should  include  more  variables.  Hedar  agreed  with  Hudson  that  wave 
period  does  not  effect  stability  and  for  non-breaking  waves,  depth  of  water  in  front  of  the 
structure  is  not  significant.  From  experimental  results  to  obtain  empirical  coefficients 
related  to  the  force  balances  in  Equations  47  and  48,  the  final  stability  formulas  were 
developed  for  uprushing  and  downrushing  waves.  The  stable  stone  size  was  determined 
to  be  a  function  of  the  angle  of  repose,  angle  of  attack,  and  the  application  point  of  the 
shear  stress.  For  uprushing  breaking  waves 


K  (d. 

up  b 

+  H) 

A  —-32 

f,  14.83h1 
log - 

2 

(tan(j)cosa  +  sina) 

6 

1  D  J 

(57) 


Hedar  stated  for  impervious  slopes  the  coefficient,  Kyp ,  was  18.0  and  the  depth  of  the 
breaking  wave,  db ,  is  assumed  to  be  approximately  1.3  times  H.  The  stability  number 
becomes 


N  =  —-0.89 
®  6 


log 


14.83  H 
D 


(tan^Jcosa  +  sina) 


(58) 
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For  downrushing  waves  at  reflection  and  at  breaking 


—  •32(tan(bcosa  -  sina) 

N  -  -  (59) 

K, 

down 

The  coefficient,  ,  was  given  as  a  function  of  (j)  minus  a  and  was  graphically 
presented  in  Hedar's  paper  (1960)  for  both  pervious  and  impervious  slopes. 

Bhowmik  (1976)  conducted  a  field  investigation  on  the  stability  of  existing  shore 
protection  on  Carlyle  and  Rend  Lakes  in  Illinois  from  wind  waves.  From  his  review  of 
the  literature,  Bhowmik’ s  recommended  design  equation  is 


W 


so 


ky,h^ 

yI  (cos  cc  -  sin  a)^ 


(60) 


where 

W50  =  stone  weight  which  is  50  percent  lighter  by  weight,  lb  (kg) 
K  =  coefficient 

An  "average"  coefficient  was  selected  as  follows: 


—  =  0.388 

•  w 


(61) 
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Reformulating  in  terms  of  the  stability  number,  this  becomes 


Nj  =  5.438  (cosa  -  sina) 


(62) 


Ahrens  (1989)  also  presented  a  formulation  relating  the  stability  number  to  the 
significant  wave  height,  .  To  account  for  instabilities  that  may  result  jfrom  waves  in 
excess  of  the  significant  wave  height,  a  safety  factor  was  employed  such  that 


=  1.14  (cot 


(63) 


This  resulted  in  a  zero  damage  stability  number  of  1.28  for  a  2H:1V  slope  and  1.37  for  a 
3H:1  V  slope.  Since  this  formulation  resulted  in  too  conservative  rock  sizes,  he  ultimately 
adopted  an  approach  which  incorporated  a  surf  parameter  (Van  der  Meer  and  Pilarczyk 
1984). 

The  Delft  Hydraulics  Laboratory  has  published  numerous  articles  and  reports 
regarding  the  design  of  bank  protection.  Van  der  Meer  and  Pilarczyk  (1984)  conducted 
over  200  model  tests  on  rubble  mound  slopes  tmder  random  wave  attack  to  determine  the 
importance  of  the  number  of  waves,  wave  period,  stone  grading,  and  core  permeability. 
Their  tests  showed  that  erosion  damage  was  definitely  dependent  upon  wave  period.  This 
result  was  contrary  to  other  author's  results  (Hudson  1957;  Hedar  1960;  and  Broderick 
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and  Ahrens  1982).  Van  der  Meer  and  Pilar2yck  (1984)  developed  the  following  rela¬ 
tionships  for  Nj  that  included  the  surf  similarity  parameter,  ^  ,  and  the  damage  level,  S  . 
The  surf  similarity  parameter  at  zero  damage,  ,  is  defined  as  follows: 


tana 


27tH. 


gT. 


\  1/2 


(64) 


is  the  wave  period  associated  with  "no  damage,"  and,  according  to  Ahrens  (1989),  is 
generally  taken  to  be  when  S  is  between  1  and  3  stones  eroded.  "Failure"  is  defined  as 
S  between  8  and  17. 


Where 


(65) 


A  is  the  area  eroded  on  the  profile  in  Figure  7.  Both  permeable  and  impermeable  slopes 
were  tested.  Soundings  and  a  profiler  were  used  to  determine  S  . 
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Figure  7.  Profile  of  bankline  eroded  by  wave  action 


Then  the  following  stability  formulae  were  used  for  impermeable  slopes: 
For  breaking  waves  with  <  2.5  -  3.5  , 


^  Ad 


=  4.4 


50 


^_S_^ 

V^, 


0.22 


'-0.S4 


(66) 
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where 


N  =  number  of  waves 


For  nonbreaking  waves  with  >  2.5  -  3.5  and  cot  a  ^  3 


_ / 

=  1.25  \/cot  a 

\ 


For  nonbreaking  waves  with  >  2.5  -  3.5  and  cot  a  ^  3 


(67) 


=  1.25  v/3 


S 

V^, 


1/6 


rO.l 


(68) 


In  these  equations,  the  stones  were  assumed  to  be  cubic  not  spherical.  Figure  8 
schematically  represents  the  surf  similarity  parameter  for  riprap  slopes  as  defined  by 
PIANC  (1987). 


Figure  8.  Approximate  values  of  the  surf  similarity  parameter  as  defined  in  PIANC  (1987) 
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The  stability  formula  given  for  Avind  waves  in  PIANC  (1987)  in  terms  of  the 
dimensionless  stability  number  is 


H 


Ad 


<  2.25^ 


-0.5 


50 


(69) 


Vessel  equations 

Fuehrer,  Romisch,  and  Engelke  (1981)  conducted  extensive  tests  regarding  bank  and 
bottom  protection  against  vessel  waves  and  propeller  jets.  They  determined  that  return 
current  was  only  important  on  the  lower  portion  of  the  slopes  particularly  on  unprotected 
banks.  For  riprap  stability  in  a  trapezoidal  navigation  canal,  they  recommended 
protection  from  the  primary  wave  system,  particularly  the  front  wave.  For  initiation  of 
motion,  the  stability  number  is 


=  3.2  (cot a)*'^ 


(70) 


The  authors  felt  a  safety  factor  was  necessary,  so  recommended 


=  1.05  (cot a) 


(71) 
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Additionally,  the  authors  found  that  relative  scour  was  a  function  of  the  number  of  runs 
(vessels  passing  the  embankment)  and  after  approximately  100  runs,  a  quasi-steady  state 
condition  was  reached. 

Blaauw  et  al.  (1984)  presented  a  total  methodology  including  formulas  for 
drawdown,  secondary  waves,  and  return  current,  as  well  as  guidance  on  the  design  of 
riprap  due  to  vessel  motion.  Tests  were  carried  out  at  prototype  and  at  reduced  scales 
mainly  concerning  push  tows  in  canals  having  a  bank  slope  of  4H:1  V  and  for  riprap  on  a 
granular  filter.  Details  regarding  the  prototype  testing  on  the  Hartel  Canal  can  be  found 
in  Pilarczyk  (1984).  Blaauw  et  al.  (1984)  conducted  their  study  specifically  to  address 
tow-induced  motion  and  riprap  design. 

Blaauw  et  al.  (1984)  discussed  navigation-induced  effects  in  three  primary  categories 
of  water  motion  already  described.  They  stated  that  the  upper  bank  revetment  is  affected 
by  the  secondary  waves  and  transversal  stem  wave,  while  the  lower  slope  and  channel 
bottom  are  affected  by  the  propeller  jet  and  return  current.  The  subsoil  is  effected  by  pore 
pressure  variations  due  to  the  front  wave,  drawdown,  and  secondary  waves. 

The  authors  evaluated  effects  of  return  current  using  two  approaches.  The  first,  a 
shear  stress  approach,  stated  that  the  maximum  shear  stress,  ,  occurs  at  maximum 
return  current,  Uj^  ,  such  that 

•'.  =  (72) 
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where 


p  =  density  of  water,  slug/ft^  (kg/m^ ) 

The  equation  for  the  coefficient  of  drag,  ,  is  given  according  to  the  Schlicting  formula 
for  rough  plates  (Schlichting  1979). 


C 


fr 


'  \  -2.S 

2.87  +  1.58  log  — 
k 

V  s/ 


(73) 


where 


and 


X  = 


(74) 
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50 


(75) 


The  bottom  roughness,  Iq ,  is  defined  as  an  "average"  value  using  the  multiple  of  four  to 
account  for  surface  irregularities.  The  velocity,  u,. ,  is  the  average  return  current,  Lj  is 
the  total  ship's  length  and  x  is  the  distance  a  water  particle  moves  due  to  u^ .  According 
to  the  authors,  the  maximum  return  current  occurs  at  0.25  to  0.35  of  the  ship's  length  from 
the  bow. 
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Using  Shield's  criteria  for  stable  rock  design  and  taking  advantage  of  the  shear  stress 
formula  in  Equation  72,  Blaauw  et  al.  (1984)  uses  the  flow  parameter,  ij; ,  given  as 


r 

ij;  =  - E - 

Pg^Djo 


(76) 


where 


k,  =  cosa 

Riprap  stability  criteria  for  this  method  was  given  as 
t|f  ^  0.03  =  no  transport  of  riprap 
t|;  ^  0.06  =  small  transport  of  riprap 
i|r  >  0.06  =  transport  increases 

For  their  recommended  design  for  return  current,  the  authors  used  the  well-known 
Isbash  equation. 


1- 


tan^g 

tan^d) 


(77) 
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Using  the  return  current  near  the  bank,  u^^  if  known,  the  slope  was  stable  for  values  of  kj 
from  1.2  to  1.5. 

An  equation  was  given  for  zero  damage  (initiation  of  motion)  against  the 
transversal  stem  wave  for  a  push  tow  sailing  near  a  bank  based  on  a  stability  number 
called  the  "characteristic  stem  wave  parameter"  by  Blaauw  et  al.  (1984). 


^  <  2.3 


Ad 


50 


(79) 


This  equation  is  valid  for  a  4H:1  V  slope  and  assumes  cubic  stones. 

Blaauw  et  al.  (1984)  found  during  testing  that  no  material  was  transported  due  to 
secondary  waves  for  the  following  condition,  where  the  cot  a  =  4  and  p  =  55  deg. 


3 

Ad 


so 


<  3.0 


(80) 


The  use  of  gradients  and  slope  supply  flow  was  only  applied  when  trying  to  predict 
transport  of  the  stones.  Blaauw  et  al.  (1984)  gave  equations  for  determining  the  lowest 
level  of  transport  below  the  undisturbed  water  by  the  stem  wave  and  the  secondary 
waves.  In  this  manner,  they  recommended  a  duel  riprap  cover  design  protection  (upper 
for  waves,  lower  for  return  current). 
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In  their  study  of  ship  waves  and  slope  protection,  Verhey  and  Bogaerts  (1989) 
examined  the  validity  of  the  Hudson  approach  to  ship  waves.  They  stated  that  ship 
waves,  propeller  jet,  return  current,  and  water  level  depression  attacked  the  banks. 
Verhey  and  Bogaerts  distinguished  ship  waves  from  secondary  waves  as  the  interference 
peaks  of  the  diverging  and  transverse  waves  and  should  be  determined  by  Equation  25. 
Adjusting  the  wave  height,  H; ,  by  the  waves’  angle  of  incidence  normal  to  the  slope, 
p,  model  studies  indicated  that  riprap  protection  was  stable  for 


H.cosp®^ 

-  <  1.8  to  2.3 
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(81) 


However,  their  recommended  equation  incorporated  the  surf  similarity  parameter,  ^  ,  to 
account  for  wave  length,  such  that 


N  <  <  2.25 (cos g  +  sing) 

Here  the  surf  similarity  parameter  is  given  as  a  function  of  incidental  wave  length, 
instead  of  period  as  in  Equation  64. 


(82) 


tang 


[h,] 

0.5 

(83) 
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Furthermore,  the  relationship  is  only  valid  for 


^  <  (0.05  cota)'°-* 


(84) 


and  was  only  confirmed  for  a  slope  of  4H:1V.  The  authors  stated,  "the  equation  seems 
too  optimistic"  for  a  slope  of  3H;  1 V.  Based  on  a  graphical  representation  for  full-scale 
tests  on  2H:1V  and  3H:1V  slopes,  values  of  (Hj  not  adjusted  by  angle  of  attack) 
between  1.6  to  2.0  resulted  in  only  a  small  displacement  of  stones. 

Verhey  and  Bogaerts  (1989)  presented  both  deterministic  and  probabilistic  design 

criteria.  These  studies,  similar  to  those  conducted  by  Blaauw  et  al.  (1984)  at  the  Delft 

Hydraulics  Laboratory,  presented  design  methodologies  for  bank  protection  based  on 

both  approaches.  According  to  the  authors. 

In  the  deterministic  approach  a  dominant  design  condition  is  selected. 

On  the  basis  of  this  condition  the  dimensions  of  the  protection  and  filter 
layer  are  determined  for  the  criterion  of  "initiation  of  motion"  (no  or 
slight  displacement  of  individual  stone  of  a  riprap  top  layer  can  be 
accepted.) 

Conversely,  a  probabilistic  approach  calculates  total  damage  during  a  time  period  based 
on  number  of  stones  transported. 

A  generalized  form  of  Equation  79  was  found  in  PIANC  (1987)  for  the  transversal 
stem  wave.  That  equation  is 
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N  <  __iE^ 

=  Ad 


<1.5  (cot  a) 


1/3 
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(85) 


At  cot  a  =  4  ,  Nj  =  2.38  ,  making  it  similar  to  that  found  above  for  a  4: 1  slope.  They 
stated  that  the  water  level  depression  can  be  considered  a  "truncated  form  of  a  long  period 
wave."  Water  level  depression  causes  excess  pore  water  pressure  because  the  soil  in  the 
banks  cannot  respond  immediately  to  the  change. 

The  secondary  wave  height,  Hj ,  of  the  interference  peaks  at  the  revetment  is 
estimated  using  Equation  25.  The  stable  riprap  size  was  determined  as 

H 

^  1.8(cosP)-“-'  (86) 

Again,  the  angle  of  incidence,  p  ,  is  assumed  to  be  55  deg. 

PIANC  (1987)  also  used  the  Isbash  approach  to  evaluate  return  current  protection. 
The  equation  presented  is  similar  to  that  given  by  Blaauw  et  al.  (1984)  in  Equation  76  but 
uses  different  coefficients  and  the  maximum  return  current,  u^  , 


^50  ^ 


0.7  u 


(87) 


The  coefficient,  k^  ,  is  defined  by  Equation  77  and  the  angle  of  repose,  (j) ,  is  assumed  to 
be  between  40  and  45  deg.  This  reference  stated,  "Current  attack  will  affect  the  entire 
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wetted  zone  of  the  bank.  However,  wave  attack  is  more  likely  to  be  the  dominant  effect 
in  the  zone  near  the  water  line." 

Sorensen  (1986)  recommended  using  the  Hudson  equation  with  a  stability  coefficient 
using  the  incident  wave  height,  H; ,  as  follows: 


=  1.37(cota)^^^ 


(88) 


He  suggested  this  approach  "until  further  insight  is  gained  on  the  impact  of  a  group  of 
vessel  waves  on  a  riprap  armor  layer."  He  further  suggested  that  layer  thickness  and 
gradation  designs  be  obtained  from  Engineering  Manual  1 110-2-1614  (Headquarters, 
U.S.  Army  Corps  of  Engineers  (HQUSACE)  1985). 

Other  studies 

Maynord  (1984)  conducted  a  site  specific  study  of  the  Sacramento  ship  channel 
because  existing  channel  protection  was  experiencing  extensive  channel  maintenance  due 
to  the  extremely  small  blockage  ratios  (typically  around  4).  The  model  section  tested  bad 
a  trapezoidal  channel  with  a  shallow  berm.  The  study  results  found  that  with  small 
increases  in  vessel  speed,  otherwise  stable  rock  failed  due  to  large  increases  in  drawdown, 
which  was  due  to  the  small  blockage  ratio. 

Prototype  gradations  of  W50  of  90  lb  and  340  lb  (40.8  kg  and  154.2  kg)  were  tested 
in  this  model.  Failure  of  the  small  stone  was  observed  for  two  blockage  ratios,  4.3  and 


50 


6.1.  V essel  speeds  for  initiation  of  motion  respectively,  were  approximately  8  mph  and 
8.5  mph  (12.9  km/h  and  13.7  km/h) .  The  drawdown  associated  with  a  6  ft  (1.8  m)  depth 
over  the  test  berm  was  approximately  2  ft  (0.6  m)  for  these  conditions.  An  increase  in 
channel  area  of  20  percent  only  affected  the  maximmn  speed  by  an  increase  of  less  than 
1  mph  (1.6  km/h).  The  larger  gradation  was  only  stable  for  increases  in  ship  speeds  of 
approximately  1  mph  (1.6  km/h). 

A  physical  model  study  was  conducted  for  the  U.S.  Army  Engineer  District,  Mobile 
by  WES  regarding  riprap  stability  on  the  Tennessee-Tombigbee  Waterway.  The  main 
objective  of  this  study  was  to  find  optimum  tow  speeds  for  existing  riprap  (Maynord  and 
Oswalt  1986).  Two  gradations,  Wjq  of  27  and  13  lb  (12.2  and  5.9  kg),  smaller  than  the 
existing  stone,  and  Wjq  of  68  lb  (30.8  kg),  were  tested  in  a  1 :20  model  study  by  driving  a 
firee-moving  vessel  past  the  embankment  under  various  tow  speeds,  water  depths,  and 
sailing  lines.  A  general  equation  was  not  developed,  but  incipient  motion  data  verified 
that  the  Hudson  equation  (Equation  56)  would  be  safe.  The  conclusion  stated  that 
operating  tow  speeds  less  than  80  to  90  percent  of  the  limiting  speed  would  be  sufficient 
to  prevent  failure  of  the  existing  stone  due  to  tow  waves.  However,  it  was  concluded  that 
propeller  jets  directed  toward  the  stone  could  easily  fail  the  stone.  Observations  indicated 
that  incipient  motion  of  the  smaller  stone  did  not  occur  for  wave  heights  less  than  2  ft 
(0.6  m). 

In  a  preliminary  evaluation  of  results  to  date  from  data  collected  for  navigation 
research  and  from  several  site-specific  studies,  including  the  Tennessee-Tombigbee  data 
Martin  (1992)  observed  in  channels  with  blockage  ratios  (generally  from  6  to  10)  that 
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stone  having  a  W50  of  13  lb  (5.9  kg)  was  stable  for  tows  traveling  at  speeds  of 
approximately  6  to  8  mph  (9.7  to  12.9  km/h).  It  was  also  observed  that  the  slope 
protection  on  the  3H:  IV  slope  could  withstand  higher  vessel-induced  forces  than  the 
2H:1V  slope. 

Gelencser’s  (1977)  conclusions  from  studies  of  the  Welland  Canal,  both  prototype 
and  model  investigations,  stated  the  following: 

a.  At  vessel  speeds  of  14  to  15  mph  (22.5  to  24.1  km/h),  a  730-ft  (222.6-m)  vessel 
failed  riprap  having  a  maximum  diameter  of  5  ft  (1.5  m). 

b.  Reducing  the  speed  to  10  mph  (16.1  km/h)  resulted  in  riprap  protection  with 
practically  no  damage  for  stones  having  a  maximum  diameter  of  30  in.  (0.76  m). 

c.  Bank  protection  could  be  accomplished  with  two  stone  sizes  due  to  the 
concentration  of  erosive  forces  on  the  upper  slope  near  the  water  line. 

In  the  prototype  data,  blockage  ratios  for  displacement-type  vessels  ranged  from  5.3  to 
22,  vessel  speed  (relative  to  water)  from  approximately  8.8  to  17  mph  (14.2  to  27.4  km/h) 
with  a  maximum  recorded  drawdown  of  2.78  ft  (0.85  m). 

Balanin  and  Bykov  (1965)  stated  in  their  introduction,  "Selection  of  rational  bank 
protection  is  the  most  complicated  problem  to  be  solved  under  navigable  canal 
conditions..."  The  philosophy  employed  in  design  of  navigation  canals  in  the  previous 
USSR  balanced  costs  of  bank  protection  with  increasing  canal  size.  They  did  not  give 
specific  criteria  for  designing  riprap  protections  but  stated  that  riprap  protection  was  a 
practical  and  economical  solution  for  navigation  canals.  Balanin  and  Bykov  also 
suggested  that  lower  slopes  can  have  lesser  protection  than  upper  slopes. 
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Tenaud  (1977)  conducted  a  number  of  model  tests  specifically  designed  to  determine 
riprap  protection  from  tow  induced  impacts.  Tenaud's  stone  size  recommendations  were 
as  follows: 

a.  2H:  1 V  slope,  66.2  lbs  (30  kg)  at  24  in.  (0.6  m)  thickness  for  the  entire  slope. 

b.  3H:1V  slope,  44.1  lbs  (20  kg)  at  20  in.  (0.5  m)  thickness  for  4.9  ft  (1.5  m)  of 
slope. 

This  "design"  was  based  on  a  ship  with  a  3.2-ft  (0.98-m)  draft  navigating  near  the  limit 
speed  of  8.7  mph  (14  km/h)  in  a  channel  section  with  a  cross-sectional  area  of  2,152  ft^ 
(200  m^).  The  ship  beam  was  approximately  37.4  ft  (1 1 .4  m),  making  the  blockage  ratio 
approximately  17. 

Tenaud  differentiated  between  drawdown  and  the  "attacking  wave"  saying  the 
attacking  waves  (secondary  waves)  can  be  larger  for  lighter  drafts.  He  stated,  "empty  or 
light  laden  boats,  as  well  as  service  boats,  travel  at  high  speed,  and  this  contributes  the 
most  to  bank  erosion."  According  to  Tenaud,  the  attacking  wave  rather  than  the  water 
level  drop  has  the  greatest  erosion  effect  on  the  bank. 

Bouwmeester  et  al.  (1977)  conducted  a  l:25-scale  model  test  to  observe  the  behavior 
of  bank  protection  and  measure  the  primary  and  secondary  waves.  They  concluded  that 
the  transverse  stem  wave,  the  transition  from  the  water  level  depression  alongside  the 
ship  and  the  region  behind  the  ship,  was  likely  the  cause  of  rock  displacements.  Slope 
supply  flow  contributed  to  the  transport  and  the  depressed  water  level  caused  ground 
water  flow.  Stones  were  displaced  no  higher  (vertically)  than  1.6ft  (0.5  m)  and  no  lower 
than  6.6  ft  (2  m)  below  the  still  water  level.  They  stated. 
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Summarizing,  it  can  be  said  that  the  sum  total  of  the  water  motions 
described  above  gives  rise  to  pressure  gradients  in  and  above  the  banks, 
and  erosive  currents,  which  may  result  in: 

a.  loss  of  stability  inside  the  slope  because  resistance  to  sliding  is  reduced; 

b.  loss  of  stability  of  individual  particles  in  and  near  the  surface  of  the  bank. 


Pilarczyk's  (1984)  paper  supported  general  conclusions  of  research  conducted  by 
both  the  Delft  and  WES.  The  prototype  tests  on  the  Hartel  Canal  were  conducted  with 
gravel  on  geotextiles  and  sandy  subsoil.  Riprap  sizes  on  one  embankment  were  1 .2-  to 
3.2-in.  (30-  to  80-mm)  gravel  and  3.2-  to  7.9-in.  (80-  to  200-mm)  gravel  on  another.  The 
canal  had  a  246-ft  (75-  m)  bottom  width  and  was  approximately  23  ft  (7  m)  deep.  The 
banks  had  4H:  1 V  slopes.  By  comparison  the  Tennessee-Tombigbee  had  a  280-ft 
(85.4-m)  bottom  width  and  was  14  to  22  ft  (4.3  to  6.7  m)  deep.  Boat  speeds  were  not 
given,  but  the  push  tows  had  4,500  hp  and  5,400  hp.  Each  barge  was  250  ft  long  by  39  ft 
wide  by  9.8  ft  draft  (76.5  by  1 1.8  by  3  m).  The  configurations  were  3  by  2  and  2  by  3. 
Blockage  ratios  for  these  tests  were  generally  less  than  10. 

The  studies  concluded,  "High  gravel  transports  were  only  observed  when  push-tows 
and  tugs  sailed  at  high  speeds  near  the  bank."  Also,  small  vessels  near  the  bank  appeared 
to  cause  a  high  transport  of  fine  gravel.  For  riprap  1 1  to  88  lbs  (5  to  40  kg)  on  geotextile, 
Pilarcyzk  (1984)  observed: 

For  this  class  riprap,  the  beginning  of  movement  was  rather  exceptional.  It  has 
been  observed  to  occur  mainly  due  to  secondary  waves  induced  by  empty  barges 
and  tugs  sailing  very  close  to  the  bank  (at  the  toe  of  the  slope).  It  can  be 
concluded  that  in  the  most  practical  cases  this,  or  a  little  higher  class  of  riprap 
would  be  satisfactory  for  the  normal  inland  fairways.  In  the  Netherlands  the 
stone  classes  used  normally  for  bank  protection  of  chaimels  with  high  shipping- 
intensity  (inch  push-tow)  are  of  the  10-60  kg  and/or  60-300  kg  types. 


54 


Variability  in  riprap  size  resulting  from  equations 


The  wide  array  of  form  and  parameters  in  the  coefficient  produces  a  design  rock 
weight  that  can  vary  considerably.  Table  1  summarizes  typical  values  of  Nj  developed 
for  waves.  These  formulas  represent  different  interpretations  of  the  force  balance, 
stability  criteria,  assumptions  regarding  nominal  diameter  (cubic  or  spherical  stones), 
design  wave  height,  and  waves  characteristics  such  as  the  surf  parameter  and  angle  of 
attack.  Figure  9  shows  the  relationship  between  three  different  wave  heights  and  the 
nominal  stone  diameter  using  the  definition  of  stability  number  in  Equation  51.  The 
specific  weight  of  the  stone  is  Yr  165  Ib/ft^  (2.12  kg/m^). 

The  stability  coefficients  determined  from  these  formulations  resulted  in  highly 
variable  stable  rock  designs.  Even  among  the. formulas  developed  for  vessel  waves,  the 
range  of  is  from  1.32  from  the  Fuehrer,  Romisch,  and  Engelke  (1981)  formulation  to 
2.38  in  the  PIANC  (1987)  criteria  for  a  2H:1  V  slope.  Going  to  Figure  9,  the  resulting 
design  stone  diameter  for  a  2-ft  (0.6-m)  wave  height  can  range  from  approximately  1 1  to 
6  in.  (0.3  to  0.15  m).  Then  resulting  stone  weights  are  from  129  to  22  lb  (58.5  to  10  kg), 
respectively.  At  a  wave  height  of  3  ft  (0.9  m),  the  stone  weights  vary  from  436  to  74  lb 
(197.7  to  33.6  kg). 

The  existing  stability  equations  produced  a  wide  variation  in  stable  riprap  size,  and 
with  the  exception  of  the  Teimessee-Tombigbee  study  (Maynord  and  Oswalt  1986)  and 
the  Hartel  Canal  study  (Pilarczyk  1984),  few  have  been  specifically  verified  at  full-scale 
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Table  1 

Values  of  the  Stability  Number,  N^,  Based  on  Different  Formulations 


Source"' 


Equation  cot  a  cot  a 

#  =2  =3 


Wind  Wave  Equations 


Iribarren  (1948) 


Hudson  (1957) 


Hedar2(1960) 


Bhowmik  (1976) 


Ahrens  (1989) 


Van  der  Meer  and  Pilarczyk  (1984) 


PIANC  (1987) 


Vessel  Equations 


Maximum  Front  Wave 


Fuehrer,  Romisch,  and  Engelke  (1981) 


Transversal  Stern  Wave 


Blaauw  et.  al.  (1984) 


PIANC  (1987) 


Secondary  Waves 


Blaauw  et.  al.  (1984) 


Verhey  and  Bogaerts  (1989) 


PIANC  (1987) 


Sorensen  (1986) 


1.32  1.51  1.67 


1.89  2.16  2.38 


'  To  obtain  values  of  N^  the  following  assumptions  were  used  depending  on  the 
formula:  (j)  =  45° ;  p  =  55°;  ^  =  3;  N  =  1000;  S  =  2. 

^down  was  obtained  for  each  slope  from  a  graph  in  Hedar  (1960). 


2 


Figure  9.  Relationship  between  and  D50 


for  commercial  tow  traffic.  Determination  of  a  generalized  equation  would  minimize  the 
need  for  site-specific  studies  and  maximize  the  safe  and  economical  design  of  riprap  in 
United  States  inland  navigable  waterways. 

Scale  Effects  in  Physical  Modeling 

Physical  model  scales  should  be  carefully  selected  to  preserve  similarity  between  the 
model  and  the  prototype.  Hydraulic  models  usually  are  geometrically  similar  (have  the 
same  length  ratios  from  prototype  to  model).  To  ensure  dynamic  similarity,  either  the 
Froude  number  or  the  Reynolds  number  must  be  the  same  in  model  and  prototype.  The 
Froude  number  approach  accurately  scales  gravitational  effects  and  the  Reynolds  niunber 
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accounts  for  viscous  effects.  Since  no  fluid  exists  to  satisfy  both  criteria  simultaneously, 
Froude  criteria  is  generally  used  and  often  an  empirical  adjustment  is  made  to  the  forces 
affected  by  viscosity  (Sharp  1981).  A  model  scale  is  selected  such  that  Reynolds 
numbers  are  sufficiently  high  to  minimize  these  effects  (Valin  1989). 

In  the  study  of  ship  dynamics,  wave  dynamics,  stone  slope  protection,  and  navigation 
effects,  researchers  have  employed  many  different  physical  model  scales.  The  selection 
of  scale  is  based  on  a  combination  of  Reynolds  numbers  and  on  experience.  When  scale 
effects  can  affect  the  measured  outcome  of  the  parameter  of  interest,  the  method  for 
compensation  of  the  effects  is  a  function  of  the  researcher  and  the  techniques  used  by  the 
researcher’s  laboratory. 

Model  scale  and  Reynolds  values 

Vessel  models.  Ship  modeling  presents  some  scaling  challenges  related  to  boundary 
conditions.  As  a  vessel  moves  through  the  water,  it  develops  a  boundary  layer  on  the 
bottom  of  the  channel  and  the  bottom  of  the  vessel  (Gates  and  Herbich  1977).  Because 
the  boundary  layer  is  a  function  of  Reynolds  number  and  viscosity,  this  value  can  be 
different  in  the  model  and  prototype.  The  boundary  layer  problem  is  compoimded  in  the 
model  due  to  low  Reynolds  numbers  which  causes  a  greater  influence  due  to 
displacement  thickness  than  actually  occurs  in  the  prototype.  In  other  words,  the  smaller 
the  model  scale,  the  greater  the  viscous  forces  and  the  displacement  thickness; 
consequently,  the  effective  size  of  the  tow  is  larger  than  reality.  The  shape  of  the  lead 
barges  can  cause  flow  separation  both  in  the  horizontal  and  vertical  planes.  In  the  vertical 
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plane  this  further  complicates  the  issue  of  boundary  layer  development  beneath  the  tow 
and  influences  its  effective  draft.  In  the  horizontal  plane,  this  separation  pushes  the  flow 
field  away  from  the  sides  of  the  barges,  effectively  increasing  its  width.  This  violates  the 
assumptions  in  the  theoretical  approach  presented  by  Jansen  and  Schijf  (1953). 

Sorensen  (1966)  foxmd  that  viscous  effects  can  cause  damping  of  very  short  waves. 
He  conducted  a  series  of  scale  effect  tests  to  determine  the  circumstances  affected  by 
viscous  damping  for  ship  studies  in  deep  water.  He  concluded  that  some  scale  effects 
were  present  at  low  Froude  numbers  (less  than  0.5)  and  with  vessel  speeds  of  less  than 
1.8  fps  (0.5  mps).  He  recommended  avoiding  the  smallest  model  tested. 

In  the  computation  of  drag  coefficients  in  towing  tank  tests.  Sharp  (1981) 
recommends  using  velocities  based  on  Freudian  criteria  and  adjusting  the  drag  forces  by 
an  empirical  calculation  of  the  viscous  forces.  He  further  stated  the  model  scale  can  be 
increased  to  make  the  differences  between  prototype  and  model  negligible. 

At  WES,  model  scales  from  as  large  as  1 : 1 0  to  as  small  as  1:120  have  been  used  to 
evaluate  navigation  effects.  The  tests  at  1 :70  and  smaller  were  influenced  by  the  scale 
effects  and  the  studies  provided  qualitative  results  at  best.  A  1 :37.5  scale  model  has  been 
used  on  several  occasions  due  to  lack  of  adequate  space  to  evaluate  return  currents  in 
large  river  sections,  but  it  was  recognized  that  the  absolute  magnitude  of  the  answers 
require  adjustment.  Maynord  (1984)  has  also  modeled  ships  and  tows  at  1 :10, 1 :20, 
and  1:30. 

The  boundary  layer  and  separation  problems  can  cause  the  effective  size  of  the 
vessel  to  change  in  the  model.  In  a  study  conducted  at  a  1 :37.5  scale  by  WES  for  the 
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U.S.  Army  Engineer  District,  Louisville,  a  comparison  was  made  between  model  data  and 
prototype  data.  It  was  found  that  for  the  channel  modeled,  the  effective  draft  of  the  tow 
was  actually  on  the  order  of  10.8  ft  not  9  ft  (3.3  m  not  2.74  m).  Also  in  a  recent 
1 :25-scale  model  study  for  the  Upper  Mississippi  River  System  where  model  values  were 
compared  to  actual  field  data,  the  model  draft  was  lightened  from  a  protot5^e  draft  of  9  ft 
(2.74  m)  to  approximately  7.5  ft  (2.29  m)  to  reproduce  field  data  (Maynord  and  Martin 
1996). 

Based  on  Reynolds  numbers  related  to  vessel  length,  the  rudder,  and  the  propeller, 
Novak  and  Cabelka  (1981)  recommended  that  self-propelled  vessels  be  reduced  at  most 
to  a  1 :15  scale.  Bouwmeester  et  al.  (1977)  conducted  their  experiments  with  a  free 
moving  vessel  at  a  1:25  scale.  Verhey  (1983)  conducted  studies  primarily  regarding 
effects  of  propeller  jets  on  channel  bottom  and  banks  and  assumed  scale  effects  were 
insignificant  at  a  1 :25  scale. 

Koster  (1975)  did  a  scale  effects  study  using  1:10.5, 1:25,  and  1:40;  ultimately 
recommending  the  use  of  a  1 :25-scale  model.  In  his  tests,  the  scale  effects  to  the  rudder 
action  and  drift  angle  seemed  to  counteract  each  other,  making  it  difficult  to  isolate  and 
compare  them  to  prototype  values. 

Scale  effects  due  to  frictional  resistance  are  sometimes  compensated  by  increasing  the 
propeller  speed.  At  the  Netherlands  Ship  Model  Basin,  a  10  percent  increase  in  propeller 
speed  was  used  over  the  Froude  scaling  criteria  (Kooman  1973).  These  studies  were 
conducted  at  a  1 :25  scale.  Instead  of  adjusting  propeller  speed,  tests  conducted  at  WES 
to  evaluate  propeller  effects  are  related  to  the  pushing  force,  or  bollard  push.  Force 
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measurements  were  made  on  the  model  tow  at  various  propeller  speeds  and  scaled  to 
prototype  values  according  to  Froude  criteria  (Toutant  1982). 

Riprap  and  wave  models.  Hudson  (1975)  suggested  the  model  Reynolds  number, 
Rq  ,  be  greater  than  3  x  in'*  for  conducting  riprap  stability  tests  under  wave  activity. 
Reynolds  is  defined  in  terms  of  the  model  wave  height,  H,  and  stone  size,  according  to 
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where 

V  =  kinematic  viscosity  of  water,  ft^/s  (mVs) 

Parameters  in  Hudson's  (1957)  small-scale  tests  to  evaluate  non-breaking  waves  on 
pervious  breakwaters  that  resulted  in  Equation  56  were  as  follows:  wave  heights  from 
0.28  to  0.69  ft  (0.09  to  0.21  m)  with  periods  of  0.88  to  2.65  sec;  water  depth  from  1.26  to 
2.00  ft  (0.38  to  0.6  m);  quarry  stone  sizes  from  0.09  to  0.31  lb  (0.04  to  0.14  kg);  and 
seven  breakwater  slopes  from  cot  a  of  1.25  to  5. 

O'Loughlin  et  al.  (1970)  discussed  scale  effects  on  rock  due  to  jets.  They  stated  that 
scour  rate  is  "practically  independent  of  Reynolds  number  at  values  above  2.5  x  10^..." 
The  authors  define  the  particle  Reynolds  number,  Rp ,  as: 
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where 

V  =  jet  efflux  velocity,  (mp/s) 

D  =  particle  diameter,  ft  (m) 

In  their  experiments,  data  corresponding  to  small  beads  (3  mm)  fell  out  of  sync  with  data 
from  larger  diameter  beads  (6  mm  to  15  mm)  indicating  a  scale  effects  problem. 

Verhey  (1983)  in  his  studies  of  the  effects  of  propeller  jets  on  the  channel  bottom  and 
banks,  also  recommends  using  the  particle  Reynolds  number  and  that  the  model  have 
values  greater  than  3,000. 

Hedar  (1960)  conducted  his  experiments  in  a  wave  tank  at  a  water  depth  of  1 .8  ft 
(0.5  m)  with  crushed  stone  on  an  impermeable  sheet  of  pl5^ood.  Bank  slopes  varied 
from  1.125H:1V  to  4.5H:1V.  The  median  diameter  ofthe  stone  was  1.87  in.  (47.5  mm). 
Wave  period  was  2  sec  and  wave  heights  were  from  0.25  to  0.75  ft  (80  to  230  mm).  The 
scale  was  1:10  and  1 :20.  There  was  no  discussion  of  scale  effects  or  Reynolds  criteria. 

Dai  and  Kamel  (1969)  performed  tests  to  evaluate  the  stability  of  rough  quarry  stone 
from  nonbreaking  waves  on  a  pervious  rubblemound  breakwater.  Two  small  rock  sizes 
were  tested  and  compared  to  a  prototype  size  stone.  According  to  them,  the  fiiction  force 
is  directly  proportional  to  the  product  of  the  coefficient  of  fiiction  and  the  weight  of  the 
armor  unit.  They  found  appreciable  scale  effects  in  tests  conducted  at  the  smallest  scale. 
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Dai  and  Kamel  attribute  these  effects  to  the  relatively  high  value  of  the  drag  coefficient 
for  the  small  model  versus  the  large  model.  The  higher  the  coefficient  the  less  stable  the 
cover  layers.  Viscous  forces  were  predominant  and  the  drag  coefficient  decreased  with 
increasing  values  of  the  Reynolds  number  until  it  reached  a  critical  value  of  3  x  10"^.  The 
Reynolds  number  was  given  as 

VrD 

Rd  =  (91) 

V 

where 

Vr  =  value  of  water  particle  velocity  parallel  to  breakwater  side  slope  at  a  distance  R 
equal  D/2. 

The  smallest  stone  had  a  weight  of  0.048  lb  (0.02  kg)  with  a  diameter  of  0.20  in. 

(5  mm).  Water  velocities  in  the  wave  near  the  breakwater  slope,  Vr  ,  were  fi'om  1.58  to 
2.8  fys  (0.48  to  0.85  mps)  and  zero  damage  wave  heights  were  0.14  to  0.22  ft  (0.04  to 
0.07  m).  The  next  size  rock  had  a  weight  of  0.38  lb  (0.17  kg)  and  a  diameter  of  0.57  in. 
(14.5  nun).  Zero  damage  wave  heights  were  from  0.35  to  0.44  ft  (0.1 1  to  0.13  m)  and 
velocities  were  from  2.6  to  3.0  fps  (0.8  to  0.9  mps). 

Madsen  and  White  (1975)  conducted  tests  for  reflection  coefficients  on  rough 
impermeable  slopes.  The  theoretical  assumptions  are  (a)  relatively  long,  normally 
incident  waves  (linear  long  wave  theory)  and  (b)  energy  dissipation  on  a  rough 
impermeable  slope  is  a  function  of  bottom  fiictional  effects.  Their  experiments  were 
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considered  valid  for  folly  rough  turbulent  flow  as  defined  by  the  following  wave 
Reynolds  number,  ,  relationship  developed  by  Jonsson  (1966): 


where 
A 
6) 

The  conditions  tested  by  Madsen  and  White  (1975)  on  an  impermeable  slope  were: 
wave  periods  of  2,  1.8,  and  1.6  sec;  a  water  depth  of  1  ft  (0.3  m);  djo  of  0.5  in.  (12.7  mm), 
1.0  in.  (25.4  mm),  1.5  in.  (38.1  mm),  and  2.0  in.  (50.8  mm);  and  wave  heights  varied 
from  0.05  to  0.17  ft  (0.02  to  0.05  m). 

Broderick  conducted  small-scale  tests  in  the  laboratory  and  compared  them  to  large 
wave  tank  tests  conducted  by  Ahrens,  essentially  representing  foil-scale  field  conditions 
(Broderick  and  Ahrens  1982).  Failure  of  stone  was  defined  as  "the  riprap  being  shifted 
enough  to  expose  part  of  the  filter  layer  which  was  removed  by  the  wave  action."  The 
authors  measured  volume  of  erosion  by  profiling  the  slope  between  tests. 

Broderick  and  Ahrens  (1982)  also  foimd  evidence  that  at  higher  Reynolds  values, 
scale  effects  diminished.  They  calculated  Reynolds  values  using  two  different  equations. 
The  first  was  a  variation  of  Equation  92  used  by  Madsen  and  White  (1975).  Fully 
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representative  wave  amplitude  or  the  runup 
radian  frequency  (27i:/T) 
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turbulent  flow  was  defined  by  Rg  as  a  function  of  wave  period,  zero-damage  runup,  and 
slope  angle  such  that 


(93) 


where 

=  the  wave  runup  associated  with  the  zero-damage  wave  height,  ft  (m) 

T  =  the  wave  period,  sec 
a  =  the  slope  angle 

The  critical  values  of  Rg  were  defined  by  a  graph,  where  Rg  was  given  as  a  function  of  a 
roughness  term.  According  to  this  criteria,  the  authors  stated  that  both  small-scale  and 
prototype  tests  were  in  the  rough  turbulent  flow  regime. 

A  Reynolds  number,  R^  was  also  computed  using  Equation  89,  but  replacing  H 
with  the  zero  damage  wave  height,  Hj. ,  and  ’withWjQ. 
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According  to  a  figure  given  in  their  paper,  scale  effects  were  less  evident  at  an 
greater  than  10^.  The  authors  evaluated  results  from  other  sources  and  foimd  that  where 
low  values  of  were  observed,  was  as  much  as  40  percent  lower  than  prototype 
values.  Whereas,  in  the  tests  Broderick  conducted,  the  stability  numbers  of  the  small 
scale  tests  were  observed  to  be  only  about  20  percent  lower  than  full-scale  values. 

Broderick  conducted  her  studies  at  a  1 ;  10  scale.  The  median  weight  of  the  model 
riprap  was  0.075  lb  (0.03  kg)  on  a  bank  slope  of  3.5H:  1 V.  The  cover  layer  had  a  specific 
gravity  of  2.71  and  was  laid  on  a  granular  filter  on  top  of  a  compacted  sand  core.  The 
still  water  depth  in  the  flume  was  1.5ft  (0.46  m),  and  the  wave  heights  varied  from  2  to 
3  in.  (50.8  to  76.2  mm).  Number  of  waves  varied  from  340  for  long  period  waves  to 
1,050  for  short  period  waves.  Wave  periods  in  the  model  were  from  about  0.9  to  3.6  sec. 

In  conclusion  this  reference  stated,  "Tests  at  1:10  Froude  scale  yield  zero-damage 
stability  numbers  about  20  percent  lower  than  prototype  tests...  Scale  effects  were  less 
severe  at  high  levels  of  damage  than  at  zero-damage  level."  Since  both  small-scale  and 
large-scale  tests  were  considered  to  be  in  the  rough  turbulent  regime,  the  20  percent 
difference  is  still  imexplained!  Other  conclusions  are  as  follows:  runup  at  zero  damage  is 
20  percent  higher  in  the  small  scale  tests;  damage  profiles  are  the  same  in  model  and 
prototype;  and  wave  period  had  little  influence  in  the  small-scale  tests. 

Van  der  Meer  and  Pilarczyk  (1984)  used  broken  stone  with  the  following 
characteristics:  W50  =  0.27  lb  (0.12  kg);  D50  =1.42  in.  (36.1  mm);  layer  thickness  = 

3.15  in.  (80  mm);  and  filter  layer  thickness  =  0.79  in.  (20.1  mm).  The  slope  was 
constructed  of  mortar  when  an  impermeable  core  was  being  tested.  A  total  of  3,000 
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waves  ranging  from  a  height  of  0.16  to  0.85  ft  (0.05  to  0.26  m),  periods  of  1.3  to  3.2  sec, 
a  water  depth  of  2.62  ft  (0.8  m),  and  slope  angles  from  cot  a  of  2  to  6  were  tested.  Scale 
effects  were  not  mentioned. 

Breakwater  studies  have  been  conducted  at  variable  scales  as  well.  Jensen  (1989) 
suggested  that  typical  breakwater  studies  are  conducted  at  a  1:30  scale.  To  maintflin 
similarity  in  the  cover  layer  and  the  core,  the  core  of  the  structure  should  be  built  of  larger 
stones  than  required  by  geometric  similitude.  Jensen  (1989)  recommended  a  lower  value 
of  the  critical  particle  Reynolds  number  (6,000  to  7,000)  for  the  cover  layer  than  Dai  and 
Kamel  (1969).  Jensen  calculated  the  Reynolds  number  using  Equation  90  and  fmmd  a 
high  degree  of  variability  in  the  minimum  or  critical  Reynolds  criteria. 

Oumeraci  (1984)  reevaluated  test  results  conducted  by  Hudson  (1957),  Broderick  and 
Ahrens  (1982),  Dai  and  Kamel  (1969),  and  others  to  determine  the  critical  Reynolds 
values  to  minimize  scale  effects.  Oumeraci  found  a  wide  variation  in  reported  critical 
values  of  Reynolds  with  Dai  and  Kamel's  representation  at  a  value  of  3  x  10'*  being  the 
most  reasonable.  He  plotted  the  stability  number  against  the  Reynolds  value  and 
suggested  a  method  for  adjusting  the  stability  numbers  in  the  model  if  the  scale  was  not 
sufficient  to  meet  the  Reynolds  criteria.  Oumeraci  also  recommended  a  mintmum  wave 
height  in  the  model  of  0.33  to  0.49  ft  (0.1  to  0.15  m)  and  unit  weights  of  stone  of  0.22  to 
0.33  lb  (0.1  to  0.15  kg).  Like  Jensen,  Oumeraci  recommended  that  core  materials  be 
geometrically  larger  in  the  model  to  minimize  pressures  inside  the  structure. 

Koyabashi  and  Greenwald  (1988)  conducted  wave  tests  in  a  flume  on  an  impermeable 
3H:  1 V  slope  for  calibrating  a  numerical  model.  Using  one  gradation  of  model  gravel 
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with  a  median  diameter  of  0.83  in.  (21  mm),  they  measured  values  of  the  specific  weight, 
friction  angle,  and  volume  coefficient  of  stone  for  input  to  their  numerical  model 
equations.  Results  were  not  compared  to  another  larger  model  gradation  or  to  the 
prototype.  Experimental  error  on  the  order  of  10  percent  was  discussed,  but  error  due  to 
scale  effects  was  not. 

Riprap  and  vessel  models.  In  studies  where  both  navigation  effects  and  riprap 
stability  were  tested,  there  have  also  been  several  size  scales  used.  Tenaud  (1977)  used  a 
scale  of  1 :20  in  a  flume  262.5  ft  (80  m)  long.  Blockage  ratios  varied  from  5.6  to  33.7. 
Model  barges  were  10.25  ft  (3.1  m)  long  by  1.87  ft  (0.6  m)  wide  and  drafted  from  0.08  to 
0.49  ft  (0.02  to  0.15  m).  Model  grains  sizes  ranged  from  0.16  to  0.87  in.  (4  to  22  mm)  for 
riprap  and  0.2  to  2  mm  sand.  Slope  angles  included  3H:  IV  and  2H:  IV.  Scale  effects 
were  not  discussed  for  this  model  but  were  determined  to  be  significant  for  a  1 :50-scale 
model. 

Gelencser  (1977)  conducted  tests  at  a  1:60  and  1:80  scale  to  determine  drawdown  in  a 
restricted  ship  canal.  The  blockage  ratio  for  these  studies  was  5  to  7.  The  maximum 
prototype  vessel  speed  of  14  to  15  mph  (22.5  to  24.1  km/h)  was  assumed  to  be  the 
limiting  speed  due  to  the  vessel  touching  the  bottom.  Bank  slopes  were  2.75H:1V  and 
2.25H:1V.  No  details  were  given  regarding  the  materials  used  in  the  bank  stability 
studies. 

Fuehrer,  Romisch,  and  Engelke  (1981)  conducted  extensive  tests  regarding  bank  and 
bottom  protection  against  vessel  waves  and  propeller  jets.  These  tests  were  conducted  at 
low  Reynolds  numbers  and  with  small  grain  sizes.  Their  scale  model  testing  found  that 
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the  limiting  stability  coefficient  was  a  function  of  the  Reynolds  number,  Rq  .  They 
defined  the  Reynolds  number  as  a  function  of  the  front  wave  height  (Equation  15)  by 


^50 
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Fuehrer,  Romisch,  and  Engelke  plotted  the  stability  number  versus  Reynolds  number  and 

concluded  that  at  initiation  of  motion  was  constant  above  a  value  of  1,000.  They 

described  the  Shield's  type  graph  stating. 

The  coefficients  of  the  limiting  stability,  determined  for  the  different  grain  sizes 
and  loading  situations  of  the  canal  slopes  by  the  ship  in  the  case  of  normal 
navigation  can  be  fitted  by  a  reference  curve,  of  which  the  functional  course 
corresponds  to  the  Shield  limiting  tractive  force  relationship  for  the  beginning  of 
bedload  movement  in  flows. 


Fuehrer,  Romisch,  and  Engelke  (1981)  used  a  1 :60-scale  model  with  model  grain 
sizes  from  0.019  to  0.106  in.  (0.48  to  2.7  mm).  Water  depths  were  from  0.69  to  1.08  ft 
(0.21  to  0.33  m)  with  draft  to  depth  ratios  from  1.45  to  1.98.  Blockage  ratios  were  from 
approximately  22  to  28  on  3H:1V  slopes. 

Maynord  and  Oswalt  (1986)  used  a  l:20-scale  model  for  the  riprap  stability  studies  on 
the  Tennessee-Tombigbee  Waterway.  Blockage  ratios  varied  from  4.6  to  16.4  in  a  model 
trapezoidal  channel  with  a  2H:1V  side  slope.  Median  diameter  stone  sizes  tested  ranged 
from  0.32  to  0.56  in.  (8.1  to  14.2  mm).  Water  depths  varied  from  0.7  to  1.1  ft  (0.21  to 
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0.34  m).  A  free-moving  model  tow  and  barges  at  two  drafts  and  several  configurations 
were  used  during  testing.  Model  Re5molds  values  were  computed  to  be  on  the  order  of 
10'^  and  determined  to  be  in  the  fully  rough  turbulent  regimen  according  to  the  criteria  by 
Madsen  and  White. 

Angle  of  Repose 

In  the  literature  discussions  of  scale  effects,  little  information  is  given  regarding  how 
the  natural  angle  of  repose  affects  model  test  results.  While  the  angle  of  repose  is  built 
into  all  the  stability  equations  presented  in  this  paper,  except  the  Hudson  formula,  this 
value  is  generally  assumed  to  be  constant  over  the  range  of  stone  materials  tested  in  the 
model  and/or  used  in  the  prototype.  The  natural  angle  of  repose  is  most  likely  not  a 
constant  over  the  range  of  sizes  tested  in  both  model  and  prototype;  therefore,  a  potential 
source  of  error  exists  in  scaling  up  results. 

Ulrich  (1987)  contended  that  there  are  two  different  angles  of  repose,  one  for 
resistance  to  sliding  and  one  for  resistance  to  wrenching.  The  resistance  to  sliding  is  met 
if  the  ratio  of  parallel  forces  to  normal  forces  is  less  than  the  tangent  of  the  support 
layer’s  angle  of  friction,  as  presented  in  Equation  34.  Ulrich  slightly  modified  the 
balance  of  forces  assuming,  for  downslope,  the  hydrostatic  pressure  is  normal  to  the 
slope.  The  resulting  modified  coefficient  for  "steady,  uniform,  fiilly-turbulent, 
downslope  flow  in  the  absence  of  seepage  and  wave  forces,"  is 
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He  also  stated  that  the  material’s  natural  angle  of  repose,  based  on  dumping  the  stone, 
was  not  the  appropriate  angle  for  the  angle  of  friction.  Ulrich  recommended  a  method  for 
determining  the  value  by  raising  a  layer  of  materials  on  a  tilting  surface  and  observing 
motion.  From  Ulrich’s  experiments,  he  suggested  using  a  constant  bearing  angle  of 
75  deg  making  the  effects  of  this  angle  negligible  in  the  force  equations.  However,  this 
value  was  determined  experimentally  over  a  very  short  length  of  slope,  L ,  (two  stone 
diameters),  even  though  he  determined  that  the  angle  decreases  as  the  L/D50  increases. 

Using  a  similar  experimental  procedure  as  Ulrich,  Maynord  (1988)  determined  the 
angle  of  repose  was  53  deg.  He  concluded  that  the  angle  of  repose  was  not  always  equal 
to  the  material’s  bulk  angle  of  repose  and  that  the  angle  was  affected  by  the  height  of  the 
revetment,  thickness,  and  placement.  As  L/D50  increases,  the  angle  decreases. 

Maynord  deals  with  the  variability  of  the  angle  empirically  as  Hudson  did  in  his 
experiments. 

Hudson  (1957)  abandoned  correlations  using  the  angle  of  repose  in  his  stability 
formula  due  to  the  wide  variation  of  experimental  results  when  attempting  to  measure  this 
value.  From  Hudson's  tests  to  determine  the  friction  coefficient,  defined  as  tan  (j) ,  and 
from  tabular  data  in  his  report,  the  average  value  of  the  friction  coefficient  of  stone 
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dumped  in  water  and  stone  stacked  in  water  for  0.1,  0.3,  and  0.62  lb  (0.05,  0.14,  and 
0.28  kg)  quarry  stone  was  1.06, 1.09,  and  1.20,  respectively. 

Froelich  and  Ogden  (1994)  concluded  from  field  studies  of  dumped  stone,  that  angle 
of  repose  was  more  correlated  to  angularity  and  gradation  rather  than  median  stone  size. 
They  qualify  their  conclusion  with  the  following:  "A  discernible  relation  between  (j)jj  and 
D50  is  not  apparent...  However,  more  data  with  a  larger  range  of  particle  sizes  might  be 
needed  to  detect  such  a  trend." 

Based  on  other  researcher's  testing,  however,  there  is  clearly  an  increase  in  the  angle 
of  repose  with  increasing  grain  size,  and  even  in  Hudson's  testing,  the  average 
experimental  values  increased  with  rock  size.  In  data  provided  by  Froehlich  and 
Ogden  (1994),  the  measured  value  of  angle  of  repose  for  small  stones  in  the  less-than- 
50-mm  range  varied  from  approximately  35.5  to  40  deg. 

Several  references  were  consulted  to  determine  the  angle  of  repose  values  for 
different  stone  sizes  as  compared  to  Hudson's  tests.  Figure  10  shows  a  composite  of  this 
research.  Data  for  this  graph  were  taken  from  Chow  (1959),  ASCE  (1977),  Henderson 
(1966),  Hedar  (1960),  and  Maynord  (1988).  In  Chow  (1959),  points  were  obtained  from 
a  graph  developed  from  Bureau  of  Reclamation  data  for  very  angular  stones.  Values  by 
White  and  Chepil  for  sand  were  obtained  from  the  ASCE  manual.  The  value  from  Hedar 
for  (j) ,  obtained  experimentally,  was  48  deg.  Since  the  reference  does  not  clearly  state  the 
stone  size  he  used  for  the  tests,  an  average  stone  diameter  was  assumed  based  on  the 
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Figure  10.  Comparison  of  angle  of  repose  from  various  sources 


range  of  model  stones  he  tested  (125  to  175  grams).  The  data  labeled  "Simons"  came 
from  a  graph  in  Henderson's  Open  Channel  Hydraulics  (1966)  and  was  after  research  by 
Simons  and  Albertson. 

In  the  curves  developed  by  the  Bureau  of  Reclamation  and  Simons,  it  can  be  seen  that 
at  diameters  greater  than  approximately  0.15  ft  (0.05  m)  the  angle  of  repose  has  an  almost 
constant  value.  Certainly  all  prototype  stones  would  fall  in  this  range  on  the  graph  and 
some  model  stone.  However,  for  small  diameter  stones,  particularly  less  than 
approximately  0.1  ft  (30.5  mm),  the  angle  of  repose  has  a  very  nonlinear  relationship  to 
diameter. 
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Iribarren's  original  formula  for  the  weight  of  armor  units  in  rubblemound  breakwaters 
that  retained  the  coefficient  of  friction  as  a  variable  was  given  in  Equation  54.  If  we  let 
represent  terms  related  to  the  internal  friction,  such  that 


(tan  (j)) 

(tan({)cosa  -  sina) 


(97) 


then,  stone  weight  is  a  function  of  k^  cubed,  and  stone  diameter,  D ,  is  directly  related 
to  k^  as 


D 


_  1/3  <P 

A 


(98) 


Using  Hudson's  experimental  values  for  (j> ,  Iribarren's  formulation  of  k^  is  calculated  in 
Table  2. 


Table  2 

Calculation  of  Using  Hudson's  Experimental  Data 

W„  lb 

4) 

tan  4) 

a 

tan  a 

mm 

0.10 

46.7 

1.06 

26.6 

2 

2.12 

0.30 

47.5 

1.09 

26.6 

2 

2.07 

0.62 

50.2 

1.20 

26.6 

2 

1.92 

0.10 

46.7 

1.06 

18.4 

3 

1.58 

0.30 

47.5 

1.09 

18.4 

3 

1.52 

0.62 

50.2 

1.20 

18.4  . 

3 

1.46 

0.10 

46.7 

1.06 

14.0 

4 

1.35 

0.30 

47.5 

1.09 

14.0 

4 

1.34 

0.62 

50.2 

1.20 

14.0 

4 

1.30 
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Since  the  stable  diameter,  D  ,  is  a  direct  function  of  ,  it  can  be  seen  that  for 
smaller  angles  of  repose,  increases;  therefore,  a  larger  stone  size  is  required.  As  the 
angle  of  repose  approaches  the  slope  angle,  k^  is  more  influenced  by  the  angle  of  repose. 
According  to  Equations  55  and  98,  is  inversely  proportional  to  k^ ,  so  as  k^ 
increases,  Nj  decreases. 

Another  example  using  the  shear  stress  approach.  Equation  76  and  77,  requires  the 
factor,  kj ,  that  is  also  dependent  upon  the  angle  of  repose,  (j) .  The  relationship  is  such 
that  as  4)  gets  smaller  so  does  kj .  Furthermore,  the  required  stable  stone  size  is 
inversely  proportional  to  k^  resulting  in  larger  stone  sizes  for  smaller  coefficients. 

Ignoring  variations  in  the  angle  of  repose  could  accoimt  for  more  observed  failures  in 
the  model  at  small  stone  sizes  than  previously  observed  in  the  prototype  for  equivalent 
riprap  sizes.  In  several  papers  where  small  model  stones  were  used  in  testing,  the  amount 
of  stone  movement  and  degree  of  failure  were  in  disagreement  with  observations  for 
equivalent  larger  stone.  Recall  in  Broderick  and  Ahrens  (1982)  that  even  though  the  flow 
regime  for  the  smaller  stone  was  considered  to  be  fully  turbulent,  stability  numbers  for 
the  small  stone  were  as  much  as  20  percent  lower  than  the  equivalent  prototype  values. 
They  had  no  explanation  for  this  difference. 

Van  Gent  (1995)  used  small-scale  modeling  to  calibrate  a  numerical  approach  for 
predicting  geomorphic  reshaping  of  permeable  breakwaters  due  to  wave  action.  He  first 
adjusted  the  lift  coefficient  in  his  numerical  model  to  match  results  of  physical  model 
tests  using  a  small  diameter  stone,  0.43  in.  (10.9  mm),  then  attempted  to  use  the  same 
coefficients  to  predict  the  geomorphic  response  with  a  larger  material,  1  in.  (25.4  mm). 
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The  lifting  process  appeared  to  be  underestimated,  meaning  the  lift  coefficient,  Cl  ,  is 
larger  for  larger  stones.  To  maintain  a  balance  between  the  lift  forces  and  the  internal 
fnction,  Cl  would  have  to  go  up  if  tan  (j)  went  up.  Dai  and  Kamel  (1969)  observed 
similar  discrepancies  with  the  smallest  stones  they  tested  and  concluded  that  the  drag 
coefficient,  Cp  must  be  higher  with  the  smaller  stones.  Again,  to  maintain  a  balance 
between  internal  friction  and  drag,  Cp  would  have  to  increase  as  the  tan  (j)  decreased.  In 
these  examples,  the  resulting  observations  could  as  easily  be  explained  by  inaccurate 
representation  of  the  friction  angle,  that  is  assuming  a  constant  value,  as  errors  in  the  drag 
or  lift  coefficients. 

The  absolute  value  of  the  angle  of  repose  is  questionable  and  determined  by  many 
means,  but  assuming  that  for  small  experimental  stones  a  factor  could  exist  regarding  this 
effect,  an  effort  has  been  made  to  quantify  such.  Clearly  many  scaling  problems  are 
directly  attributed  to  the  Reynolds  number,  while  some  are  a  function  of  the  angle  of 
repose. 
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3  Methodology-Physical  Model  Tests 


Physical  model  approaches  are  often  used  in  navigation  effects  studies  due  to  the 
highly  3-D  effects  near  the  vessel,  the  complexity  of  the  phenomena,  and  the  lack  of 
verified  numerical  solutions.  Physical  modeling  cannot  answer  all  questions  regarding 
the  forces,  nor  can  it  address  the  manner  in  which  nature  couples  these  forces,  but  it  is  a 
viable  means  of  systematically  evaluating  them.  Simple  analytical  approaches  based 
largely  on  empiricism  are  slowly  evolving  into  more  complex  and  intricate  solutions  as 
modeling  capabilities  regarding  instrumentation  are  improved  and  as  the  need  increases 
to  define  more  precisely  the  interaction  of  these  forces  with  the  surrounding  environment. 

Similitude 

Ideally  Reynolds  values  would  dictate  model  size,  but  often  available  space  and  cost 
govern  the  selection.  Testing  conducted  for  this  study  varied  from  a  1 :37.5  to  1 :20  model 
scale  based  on  (a)  the  literature  review,  (b)  results  of  previous  navigation  effects  studies, 
and  (c)  the  size  of  the  available  flume.  This  study  used  previous  riprap  stability  data 
collected  for  the  Tennessee-Tombigbee  study  at  a  1:20  scale.  Additional  stability  data 
were  collected  at  a  1 :25  scale.  Hydrodynamic  data  were  collected  and  analyzed  from 
three  different  model  scales:  1  ;37.5, 1 :25,  and  1 :20. 
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The  equations  of  hydraulic  similitude  for  Froude  number  scaling  are  provided  in 
Table  3. 


Table  3 

Similitude  Relationships 

Characteristic 

Dimension 

Scale  Relations 

Model :  Prototype 

Length 

L=L 

1:20 

1:25 

1:37.5 

Area 

A=L,2 

1:400 

1:625 

1:1,406.3 

Velocity 

1:4.472 

1:5.000 

1:6.124 

Time 

T=4'^2 

1:4.472 

1:5.000 

1:6.124 

Discharge 

1:1,788.9 

1:3,125.0 

1:8,611.5 

Weight 

W=L,3 

1:8,000 

1:15,625 

1:52,734 

Roughness 

Coefficient 

N=Lr'^6 

1:1.648 

1:1.710 

1:1.830 

Physical  Model  Tests 

Riprap  stability  tests 

Data  regarding  riprap  stability  came  from  four  physical  model  studies  conducted  at 
WES,  the  Tennessee-Tombigbee  Waterway  (Maynord  and  Oswalt  1986),  the  Gallipolis 
Lock  approach  study  for  the  U.S.  Army  Engineer  District,  Himtington,  and  two  separate 
data  sets  from  research  sponsored  by  the  Headquarters  U.S.  Army  Corps  of  Engineers, 
Washington,  DC,  in  1990  and  1993.  The  primary  objective  of  these  tests  was  to  identify 
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the  most  severe  operational  conditions  in  which  a  specific  gradation  of  riprap  remained 
stable.  Secondarily,  a  relationship  was  sought  to  link  these  operations  (sailing  speed  and 
blockage  ratio)  to  their  effects  (the  return  current  and/or  waves). 

The  testing  program  covered  a  range  of  conditions.  Underway  test  conditions 
represented  normal  navigation  operations  and  assessed  the  stability  of  stone  slope 
protection  in  a  confined  channel  for  tows  traveling  at  various  speeds,  drafts,  and  distances 
from  the  bank.  All  tests  were  conducted  under  quiescent  pool  conditions;  therefore  all 
flow  disturbances  were  those  created  by  the  movement  of  the  tow.  A  summary  of  details 
regarding  the  Teimessee-Tombigbee,  Gallipolis,  and  the  1990  data  set  are  found  in  Martin 
(1992). 

In  general,  various  combinations  of  the  following  prototype  conditions  were  tested 
as  follows: 

a.  Depth  of  pool:  14  to  22  ft  (4.3  to  6.7  m). 

b.  Sailing  line  (distance  from  tow  C/L  to  toe  of  slope):  52.5  to  300  ft  (16.0  to 
91.4  m). 

c.  Speed  of  tow:  5  to  1 5  J^s  (1 .5  to  4.6  mps). 

d.  Configuration  of  barges:  2-  and  3-barge  widths,  1-  to  5-barge  lengths. 

e.  Lead  barges:  raked  and  square. 

/  Bank  slopes:  2H:1V,  3H:1V. 

g.  Equivalent  spherical  diameter  of  stone  for  which  50  percent  is  lighter  by  weight, 
D50 :  3.25  to  13  in.  (1.3  to  5.1  cm). 
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h.  Weight  of  stone  for  which  50  percent  is  lighter  by  weight,  W50:  1 .7  to  1 1 1  lb 
(0.8  to  50.3  kg).. 

The  cross  sections  in  each  test  series  and  other  pertinent  geometric  data  in  prototype 
dimensions  are  shown  in  Figures  1 1-14  for  the  Tennessee-Tombigbee,  the  Gallipolis,  and 
two  sets  of  navigation  research  studies,  respectively.  The  data  sets  are  identified  as: 

NR93  for  navigation  research  conducted  in  fiscal  year  1993;  NR90  for  navigation 


Figure  11.  Pertinent  cross-section  and  test  data,  Tennessee-Tombigbee  study 
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Figure  12.  Pertinent  cross-section  and  test  data,  Gallipolis  study 


Figure  13.  Pertinent  cross-section  and  test  data,  navigation  research,  FY90 
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Figure  14.  Pertinent  cross-section  and  test  data,  navigation  research,  FY93 


research  conducted  in  fiscal  year  1990;  GAL  for  Gallipolis  lock  approach  studies;  and  TT 
for  Tennessee-Tombigbee  model  studies.  Eight  different  gradations  of  riprap  v^ere  tested. 
Tables  4-7  contain  the  prototype  gradations  tested  in  all  the  studies.  Table  8  provides  a 
summary  of  the  model  dimensions  for  the  eight  gradations. 

A  nonporous  slope  condition  was  modeled  so  no  seepage  could  occur.  The  model 
riprap  was  placed  on  filter  fabric  at  a  thickness  of  approximately  Djoo  •  Model  stone  sizes 
were  selected  based  on  gradations  recommended  in  the  Engineering  Technical  Letter 
(ETL)  120-2-120,  (HQUSACE  1971).  Actual  model  gradations  in  relationship  to  the 
ETL  guidance  are  shown  in  Figures  15-18. 
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Table  4 

Gradations  for  the  Tennessee-Tombigbee  Study 


Thickness,  in. 


D50,  in. 

Yr=  166  pcf 


Percent 
Finer  by 
Weight 


Type  II 
Gradation 


TypeX 

Gradation 


Type  II 


Weight,  Size* 

lb  in. 


Type  X 


Weight  Size 

lb  in. 


0 

0 

360.0 

19.3 

170.0 

80 

170.0 

15.0 

75.0 

60 

92.0 

12.2 

37.0 

50 

68.0 

11.1 

27.0 

40 

50.0 

10.0 

21.0 

30 

33.0 

8.7 

13.5 

20 

21.0 

n 

8.8 

10 

9.0 

5.6 

5.6 

0 

3.6 

4.2 

3.6 

Type  Y 
Gradation 


Type  Y 


Weight  Size 

lb  in. 


’  Equivalent  spherical  stone  diameter  using  Equation  52. 


Table  5 

Gradations  for  the  Gallipolis  Study 

Original 

Proposed 

Thickness,  in. 

24.0 

12.0 

D50,  in. 

Yr=  167  pcf 

13.0 

6.0 

Percent 

Finer  by 

Weight 

Original 

Proposed 

Weight,  lb 

Size,  in. 

Weight,  lb 

Size,  in. 

100 

347.0 

19.0 

41.7 

80 

217.0 

16.3 

25.0 

60 

139.0 

14.0 

14.0 

6.5 

50 

111.0 

13.0 

10.9 

6.0 

40 

87.0 

12.0 

8.3 

5.5 

30 

64.0 

10.8 

6.2 

5.0 

20 

40.0 

9.2 

4.8 

4.6 

10 

25.0 

7.9 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

0 

10.9 

6.0 

0.6 

2.3 
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Table  6 

Gradations  for  the  Navigation  Research,  FY90 


Type  1 


Type  2 


Thickness,  in. 

9.0 

9.0 

D50,  in. 

Yr=  167  pcf 

3.3 

5.0 

Percent 
Finer  by 
Weight 


Type  1 


Type  2 


Type  3 


12.0 


5.8 


T5Tpe  3 


Weight  Size 

Ib  in. 


3.1 

3.9 

1.7 

3.3 
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Table  7 

Gradations  for  the  Navigation  Research,  FY  93 
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Table  8 

Characteristics  of  Model  Riprap 

Model  Study 

Type 

Name 

Yr 

Ib/tf 

D50 

ft 

Wso 

lb 

Thickness 

ft 

NR90 

1 

167 

0.011 

1.2x10-^ 

0.03 

NR90 

2 

167 

0.017 

4.3x10-^ 

0.03 

NR90 

3 

167 

0.019 

5.99x10-^ 

0.04 

GAL 

orig 

167 

0.043 

0.08 

GAL 

prop 

167 

0.020 

6.99x10-^ 

0.04 

XT 

II 

166 

0.047 

9x10-^ 

0.075 

XT 

X 

166 

0.034 

3.4x10-^ 

0.063 

TX 

Y 

166 

0.027 

1.7x10-^ 

0.042 

NR93 

GI 

167 

0.011 

1.2x10-^ 

0.03 

NR93 

G2 

167 

0.019 

5.99x10"* 

0.04 

NR93 

G3 

167 

0.027 

1.7x10-^ 

0.05 

NR93 

G4 

167 

0.037 

4.4x10-^ 

0.06 
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Figure  15.  Model  stone  gradations  used  in  Tennessee-Tombigbee  study 


PERCENT  PASSING  |  d  |  PERCENT  PASSING 


re  17.  Model  stone  gradations  used  in  navigation  research  testing,  FY90 


Figure  18.  Model  stone  gradations  used  in  navigation  research,  FY93 
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The  stability  tests  were  conducted  by  setting  up  the  stone  slope  protection  in  a 
specific  channel  cross  section  and  water  depth.  The  tow  was  then  configured  and 
"driven"  or  towed  past  the  slope  protection  at  a  specific  sailing  line  and  speed.  Stability 
of  the  stone  was  observed  after  numerous  passages  of  the  tow  at  these  conditions.  In  fact, 
each  test  series  was  repeated  generally  100  times  unless  failure  occurred  prior  to 
completion  of  the  repetitions.  Conducting  this  number  of  tests  was  recommended  by 
Fuehrer,  Romisch,  and  Engelke  (1981)  in  their  slope  stability  studies  to  fully  expose  the 
riprap  to  a  sufficient  number  of  waves.  This  finding  corresponded  with  coastal  studies  of 
breakwaters  where  a  minimum  of  1,000  waves  to  a  maximum  of  7,000  are  used  during 
testing.  The  wave  train  accompanying  and  following  the  boat  had  multiple  waves 
depending  on  vessel  length  and  energy,  making  the  number  of  waves  impacting  the  slope 
after  100  vessel  passages  the  same  order  of  magnitude  as  in  coastal  studies.  This  implied 
that  a  minimum  of  1 00  tests  were  required  to  establish  one  data  point  in  the  analysis  of 
stable  riprap. 

Test  conditions  for  the  riprap  stability  studies  are  found  in  Appendix  B.  Each  line 
represents  a  set  of  conditions  including  barge  characteristics,  channel  characteristics, 
bank  slope,  stone  size,  and  the  vessel  speed  at  which  the  model  riprap  embankment  was 
repeatedly  passed  by  the  vessel.  Stability  is  determined  through  observation  of  the  stone 
on  the  embankment  during  testing.  Failure  (F),  also  called  incipient  failure,  is  defined,  m 
this  study  and  others,  as  the  condition  in  which  the  filter  cloth  is  exposed  (Maynord,  Ruff, 
and  Abt  1989).  Stable  (S)  indicates  that  while  some  individual  stones  may  have 
overturned  or  moved,  the  general  thickness  of  the  riprap  was  not  affected.  A  category. 
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marginal  (M),  was  added  for  the  research  work  unit.  During  these  tests,  relatively  small 
and  previously  untested  stone  sizes  were  used.  Although  the  filter  fabric  was  not 
exposed,  the  riprap  thickness  showed  signs  of  thinning  in  the  zone  of  wave  action.  A 
safety  factor  will  be  added  back  into  the  stone  size  for  stable  design  conditions. 

Other  tests 

In  addition  to  the  tests  conducted  for  stability,  large  amounts  of  hydrodynamic  data 
were  collected  to  determine  suitable  predictive  equations  for  the  navigation  effects.  In 
addition  to  data  collected  for  the  Tennessee-Tombigbee  and  the  Gallipolis  study,  wave 
and  velocity  data  collected  included  the  following: 

a.  Tests  were  conducted  for  a  prototype  rectangular  channel,  400  ft  (122  m)  wide 
and  20  ft  (6. 1  m)  deep,  at  three  different  model  scales  (1 :37.5, 1 :25, 1 :20)  to 
determine  associated  problems  with  smaller  scales.  Figure  19  is  a  schematic  of 
the  prototype  dimensions  tested  in  the  model.  Wave  and  velocity  data  were 
collected  for  variable  boat  speeds  at  approximately  three  lateral  locations. 
Velocity  data  were  also  collected  in  three  vertical  positions.  Hydrodynamic 
data  included  boat  speed,  lateral  location  of  the  measurement  from  the  bank, 
vertical  location  of  the  measurement  from  the  bottom,  the  absolute  value  of  the 
maximum  return  current  (x-direction),  the  absolute  values  of  the  lateral  peak 
velocities  at  bow  and  stem,  and  the  maximum  drawdown.  This  test  series  is 
referred  to  as  the  scale  tests  (SCT). 
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b.  In  the  fiscal  year  1993  navigation  research  program,  wave  and  velocity  data 
were  collected  for  most  conditions  in  which  riprap  stability  were  evaluated,  and 
for  some  conditions  (different  boat  draft,  depth,  etc.)  not  relevant  to  actual 
stability  tests.  These  data  were  used  to  develop  regression  equations  for 
obtaining  drawdown  and  return  current  values  for  each  test  scenario  in  the 
stability  analysis.  Further  details  will  follow  in  the  next  chapter. 

Model  Description  and  Appurtenances 

The  model  fixime  dimensions  for  all  tests  except  the  Tennessee-Tombigbee  tests  are 
shown  in  Figure  20.  The  maximum  water  depth  that  could  be  tested  in  the  flume  was 
1  ft.  Data  were  collected  in  an  approximately  100-ft-  (30.5-m-)  long  area  where  the  floor 
was  level.  The  flume  was  not  equipped  with  a  recirculating  system;  therefore,  all  tests 
were  conducted  in  quiescent  conditions.  The  flume  was  equipped  with  a  winch  and 
pulley  system  for  the  fiscal  year  93  navigation  tests.  The  sailing  line  was  adjusted  by 
moving  the  winch  and  cable  assembly  across  the  width  of  the  usable  flume  area.  In  some 
cases,  the  location  of  columns  in  the  facility  limited  the  number  of  sailing  lines  that  could 
be  tested. 
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Figure  20.  Dimensions  of  model  flume  area  for  navigation  effects  testing,  plan  view 


Riprap  stability 

The  riprap  gradations  used  in  the  1993  tests,  were  based  on  previous  tests  and 
observed  failure  characteristics.  Four  gradations  were  selected  ranging  in  model  size 
from  a  D50  of  0.0108  to  0.0370  ft  (3.3  to  1 1.3  mm)  with  corresponding  prototype 
gradations  as  specified  in  Table  7.  A  layer  of  sand  was  placed  on  the  slope  and  covered 
with  filter  fabric.  The  scaled  riprap  was  placed  at  the  specified  thickness  from  the  toe  of 
the  slope  to  well  above  the  waterline.  Observations  of  movement  were  recorded  for  areas 
unaffected  by  boundary  edges.  Figure  21  shows  the  layout  of  the  slopes  and  gradations 
observed  for  this  test  series  in  model  dimensions.  The  riprap  was  lightly  painted  after 
placement  to  help  identify  the  direction  of  motion  and  the  origin  of  displaced  stones. 
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Figure  21 .  Model  dimensions  and  layout  of  stability  tests 


Hydrodynamic  data 


Data  collected  for  this  study  included  the  tow's  speed,  the  water-surface  movement, 
and  the  water  velocity  at  various  locations  adjacent  to  the  tow.  All  analog  data  were 
collected,  converted,  and  saved  digitally  on  DOS-based  personal  computers.  The  tow 
speed  was  monitored  through  the  test  section  by  interrupting  a  retro-reflective  light  beam 
with  1-ft  panels  attached  to  the  top  of  the  moving  barges.  All  wave  data  were  collected 
using  capacitance  wave  rods.  No  more  than  two  rods  were  used  at  once.  Sampling 
frequency  was  generally  25  hz.  The  location  of  the  instruments  varied  with  testing. 
Velocity  data  over  the  course  of  these  studies  have  been  collected  with  a  variety  of 
instruments.  An  electromagnetic  meter  measured  return  currents  in  the  Tennessee- 
Tombigbee  study.  An  impeller  type  probe  was  used  in  the  Gallipolis  study  to  pick  up  a 
few  velocity  measurements.  Recent  testing  for  the  Navigation  Hydraulics  program 
included  using  a  two-dimensional  (2-D)  laser  doppler  velocimeter  (LDV)  and,  near  the 
end  of  testing,  newly  developed  3-D  acoustic  doppler  velocity  (ADV)  probes.  The 
location  and  depth  of  velocity  measurements  were  varied  during  testing. 

The  vessel 

The  tests  conducted  for  the  Tennessee-Tombigbee,  the  Gallipolis,  and  the  first  set  of 
navigation  research  tests  used  a  free  moving  model  tow.  The  tow  in  the  Tennessee- 
Tombigbee  model  was  a  1 : 20-scale  self-powered  vessel  operated  by  a  technician  on 
board;  in  the  Gallipolis  and  NR90  studies,  it  was  a  1 :25-scale  radio-operated  tow.  For  the 
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NR93  tests,  a  towing  cable  and  winch  pulled  the  tow  model  through  the  test  section.  This 
included  pulling  the  barges  (three  wide  by  two  long)  and  the  model  towboat  while 
operating  the  propellers  during  the  stability  tests. 

The  operation  of  the  propellers  was  based  on  tow  thrust  calculated  as  a  percentage 
of  a  boat's  bollard  push.  Toutant  (1982)  defines  bollard  push  as  the  dock  or  zero-speed 
force  of  the  vessel  and  presented  formulas  for  calculating  this  force  as  a  function  of 
horsepower.  The  equations  presented  by  Toutant  were  used  to  determine  an  appropriate 
bollard  push  for  the  prototype  towboat  being  modeled.  To  obtain  a  maximum  bollard 
push,  the  prototype  was  assumed  to  have  5,600  hp;  therefore,  the  maximum  bollard  push 
was  105,500  lb  (4.7  x  10^  N)  or  6.75  lb  (30  N)  in  the  model.  The  model  towboat  was 
built  as  a  replica  of  the  U.S.  Army  Engineer  District,  Vicksburg,  working  boat,  the 
Benyaurd.  According  to  conversations  with  personnel  at  the  U.S.  Army  Engineer 
District,  Huntington,  and  the  pilot  of  the  Benyaurd,  maximum  rpm's  for  this  size  towboat 
is  200  to  220  which  equals  approximately  1,000  to  1,100  rpm's  in  the  model.  Based  on 
the  same  conversations  with  the  Huntington  District,  operators  in  open  river  use 
approximately  75  to  80  percent  throttle.  Then  at  maximum  operation,  the  bollard  push 
should  be  approximately  0.75  x  6.75  =  5.1  lb  (23  N).  From  graphs  obtained  from  testing 
the  force  produced  by  the  tow  at  various  rpm's,  the  corresponding  operation  for  the  model 
boat  should  be  approximately  860  rpm  (172  rpm  prototype).  Likewise  at  50  percent 
throttle,  the  bollard  push  would  be  3.38  lb  (15  N),  and  the  corresponding  rpm  is  720  (144 
prototype).  Operation  of  the  model  tow  using  the  prescribed  method  for  determining  the 
appropriate  rpm,  resulted  in  tow  speeds  that  were  too  slow  due  to  viscous  effects  in  the 
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model.  Operating  the  free  moving  towboat  and  barges  using  only  the  boat’s  self 
propulsion  resulted  in  the  prototype  values  of  boat  speed  shown  in  Table  9. 


Table  9 

Vessel  Thrust  and  Force 

Thrust 

% 

Force 

lb 

Tachometer 

rpm 

Boat  Speed 
fps 

10 

10,625 

74 

2.50 

25 

26,250 

112 

4.25 

50 

52,813 

152 

6.00 

75 

79,688 

172 

7.00 

To  achieve  the  testing  program  objectives,  vessel  speeds  needed  to  approach  their 
limiting  speed.  In  the  fiscal  year  1993  test  channel  at  a  20-ft  (6.1-m)  prototype  water 
depth,  the  prototype  limit  velocity  is  approximately  14.9  J^s  (4.5  mps).  All  propeller 
speeds  were  set  proportional  to  the  vessel  speed  used  in  testing.  To  overcome  the  viscous 
effects  and  get  the  vessel  to  operate  closer  to  the  higher  speeds,  the  tovsdng  cable 
assembly  was  used  for  the  folloiving  two  basic  reasons: 

a.  It  improved  the  repeatability  of  the  testing  particularly  regarding  the  sailing  line. 
Also,  it  helped  maintain  a  constant  tow  speed  that  was  not  influenced  by 
depleting  battery  supplies,  interference  in  radio  transmitted  tachometer  signals, 
and  the  like. 
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b.  It  allowed  a  wider  range  of  tow  speed  operations  without  unduly  overpowering 
the  propulsion  from  the  propellers. 

The  propellers  were  operated  to  insure  that  any  influence  of  propeller  turbulence 
was  not  disregarded  as  pertains  to  slope  stability.  As  it  turns  out,  the  primary  area  where 
the  propellers  could  influence  the  bank  stability  is  when  the  tow's  edge  operates  over  the 
toe  of  the  slope.  In  this  area,  the  turbulence  from  the  propellers  has  a  localized  influence 
on  stone  movement  at  the  slope  bottom  due  to  the  proximity  of  the  jet  to  the  channel 
bottom  and  the  potential  for  steering  maneuvers  to  keep  away  from  the  bank.  However, 
under  most  sailing  line  conditions,  the  hull  displacement  generates  the  critical  forces. 

Barges  for  all  studies  were  made  of  sheet  metal,  ballast  with  weights,  and  connected 
with  C-clamps.  Each  barge  used  in  the  l:20-scale  Termessee-Tombigbee  study  and  for 
the  scale  tests  was  9.75  ft  (3.0  m)  long  by  1.75  ft  (0.8  m)  wide  and  had  a  curved  rake  at 
the  bow.  Two  different  configurations  were  used  in  the  Tennessee-Tombigbee  study;  a 
two-wide-by-one  long  unloaded  and  a  three-wide-by-two-long  loaded  configuration.  The 
configuration  for  all  the  scale  tests  was  three  wide  by  two  long.  The  rakes  on  these 
barges  varied  according  to  the  scale  barges  used.  Barges  used  at  the  1 :37.5  scale  had  a 
configuration  of  three  wide  by  two  long  with  a  rake  that  made  an  angle  of  26.5  deg  with 
the  horizontal.  Barges  in  the  Gallipolis  study  had  a  square  bow  and  stem,  and  the 
configurations  varied  from  a  three  wide  by  two  long  to  a  three  wide  by  five  long.  Barges 
in  the  1 :25-scale  riprap  tests  were  7.8  ft  (2.4  m)  long  by  1.4  ft  (0.4  m)  wide.  The 
configuration  used  for  these  tests  was  three  wide  by  two  long.  The  lead  barges  had  a  rake 
that  made  a  45  deg  angle  with  the  horizontal.  The  rear  barges  had  square  ends. 
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4  Analysis  and  Results 


Data  Collection  and  Uncertainty 


Problems  inherent  to  physical  model  testing  include  scale-effects  phenomena, 
difficulties  with  the  flume  itself,  accuracy  of  the  instrumentation,  location  of  the 
instrumentation,  and  whether  the  tow  is  free-moving  or  towed,  etc.  The  variabilities  that 
arise  from  these  problems  can  cause  uncertainties  in  the  data  collection.  Every  effort  was 
made  to  correct  for  these  problems  and  minimize  their  effects. 

Flume  limitations 

The  length  of  the  test  facility  (Figure  20)  used  in  the  NR93  testing  limited  the  tow 
size  that  could  be  modeled,  the  length  of  time  available  to  capture  the  event,  and  the 
ability  to  get  the  vessel  operating  at  full  speed.  The  maximum  length  of  tow  modeled 
was  two  barge  lengths.  As  the  tow  starts,  particularly  when  towed,  a  celerity  wave  is 
generated  that  reflects  off  the  ends  of  the  flume.  The  travel  time  of  this  wave  can 
interfere  with  the  time  history  of  the  data  collection  since  it  is  a  function  of  the  vessel 
speed  and  water  depth.  The  startup  time  and  distance  to  testing  apparatus  was  evaluated 
for  each  test  condition  so  that  the  longest  time  history  response  could  be  obtained  prior  to 
the  bow  passing  the  instrument  and  following  the  stem.  Interference  due  to  reflective 
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waves  was  unavoidable  in  some  conditions.  Data  were  collected  near  the  center  of  the 


flume  where  the  vessel  was  operating  at  a  relatively  uniform  speed  to  minimize  potential 
accelerations  and  decelerations  of  the  vessel. 

Bow  wave  phenomena 

The  amplitude  of  the  bow  wave  could  be  effected  by  both  the  timing  of  the  celerity 
wave  and  the  influence  of  the  towing  system.  In  the  examination  of  wave  data,  the 
reflective  wave  and  the  bow  wave  were  often  indistinguishable.  The  tow  travel  time  and 
starting  position  were  adjusted  when  possible  so  that  the  celerity  wave  could  be  easily 
distinguished  from  the  bow  wave.  In  Tenaud's  testing  (1977)  conducted  at  a  1:20  scale, 
the  tow  was  attached  to  a  "tractor,"  or  mechanical  wheel  on  the  floor  that  guided  the  boat 
along  the  sailing  line,  and  was  pulled  by  a  cable.  He  observed  a  substantial  wave  at  the 
bow,  also  observed  in  subsequent  studies  conducted  at  WES  using  both  a  towing  cable 
and  a  towing  carriage.  This  is  likely  an  artifact  of  the  rigidity  of  the  towing  system  and 
its  acceleration.  It  is  also  likely  in  Fuehrer,  Romisch,  and  Engelke's  (1981)  testing  that 
the  towing  apparatus  in  their  1 :60  scale  model  caused  an  exaggeration  of  the  bow  wave. 
In  fact,  their  stability  equation  was  based  on  this  phenomena.  Conversely,  bow  waves  in 
tests  of  free-running  model  tows,  such  as  that  used  on  the  Teimessee-Tombigbee  study 
(Maynord  and  Oswalt  1986),  and  prototype  data  collected  on  the  Illinois  River  and  Upper 
Mississippi  River  indicate  that  a  rise  preceding  drawdown  at  the  bow  is  imperceptible. 
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Drawdown  discrepancies 


When  Tenaud  (1977)  compared  his  measurements  of  drawdown  to  Schijf  theory 
(Jansen  and  Schijf  1953),  he  fovmd  them  close  for  boat  speeds  less  than  90  percent  of  the 
limiting  speed.  For  higher  speeds  and  drafts,  the  water  level  drops  that  were  measured 
were  lower  than  those  calculated.  A  similar  observation  was  made  during  scale  testing  at 
WES.  At  the  1  ;25  scale  and  1 :37.5  scale,  model  data  were  higher  than  Schijf  values  and 
at  the  1 :20  scale,  values  were  lower  at  higher  speeds.  Without  any  other  reasonable 
explanation,  the  towing  cable  and  assembly  may  be  responsible  for  this  discrepancy. 
Tenaud  also  observed  that  measured  values  are  greater  than  the  estimated  values  for  all 
tests  using  a  lighter  draft. 

Variability  of  data  scatter 

Variability  in  the  data  can  occur  for  multiple  reasons;  therefore,  quality  control  was 
exercised  to  minimize  those  variabilities.  One  quality  control  method  is  to  replicate 
tests.  To  improve  consistency  of  results  and  evaluate  variability,  tests  were  replicated 
during  the  collection  of  wave  and  velocity  data  at  each  measuring  point.  Figures  22 
and  23  show  typical  scatters  of  return  currents  and  drawdowns  for  multiple  boat  speeds 
and  repetitions  collected  during  the  NR93  testing.  The  prototype  test  conditions  for  these 
measurements  were  tow  sailing  on  the  centerline  of  the  channel,  depth  of  water,  20  ft  (6.1 
m),  and  a  2H:  1 V  bank  slope.  Another  method  to  insure  repeatability  of  tow  passage  for 
both  the  vessel  speed  and  sailing  line  is  to  use  the  towing  cable.  In  spite  of  efforts  to 
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Figure  22.  Variability  in  measuring  return  currents  in  model 


Figure  23.  Variability  in  measuring  drawdown  in  model 
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minimize  irregularities  between  vessel  passage,  some  problems  still  occurred  with  tension 
on  the  cable,  motor  speed  controllability,  operator  variations,  etc.  Since  velocity 
measurements  were  taken  with  different  instruments  over  the  course  of  the  testing 
program,  the  data  results  were  compared  and  found  favorable.  Figure  24  shows  the 
comparison  of  values  taken  for  the  same  test  conditions  as  Figure  22  using  both  the  LDV 
and  ADV  to  measure  velocity.  Although  propeller  influence  was  determined  to  be 
minimal  unless  the  sailing  line  is  very  close  to  the  toe  of  the  slope,  the  propellers  were 
operated,  whenever  possible,  to  more  closely  represent  the  actual  vessel  passage. 
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Figure  24.  Comparison  of  measurements  taken  with  LDV  and  ADV 
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Lateral  variability  of  drawdown  and  return  current 

Retim  current  and  drawdown  data  were  collected  halfway  between  the  edge  of  the 
tow  and  the  toe  of  the  slope.  It  was  assumed  that  these  variables  were  laterally  uniform 
from  the  vessel  to  the  shore.  Recall,  this  is  true  for  uniform  channels  of  small  blockage 
ratios  (Jansen  and  Schijf  1953).  Since  riprap  is  generally  placed  in  uniform  channels 
such  as  narrow  canals  or  lock  approaches,  this  assumption  is  valid.  Tenaud  (1977)  also 
confirmed  this  assumption  stating  that  measured  wave  heights  at  several  locations  gave 
"practically  identical  results."  Results  of  hydrodynamic  data  collected  at  the  three  lateral 
locations  shown  in  Figure  19  during  the  l:20-scale  tests  (SCI)  support  this  assumption. 
The  data  collected  at  1 10.75  ft  (33.76  m),  73.75  ft  (22.48  m),  and  37  ft  (1 1.28  m)  in 
prototype  dimensions  from  the  vertical  wall  of  the  channel  are  referred  to  as  “near  tow,” 
“mid-chaimel,”  and  “near  bank,”  respectively.  In  Figures  25  and  26,  where  boat  speed  is 
plotted  against  return  current  and  drawdown,  respectively,  data  scatter  does  not  appear  to 
follow  any  particular  trend  related  to  the  three  locations  between  the  vessel  and  the  bank. 

Vertical  variability  of  velocity 

An  assumption  typically  made  regarding  return  currents  is  that  the  velocity  profile  is 
uniform  over  the  depth  of  the  water  column.  This  is  more  pronounced  since  the  boundary 
layer  has  not  fully  developed.  Figure  27  shows  velocities  taken  at  three  vertical  locations 
in  a  20-ft  (6.1-m)  prototype  channel  during  the  l:25-scale  testing.  More  data  are 
currently  being  collected  to  verify  this  assumption  particularly  in  nommiform  channels. 
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Figure  25.  Variability  in  the  lateral  distribution  of  return  currents 


Figure  26.  Variability  in  lateral  distribution  of  drawdown 
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Figure  27.  Vertical  variability  in  return  currents,  mid-channel,  depths  are  from  bottom 


Hydrodynamic  Data  Analysis 

Selection  of  peak  values 

Time  history  hydrodynamic  data  of  waves  and  currents  collected  during  navigation 
effects  testing  were  analyzed  and  peak  values  were  extracted.  The  test  conditions,  peak 
velocity,  and  wave  data  are  found  in  Appendixes  C  and  D  for  velocities  and  waves, 
respectively.  They  contain  the  test  series,  model  scale,  channel  dimensions,  sailing  line, 
vessel  width,  draft,  speed,  location  of  measurements,  and  all  peak  data. 
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The  maximum  return  current,  u^^  ,  in  the  x-direction  velocity  was  taken  from  the  data 
and  is  defined  according  to  the  definition  in  the  literature  given  by  PIANC  (1987)  and 
Maynord  (1996).  A  comparison  of  calculated  values  of  return  current  using  the  Schijf 
approach  with  a  coefficient,  a; ,  of  1  to  measured  values  of  u,^  from  the  SCT  and  NR93 
data  sets  is  found  in  Figure  28.  Schijf  values  are  included  in  the  summary  table  in 
Appendix  C.  Comparisons  of  TT  velocity  data  to  theory  are  found  in  Maynord  and 
Oswalt  (1986). 


Figure  28.  Comparison  of  measured  values  of  return  current  to  computed  values  using 
the  Schijf  approach 
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Since  velocities  v^ere  collected  in  two  dimensions,  plan  view,  in  the  NR93  and  SCT 
data,  the  maximum  lateral  velocities  were  also  recorded.  An  absolute  peak  value  away 
from  the  bow,  Uyj, ,  and  one  toward  the  stem,  Uyj ,  are  included  in  Appendix  C  in 
prototype  dimensions.  The  magnitude  of  these  velocities  is  small  relative  to  the  peak 
return  current  as  shovra  in  Figure  29. 


Figure  29.  Comparison  of  longitudinal  currents  to  lateral  currents  produced  by  tow 


Four  different  pieces  of  peak  wave  data  were  extracted  from  the  time  history  plots, 
the  bow  wave,  hf ,  the  maximum  drawdown,  2^^ ,  the  maximum  peak  of  the  secondary 
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waves,  Hj ,  and  a  variation  of  the  transversal  stem  wave,  (defined  as  the  trough  to 

following  peak  of  the  transversal  stem  wave).  differs  slightly  from  in  that  it 
represents  the  actual  height  of  the  wave  at  the  stem  rather  than  the  height  of  the  trough  of 
the  water  level  depression  back  to  the  still  water  level,  as  shown  in  Figure  4.  The  bow 
wave  data  were  not  used  in  any  of  the  analysis  due  to  the  potential  error,  already 
discussed,  imposed  by  the  towing  cable  and  assembly. 

To  separate  long  period  responses  fi-om  the  NR93  time  history  data,  a  moving  average 
was  used  to  filter  out  the  high  frequency  response  of  the  short  period  waves.  The 
maximum  drawdown,  ,  and  maximum  wave  height  at  the  stem  of  the  vessel,  H^ax  > 
were  determined  from  plots  of  the  smoothed  time  histories.  Short  period  waves  were 
obtained  using  a  Fourier  transform  method.  Essentially,  the  data  are  transformed  firom 
the  time  domain  to  the  frequency  domain;  low  frequency  data  are  filtered  out;  and  the 
remaining  high  frequency  data  are  converted  back  to  the  time  domain  (Press  et  al.  1986). 
Hj  is  the  maximum  short  period  response  in  the  high  frequency  wave  spectrum.  Figure  5 
shows  typical  filtered  and  unfiltered  wave  responses.  The  TT  data  set  was  not  filtered, 
and  Hn,ax  actually  represented  either  the  maximum  trough  to  following  peak  of  the 
transversal  stem  wave,  if  present,  or  a  maximum  secondary  wave. 

Figure  30  shows  how  the  secondary  waves  from  the  NR93  data  and  the  maximum 
wave  data  for  all  data  sets  relate  to  the  maximum  measured  drawdown.  As  can  be  seen, 
both  the  secondary  wave  and  the  maximum  wave  are  dispersed  on  either  side  of  the  line 
of  perfect  agreement  with  slightly  more  points  below  the  line.  This  graph  tends  to  show 
that  the  value  of  conservatively  represents  the  peak  wave  in  the  wave  spectrum  for 
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Figure  30.  Relationship  between  secondary  waves,  maximum  wave,  and  maximum 
drawdown 


commercial  tows  in  confined  channels  with  blockage  ratios  less  than  approximately  16. 
For  low  vessel  speeds,  less  than  approximately  8  fys  (2.4  mps),  Hj  was  negligible.  In  an 
unconfined  channel,  for  vessels  traveling  at  high  speeds,  the  secondary  or  hull-formed 
waves  would  dominate  the  wave  spectrum. 

The  measured  maximum  drawdown  for  the  SCT  and  NR93  data  is  compared  to 
values  obtained  by  the  Schijf  method,  again  using  an  aj  of  1,  in  Figure  31.  Comparisons 
of  TT  wave  data  to  theory  are  foimd  in  Maynord  and  Oswalt  (1 986).  The  measured 
secondary  waves  are  compared  to  calculated  values  based  on  Equation  25  in  Figure  32. 


Ill 


3.5 


Figure  3 1 .  Comparison  of  measured  values  of  to  computed  values  of 
drawdown  using  the  Schijf  approach 


Figure  32.  Comparison  of  measured  values  of  Hj  to  computed  values 
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The  coefficients  and  a2,  respectively,  taken  from  Blaauw  et  al.  (1984)  are  0.8  and  2.67, 
and  from  Verhey  and  Bogaerts  (1989)  are  1 .0  and  4.0.  Appendix  D  contains  the  test 
conditions  and  peak  wave  data  from  these  tests  and  Schijf  values  of  drawdown. 

Mapping  hydrodynamic  data  to  stability  tests 

Each  set  of  stability  tests  was  conducted  by  operating  the  vessel  at  a  constant  speed 
for  up  to  100  passes.  Since  wave  heights  and  velocities  were  not  measured  for  each  of 
these  passes,  a  representative  value  corresponding  to  that  test  series  was  needed  to 
analyze  the  relationship  between  navigation  forces  and  bank  stability.  The  hydrodynamic 
data  provided  in  Appendixes  C  and  D  had  to  be  mapped  to  the  test  conditions  presented 
for  the  stability  tests  in  Appendix  B. 

For  most  of  the  data,  regression  analysis  was  conducted  with  the  hydrodynamic  data 
to  relate  model  tow  speed  to  return  current  and  drawdown.  A  return  current  and 
drawdown  were  determined  from  the  regression  corresponding  to  the  channel  conditions, 
sailing  line,  configuration,  and  boat  speed  used  during  each  set  of  stability  tests.  The 
maximum  return  current,  drawdown,  maximum  wave,  and  secondary  wave  for  each  test 
series  were  determined  as  follows: 

a.  If  actual  measurements  were  taken  for  fixed  conditions  of  water  depth,  channel  - 
section,  and  sailing  line,  a  regression  analysis  was  done  of  boat  speed  versus  each 
parameter.  The  best  correlations  were  found  with  second  order  regressions  and 
values  corresponding  to  the  test  were  taken  from  the  regression.  Where  the  curves 
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had  low  regression  constants,  engineering  judgment  was  used  to  select  the  best 
estimate  of  these  values  corresponding  to  the  data. 

b.  For  certain  sailing  lines  and  depths,  data  were  not  collected  on  both  sides  of  the 
tow.  Actual  measured  values  were  weighted  according  to  distributions  from 
Schijf  equations  with  skew  coefficients  (Maynord  1990). 

c.  Not  many  velocities  or  wave  heights  were  collected  for  the  1 990  navigation 
research  data  nor  the  Gallipolis  study.  Values  of  return  current  and  drawdown  for 
these  tests  were  based  on  the  Schijf  equations.  The  values  resulting  from  this 
method  are  low  (as  discussed  in  scale  effects).  Due  to  the  inadequacies  of  these 
measurements,  data  points  were  not  used  in  the  stability  analysis  from  the  NR90 
and  GAL  data  sets. 

d.  In  the  Tennessee-Tombigbee  study,  maximum  wave  height,  j  was  given 
corresponding  to  failure  or  stability.  Boat  speed  data  were  not  given  for  those 
points  but  were  backed  out  from  graphs  in  the  report  relating  boat  speed  to  wave 
height  (Maynord  and  Oswalt  1986).  Drawdown  values  corresponding  to  each 
boat  speed  and  test  condition  were  extracted  from  curves  through  actual 
measurements. 

The  peak  hydrodynamic  data  corresponding  to  each  set  of  stability  tests  are  summarized 
in  Appendix  B.  These  data  are  used  in  the  original  stability  analysis  to  determine 
relationships  between  stone  failure  and  hydrodynamic  characteristics  such  as  waves  and 
currents. 
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Dimensional  Analysis 


Dimensional  analysis  is  often  used  in  experimental  studies  to  develop  predictive 
equations  not  easily  described  by  the  laws  of  nature.  When  coupled  with  experimental 
data,  the  analysis  can  provide  quantitative  results  in  the  form  of  predictive  equations.  The 
complexity  of  the  problem  presented  in  this  study  leads  to  many  dimensionless 
parameters.  Riprap  design  for  tow-induced  forces  is  based  on  essentially  three  sets  of 
variables:  those  related  to  ship  waves,  displacement  motion  of  the  tow,  and  stability  of 
the  stone  on  the  slope  embankments. 

Sorensen  and  Weggle  (1984)  used  dimensional  analysis  to  develop  four 
dimensionless  parameters  related  to  ship  waves.  The  dimensionless  variables  were: 

(1)  The  vessel  Froude  number ,  already  defined  as  Vj  /V (gh) 

(2)  The  dimensionless  wave  height,  H/  V  ,  dependent  variable 

(3)  Dimensionless  distance  from  sailing  line, 

(4)  Dimensionless  depth, 

where  V  is  the  displacement  volume  of  the  ship,  and  S  is  the  distance  from  the  sailing 
line  to  the  point  where  the  wave  height  is  measured. 

Using  dimensional  analysis,  the  author  derived  the  following  terms  for  vessel 
displacement  motion  based  on  the  variables  in  the  Schijf  equations. 

(1)  (Hg)Ws^,  dependent  variable 

(2)  Vj// (gh)  ,  vessel  Froude  number 

(3)  A  j/Ab ,  blockage  ratio 
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Hudson  (1957)  stated  in  his  analysis  of  pervious  breakwaters,  where  k  is  wave 
length, 


Characteristics  of  the  motion  of  water  particles  when  short-period  wind  waves 
encounter  a  rubble-mound  breakwater  are  determined  by  the  wave  steepness 
(H/A),  the  relative  depth  (h/A),  the  relative  height  (H/h),  the  depth  at  the 
breakwater  slope  (h),  the  angle  of  the  beach  slope  seaward  of  the  breakwater, 
angle  of  seaside  slope  of  the  breakwater  with  the  horizontal  (a),  the  angle  of 
obliquity  of  the  attacking  waves  (P),  and  the  shape,  thickness,  and  porosity  of  the 
cover  layer  and  underlayer  materials. 

Hudson  presented  his  dimensional  analysis  of  rubble-mound  structures  in  Hudson 
(1979).  As  a  result  of  both  the  analysis  of  the  problem  dynamics  and  the  dimensional 
analysis,  and  based  on  several  assumptions,  he  reduced  the  basic  equation  to  the 
foUoAving: 
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a,  shape,  placement,  %  damage 
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(99) 


The  left  side  of  the  equation  is  the  already  defined  dimensionless  parameter,  the  stability 
number,  Hudson  assumed  variations  of  H/h  and  the  Reynolds  number  to  be  minimal 
during  testing  and  did  not  include  them  in  his  analysis.  He  dismissed  the  use  of  repose 
angle  altogether  for  reasons  already  discussed. 

Other  dimensionless  variables  included  in  studies  of  stone  stability  and  vessel  effects 
were  the  repose  angle,  surf  similarity  parameter,  and  draft  vessel  Froude  number.  In  this 
study  the  stability  number  was  plotted  against  many  dimensionless  variables.  The 
Reynolds  number,  repose  angle,  and  bank  slope  were  determined  to  influence  the 
stability  number. 
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Scale  Effects  Analysis 


Scale  effects  can  exist  both  in  the  ability  to  reproduce  the  tow-induced  forces  and  in 
modeling  the  graded  stone.  Boundary  layer  problems,  low  Reynolds  values,  and  the 
angle  of  repose  are  addressed  in  an  attempt  to  minimize  scale  effects  and  make  model 
adjustments.  Corrections  to  Reynolds  or  angle  of  repose  would  leave  fewer 
dimensionless  variables  to  analyze  for  stability. 

Boundary  layer  effects 

It  has  been  discovered  through  recent  navigation  effects  modeling  at  WES,  that  as  a 
result  of  the  misrepresentation  of  viscous  forces  at  the  model  scale,  the  moving  tow 
actually  behaves  as  a  larger  tow.  Due  to  the  boundary  layer  growth,  model  tows  have 
apparent  drafts  much  larger  than  their  equivalent  prototype  tows.  Therefore,  prototype 
vessels  going  the  same  respective  speed  as  the  Froudian-scaled  model  boat  vdll  produce 
lesser  effects;  i.e.,  the  prototype  return  current  and  drawdown  will  be  smaller  for  the  same 
given  tow  configuration  and  speed.  Results  of  scale  effects  testing  revealed  differences 
in  model  data,  best  shown  graphically  in  Figures  33  and  34,  by  comparing  peak  return 
current  and  drawdown  data  between  three  model  scales.  For  equivalent  prototype  vessel 
speeds,  the  maximum  variation  in  return  current  between  a  l:37.5-scale  model  data  and 
1 :20-scale  model  data  was  as  much  as  1.5  fys  (0.46  mps).  For  drawdown,  the  variation 
was  as  high  as  0.57  ft  (0.17  m).  One  would  expect  these  scale  effects  to  be  even  greater 
when  extrapolating  from  say  a  1:20  scale  to  a  1:1  scale  (prototype). 
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Figure  33.  Results  of  scale  effects  testing,  return  current 
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Figure  34.  Results  of  scale  effects  testing,  drawdown 
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For  the  navigation  research,  vessel  speeds  in  the  model  varied  approximately  from 
0.6  ^s  (0.18  mps)  to  2.8  ^s  (0.85  mps),  and  for  the  Tennessee-Tombigbee  study, 

1.3  fps  (0.40  mps)  to  3.6  fps  (1.10  mps).  According  to  Sorensen  (1966),  only  boat  speeds 
greater  than  approximately  1.8  fps  (0.5  mps)  should  be  tested  to  minimize  scale  effects. 

As  a  result  of  these  discrepancies,  design  guidance  regarding  riprap  stability  as  a 
direct  function  of  boat  speed  has  been  avoided  in  this  analysis.  Relationships  in  the 
analysis  were  primarily  sought  relating  the  actual  model  values  of  tow-induced  currents 
and  waves  to  stability.  Boundary  layer  scale  effects  also  prohibited  attaining  limiting  or 
critical  prototype  vessel  speeds  in  the  model;  thus  making  it  essentially  impossible  to 
create  conditions  that  would  fail  the  larger  stones. 

Reynolds  criteria 

Assuming  that  stone  size  is  a  function  of  wave  height  cubed,  as  in  the  literature 
formulations,  implies  that  the  stability  number,  Nj,  is  essentially  constant  for  a  given 
slope.  Some  authors  included  the  angle  of  repose,  surf  similarity  parameter,  or  the 
number  of  waves  but  most  suggest  a  constant  value  of  .  Observations  during  the  NR93 
testing  suggested  smaller  stone  sizes  exhibited  greater  instabilities,  consequently  smaller 
stability  numbers.  These  observations  may  be  attributed  in  part  to  low  model  Reynolds 
values.  Fuehrer,  Romisch,  and  Engelke  (1981)  also  found  Nj  to  vary  with  Reynolds, 
while  conducting  tests  with  very  small  diameter  stones  for  vessel-induced  wave  heights. 
Oumeraci  (1984)  determined  to  be  smaller  at  low  Reynolds  and  recommended 
adjusting  values  of  Nj  to  compensate. 
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One  way  of  evaluating  low  Reynolds  values  is  to  compare  computed  numbers  in  the 
model  to  recommended  values  from  the  literature.  In  summary,  recommended  critical 
values  of  particle  Reynolds  numbers  in  a  wave  environment  varied  according  to  the 
source  as  follows: 

a.  Fuehrer,  Romisch,Engelke  (1981):  1000 

b.  Jensen  (1989):  6000  -  7000 

c.  Dai  and  Kamel  (1969):  3  x  lO'^ 

d.  Hudson  (1957):  3  x  10"* 

e.  Oumeraci  (1984):  3  x  10"^ 

/  Broderick  and  Ahrens  (1982):  10^ 

Particle  Reynolds  values  for  maximum  and  minimum  stone  sizes  and  wave  heights 
were  computed  based  on  Equation  89  for  a  kinematic  viscosity  of  water  at  60“  F  (16“  C). 
The  equation  is  presented  in  several  forms  in  Chapter  2,  and  Table  10  presents  these 
computation  results.  The  largest  computed  Reynolds  value  for  model  stone  and  wave 
height  was  used  in  the  Tennessee-Tombigbee  study.  With  a  model  diameter  of 
0.046  ft  (0.014  m)  and  a  maximum  model  wave  height  of  0.175  ft  (0.046  m),  the 
computed  Reynolds  value  is  9  x  10^.  The  smallest  Reynolds  value  was  less  than  400. 
According  to  these  calculations,  none  of  the  data  would  meet  Dai  and  Kamel's  (1969), 
Hudson's  (1957)  or  Broderick  and  Ahren's  (1982)  criteria.  To  meet  Jensen's  criteria,  the 
minimum  model  wave  height  for  the  smallest  stone  would  have  to  be  1 .42  ft  (0.43  m)! 
All  data,  except  wave  heights  less  than  0.04  ft  (0.01  m)  imposed  on  the  smallest  model 
stone,  have  Reynolds  values  exceeding  1,000. 
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Table  10 

Computed  Values  of  Model  Particle  Reynolds 


IL 


^50 

Prototype 

in. 

■ 

U 

^max 

Model 

ft 

R 

Model 

Navigation  Research 

3.25 

0.011 

0.005 

3.6x102 

3.25 

0.011 

0.160 

1.9x102 

5.8 

0.019 

0.005 

6.3  X  102 

5.8 

0.019 

0.160 

3.5x102 

8.1 

0.027 

0.005 

1.5x102 

8.1 

0.027 

0.160 

5.0x102 

11.1 

0.037 

0.005 

1.2x102 

11.1 

0.037 

0.160 

6.9x102 

Tennessee-Tombigbee  Study 

6.3 

0.026 

0.005 

8.6  X  102 

6.3 

0.026 

0.175 

5.1  X  102 

8.1 

0.034 

0.005 

4.0x102 

8.1 

0.034 

0.175 

6.6x102 

11.1 

0.046 

0.005 

1.5x102 

11.1 

0.046 

0.175 

8.9  X  102 

Reynolds  values  can  be  computed  independent  of  particle  size.  Using  Equation  92 
developed  by  Jonsson  (1966),  the  Reynolds  values  related  to  w^ave  runup  should  be 
greater  than  10'^.  Varying  the  model  time  periods,  from  0.5  to  2  seconds,  and  solving 
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Equation  92  for  A,  Table  1 1  gives  the  minimum  model  values  of  runup  according  to  this 
criteria.  The  actual  wave  periods  measured  in  the  model  ranged  from  about  0.5  to  0.8 
seconds.  Computed  minimum  values  of  runup  for  these  periods  on  a  2H.T  V  or  3H:1  V 
slopes  are  from  0.033  ft  (0.01  m)  to  0.062  ft  (0.019  m).  As  an  approximation,  runup  is  on 
the  order  of  3.0  times  the  actual  wave  height.  The  Shore  Protection  Manual  (1984)  and 
Broderick's  tabular  data  (Broderick  and  Ahren  1 982)  suggest  that  this  multiplier  is 
appropriate.  Calculations  regarding  wave  runup  will  be  discussed  later.  Therefore,  the 
minimum  acceptable  range  for  wave  heights  in  the  model  is  from  0.01 1  to  0.02  ft  (0.003 
to  0.006  m).  The  equivalent  prototype  runup  conditions  considered  valid  in  the  model, 
according  to  this  criteria,  for  a  3 -second  prototype  wave  period  are: 

a.  Tennessee-Tombigbee  -  For  2H:1V  slope,  H  >  1.16  ft  (0.35  m). 

b.  Navigation  Research  -  For  2H.T  V  slope,  H  >  1 .35  ft  (0.41  m)  and  for  3H.T  V 


slope,  H  >  0.89  ft  (0.27  m). 


Table  11 

Minimum  Values  of  Model  Wave  Runup  to  Satisfy  Reynolds  Criteria 

Model  Wave 
Period,  seconds 

Runup,  ft 
for  cot  a  =  2 

Runup,  ft 
for  cot  a  =  3 

Runup,  ft 
for  cot  a  =  4 

0.5 

0.049 

0.033 

0.025 

0.6 

0.054 

0.036 

0.027 

0.7 

0.058 

0.039 

0.029 

0.8 

0.062 

0.042 

0.031 

0.9 

0.066 

0.044 

0.033 

1.0 

0.070 

0.046 

0.035 

1.2 

0.076 

0.051 

0.038 

1.5 

0.085 

0.057 

0.043 

2.0 

0.098 

0.066 

0.049 
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Since  recommended  critical  values  of  Reynolds  numbers  range  from  1,000  to 
100,000,  a  minimum  value  was  sought  using  the  model  data  as  it  relates  to  the  stability 
number.  was  computed  using  the  maximum  wave  height,  ?  observed  during 
each  stability  test.  The  graph  in  Figure  35  shows  the  computed  values  of  particle 
Reynolds  numbers  for  all  TT  and  NR93  stability  data  for  a  2H:1V  slope. 


Figure  35.  Stability  number  as  a  function  of  particle  Reynolds  number,  2H:1V  slope 
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Nj  is  plotted  as  a  function  of  the  wave  Reynolds  number  in  Figure  36.  On  both  graphs, 
more  failed  points  are  observed  at  lower  Reynolds  values.  On  either  graph  a  clear 
relationship  does  not  appear,  such  that  could  be  adjusted  based  on  a  Reynolds  value 
as  Oumeraci  (1984)  suggested. 
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Figure  36.  Stability  number  as  a  function  of  wave  Reynolds  number,  2H;  IV  slope 


Angle  of  repose 

According  to  the  literature,  the  angle  of  repose  decreases  for  small  stone  sizes  and  as 
the  relative  height  of  the  revetment,  L/Djo ,  increases,  particularly  for  those  less  than 
approximately  2  in.  The  average  diameter,  D50 ,  of  the  model  stone  sizes  ranged  from 
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0.13  to  0.56  in.  (3.3  to  14  mm),  and  for  full  range  of  gradations  the  sizes  were  as  small  as 
0.092  in.  (2.3  mm)  and  as  large  as  1.0  in.  (25.4  mm).  These  variations  equate  to 
differences  in  (j)  from  the  smallest  to  largest  model  stones  of  3  to  5  deg  based  on  either 
the  "Bureau"  or  "Simons"  curve  in  Figure  10.  Differences  from  model  to  prototype  could 
be  as  high  as  10  deg. 

The  actual  length  of  the  model  slope  ranged  from  approximately  1.6  to  2.4  ft  (0.5  to 
0.7  m).  This  resulted  in  relative  revetment  heights  in  the  model  of  approximately  35 
to  220.  According  to  a  figure  relating  angle  of  repose  to  this  ratio  in  Maynord  (1988),  the 
angle  of  repose  value  at  L/D50  equal  to  35  was  45  deg  and  decreased  to  a  minimum 
datum  of  38.4  deg.  Values  were  not  given  for  ranges  of  L/D50  greater  than 
approximately  42. 

The  variation  in  ({)  between  the  largest  and  smallest  gradations,  from  either 
perspective,  equates  to  potential  impacts  from  the  angle  of  repose  on  stone  slope  stability. 
To  evaluate  these  impacts,  a  value  of  (j)  was  assigned  to  each  model  stone  gradation 
based  on  the  "Simons"  curve  in  Figure  10.  The  prototype  diameter,  model  diameter  and 
corresponding  angle  ([)  are  found  in  Table  12. 

The  stability  number  was  evaluated  as  a  function  of  the  angle  of  repose,  in  the  same 
manner  as  the  Reynolds  number.  Figure  37  shows  the  angle  of  repose  for  each  respective 
model  gradation  versus  the  stability  number  calculated  using  for  a  2H:  IV  slope. 

As  can  be  seen,  the  smaller  angle  of  repose  values  shows  more  failed  points  at  lower 
stability  numbers.  A  clear  delineation  or  relationship  between  (j)  and  cannot  be 
developed  without  more  data. 
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Table  12 

Model  Angle  of  Repose 


Figure  37.  Stability  number  as  a  function  of  the  angle  of  repose,  2H:1V  slope 
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Three  methods  were  evaluated  to  calculate  a  factor  for  adjusting  the  stone  diameter 
based  on  the  angle  of  repose  variations  between  gradations  in  the  model  stone.  Method  1 
is  based  on  the  shear  approach  presented  in  Equation  76  from  Blaauw  et  al.  (1984),  where 
the  diameter  is  a  function  of  kj .  In  method  2,  the  coefficient  is  based  on  Hedar’s 
formulation  for  a  downrushing  wave.  Equation  59,  where  D  is  a  function  of  K^own  >  ^ 
(j) ,  and  a  .  Finally,  method  3  employs  a  coefficient  from  Iribarren's  formulation  in 
Equations  97  and  98  containing  4>  •  These  methods  are  used  to  determine  if  an 
adjustment  can  or  should  be  made  to  the  model  stone  diameters. 

In  the  shear  approach  used  in  Method  1,  K,  is  assumed  to  equal  kj ,  and  computed 
according  to  Equation  77.  The  stone  diameter,  D ,  is  related  to  Kj  as 


(100) 

In  method  2  for  a  downrushing  wave  on  an  impervious  straight  slope,  the  stable  stone 
size  was  related  to  K^own 


D 


H 


f  ^ 


1 

16TC 

3 


K 


down 


(tan(j)cosa  -  sina) 


(101) 


Values  for  were  taken  from  Figure  65.3  in  Hedar  (1960). 
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Letting 


then, 


^2  = 


K 


down 


(tancf)  cosa  -  sina) 


D  «  K, 


(102) 


(103) 


In  method  3,  letting  K3  equal  in  Equation  97,  then 


D  «  K3 


(104) 


Table  13  contains  calculated  values  for  each  K  coefficient  at  two  slope  angles  and 
for  a  range  of  repose  angles  between  30  to  50  deg.  Generally,  the  K  factor  increases 
with  decreasing  angles  of  repose  and  increasing  bank  slope.  Since  these  values  are  pre¬ 
sented  as  directly  proportional  to  stone  diameter,  one  can  conclude  that  regardless  of 
method,  the  angle  of  repose  can  impact  effective  diameter.  It  is  also  shown  that  there  are 
larger  differences  in  the  K  coefficients  at  smaller  angles  of  repose,  for  example  between 
30  and  35  deg,  than  at  higher  angles,  45  to  50  deg.  Variations  in  K  are  smaller  for 
flatter  bank  slopes  than  for  steeper  slopes.  It  can  also  be  seen  that  the  shear-stress-related 
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Table  13 

Calculation  of  K  Values  for  Wide  Variation  in  Repose  Angle 

cot  a 

4) 

Kr’ 

K2 

Ks 

2 

30 

2.247 

52.57 

8.431 

2 

35 

1.600 

26.92 

3.926 

2 

40 

1.394 

19.83 

2.774 

2 

45 

1.292 

16.13 

2.240 

2 

50 

1.233 

13.60 

1.929 

3 

30 

1.289 

24.12 

2.487 

3 

35 

1.187 

19.21 

2.008 

3 

40 

1.148 

16.44 

1.746 

3 

45 

1.118 

14.37 

1.579 

3 

50 

1.098 

11.04 

1.462 

factors  are  smaller  than  the  factors  from  the  wave-related  methods.  Then  for  the 
experiments  where  (j)  varies  between  approximately  34  and  38  deg,  and  particularly  for 
the  2H:1V  slope,  the  angle  of  repose  should  be  considered  as  contributing  to  the  effective 
diameter  of  the  model  stone,  and  consequently,  to  its  stability. 

Of  the  three  methods,  it  was  assumed  that  methods  2  and  3  would  be  more  valid  since 
waves  were  considered  the  dominating  force  in  these  tests.  In  evaluating  the  factors 
based  on  Hedar's  formulation  (method  2),  the  factors  did  not  give  intuitive  results.  For 
instance,  for  angles  of  repose  from  45  to  50  deg,  the  difference  in  K  values  for  a  2H:1V 
slope  were  smaller  than  for  a  3H:  1 V  slope.  The  steeper  slope  should  experience  more 
instabilities  than  the  flatter  slope  for  the  same  change  in  (J> .  Therefore,  Iribarren's 
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formulation  (method  3)  was  assumed  appropriate  for  this  application.  Table  14  contains 
computed  values  of  K3  for  each  model  gradation  and  slope  tested. 


Table  14 

Values  of  K3  for  Each  Model  Condition 

Study 

^act 

K3 

cot  a=  2 

K3 

cot  a  =  3 

TT 

0.0463 

37.5 

3.22 

1.86 

NR 

0.0370 

36.6 

3.43 

1.91 

TT 

0.0338 

36.4 

3.49 

1.92 

NR 

0.0270 

35.9 

3.63 

1.95 

TT 

0.0263 

35.8 

3.66 

1.96 

NR 

0.0210 

35.6 

zni 

1.97 

GAL 

0.0200 

35.4 

3.79 

1.98 

NR 

0.0193 

35.2 

3.86 

1.99 

NR 

0.0108 

34.2 

4.26 

2.06 

To  compensate  for  differences  resulting  from  the  angle  of  repose,  the  model  diameters 
were  normalized  based  on  an  adjustment  factor  related  to  the  computed  values  of  K3 . 
Two  different  factors  were  generated,  one  for  adjusting  the  model  rock  sizes  to  the  largest 
model  size  and  one  for  adjusting  all  model  sizes  to  an  approximate  prototype  size.  Model 
diameters  were  adjusted  according  to 


K 

D  =  — 2E.  X  D 

CQ  2Ct 


K 


act 


(105) 
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where 


Deq  =  equivalent  stone  diameter,  ft  (m) 

Kmax  =  n^aximum  K3  value  for  normalizing  data 

Kact  “  actual  K3  value  corresponding  to  each  model  test  condition 

Dact  =  actual  stone  diameter  used  in  each  model  test,  ft  (m) 

The  equivalent  prototype  diameters  normalized  to  the  largest  model  stone  are 
represented  in  the  Table  15.  K„,ax  corresponds  to  a  37.5-deg  angle  of  repose.  The  actual 
diameters,  column  2,  represent  the  geometrically  scaled  gradations  used  in  the  model 
tests.  The  factors  in  columns  4  and  5  are  the  ratio  K„,ax/Kact  where  Ka,ax  for  a  2H:  IV 
slope  at  4)  equals  37.5  deg  is  3.22,  and  for  a  3H:1V  slope  is  1.86.  Values  for  ^nd 
Kact  ate  calculated  based  on  values  of  K3  in  Table  14.  Equivalent  diameters,  colunms  6 
and  7,  are  normalized  by  multiplying  column  2  by  columns  4  or  5,  respectively,  for  a 
2H:1V  slope  or  3H:1V  slope. 

Table  16  represents  adjusted  values  of  the  D50  stone  diameters  to  prototype  stone 
size  based  on  an  assumed  42-deg  angle  of  repose.  The  values  in  this  table  are  obtained 
using  the  K3  formulation  and  adjusting  the  values  accordingly  as  in  Table  15.  K,„ax  for 
a  2H:1V  slope  at  ^  equals  42  deg  is  2.52,  and  for  a  3H:1V  slope  is  1.671. 
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Table  15 

Equivalent  Prototype  D50  Normalized  Based  on  (j)  of  37.5  deg 


1 

Study 


^act 

prototype, 

in. 


4 

/Kact 

2H:1V 


5 

^^^ax 

3H:1V 


Deq 

2H:lVin. 


•^eq 

3H:1V 


Table  16 

Equivalent  Prototype  D50  Normalized  Based  on  (j)  of  42  deg 


3 

4 

5 

4> 

deg. 

/Kact 

^ct 

2H:1V 

3H:1V 

42.0 

1.00 

1.00 

37.5 

0.78 

0.90 

36.6 

0.73 

0.87 

36.4 

0.72 

0.87 

35.9 

0.69 

0.86 

35.8 

0.69 

0.85 

35.6 

0.68 

0.85 

35.4 

0.66 

0.84 

35.2 

0.65 

0.84 

34.2 

0.59 

0.81 

1 

Study 


'^act 

prototype, 

in. 


^eq 

2H:lVin. 


Deq 

3H:1V 
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Final  criteria  for  adjustment  of  model  data 


Several  approaches  could  be  taken  to  compensate  for  the  angle  of  repose  variations, 
Reynolds  effects,  and/or  boundary  layer  problems.  They  are  as  follows: 

a.  Data  could  be  eliminated  below  a  critical  Reynolds  criteria. 

b.  Data  associated  with  the  smallest  model  stone  gradation  could  merely  be  thrown 
out. 

c.  Model  stone  sizes  could  be  normalized  to  either  the  largest  model  size  or  to  an 
approximate  prototype  size. 

d.  The  model  stone  sizes  could  be  adjusted  and  then  a  Reynolds  criteria  applied. 

Based  on  the  investigation  of  Reynolds  criteria,  it  is  obvious  that  a  wide  range  of 

criteria  exist  for  evaluation  of  model  scale  effects  due  to  low  Reynolds  values.  Without 
doubt,  the  larger  scales  result  in  larger  particle  Reynolds  numbers  and  wave  Reynolds 
numbers.  LikeBroderickand  Ahrens  (1982),  in  this  research,  was  found  to  be  lower 
at  small  wave  heights  and  rock  sizes.  Since  Reynolds  criteria  did  not  explain  all  the 
20  percent  differences  observed  in  their  testing  and  was  uncertain  in  this  testing,  this 
method  was  not  considered  effective  by  itself  to  compensate  for  scale  effects. 

Removing  the  smallest  stone  has  a  two-fold  advantage.  It  leaves  larger  stones  that 
have  higher  angles  of  repose  (closer  to  prototype)  and  retains  threshold  points  in  the  data 
that  are  subject  to  higher  Reynolds  numbers.  Remember  that  the  smallest  diameter  model 
stones  were  not  valid  for  many  test  conditions  using  any  of  the  Reynolds  criteria. 
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However,  simply  removing  the  small  rock  data  does  not  eliminate  scale  effects  for  larger 
model  rock  sizes  and  removes  some  data  points  that  meet  Reynolds  criteria. 

Normalizing  diameters  to  the  maximum  model  size  indicates  that  observations  of 
failure  in  the  model  tests  provide  conservative  solutions.  Furthermore,  diameters 
normalized  to  a  "prototype"  value  would  indicate  that  the  size  of  the  prototype  stones 
could  be  as  much  as  40  percent  smaller  than  their  model  counterparts.  Due  to  insufficient 
data  regarding  model  to  prototype,  differences  in  the  angle  of  repose  were  not  adjusted  to 
"prototype"  values.  This  means  that  the  actual  stones  sizes  are  still  conservatively 
represented  in  the  model.  Adjusting  the  model  diameters  alone,  does  not  eliminate 
potential  Reynolds  effects,  it  simply  compensates  for  a  property  of  the  material  that  is  not 
generally  a  factor  at  the  prototype  scale. 

The  final  method  was  adopted  for  compensation  of  scale  effects.  First,  the  stone 
diameters  were  normalized  to  the  maximum  model  stone  size  as  in  Table  15  and  second, 
a  wave  Reynolds  criteria  was  applied  to  all  the  data.  All  stability  tests  that  corresponded 
to  Zjnaxj  H;  with  magnitudes  approximately  less  than  0.6  in.  (1 .5  cm)  in  the 

model  were  discarded.  If  any  of  the  three  model  wave  heights  met  the  criteria,  they  were 
kept  in  the  data  set.  Table  17  contains  the  final  adjusted  data  set,  in  prototype 
dimensions,  used  in  the  stability  analysis  for  the  2H:1  V  slope.  Table  18  contains  the  final 
adjusted  data  set  for  the  3H:  1 V  slope. 
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Table  17  (Continued) 
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Table  17  (Concluded) 


Table  18 

Final  Data  Set  for  3H:1V  Slope  in  Prototype  Dimensions 
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(Continued) 


Table  18  (Concluded) 


stability  Analysis 


Much  effort  has  been  made  to  minimize  scale  effects  due  to  both  low  Reynolds  values 
and  repose  angle  of  the  model  stone.  It  is  important  to  emphasize  in  the  deterministic 
design  of  bank  protection  for  tow-induced  impacts  that  the  critical  design  parameters  lie 
on  the  upper  bound  of  impacts.  Model  data  used  for  analysis  in  that  range  of  impacts  are 
the  least  effected  by  any  scale  problem.  Based  on  the  critical  examination  of  scale  effects 
and  the  adjustments  made  accordingly,  this  data  set  is  believed  to  be  valid  for  evaluating 
stone  slope  stability  due  to  tow-induced  impacts. 

The  experimental  approach  adopted  for  this  study  of  riprap  protection  for  tow- 
induced  forces  leads  us  to  an  integrated  solution.  That  is,  the  stable  bank  protection 
observed  from  this  testing  is  exposed  to  the  total  energy  imparted  by  the  passing  vessel 
including  pore  water  pressures  due  to  drawdown,  shear  stresses  due  to  velocities,  and  lift 
and  drag  forces  due  to  waves.  The  issue  then  becomes  what  general  equation  for  riprap 
design  can  be  used  considering  all  the  different  causal  mechanisms  acting  on  the  bank? 

There  are  basically  three  approaches  one  could  take  in  devising  a  methodology  for  the 
design  stone  slope  protection  from  the  integrated  testing.  First,  considering  that  the 
indirect  causes  of  bank  failure  are  the  vessel  speed,  its  hull  form,  the  shape  of  the  channel, 
and  the  sailing  line,  the  stability  of  stone  slope  protection  could  be  directly  related  to 
these  mechanisms.  The  advantage  of  this  approach  is  that  it  is  an  absolutely  integrated 
solution.  The  effects  are  related  to  characteristics  of  the  waterway  that  are  essentially 
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given.  There  is  no  added  error  due  to  prediction  of  measurement  of  forces  or  in  coupling 
the  causal  mechanisms.  The  problem  with  this  approach  is  that  considerably  more  data 
are  needed  to  develop  relationships  between  stability  and  say  blockage  ratio  or  boat 
speed.  This  approach  would  also  be  more  empirically  based  and  is  not  an  accepted 
methodology  like  the  Hudson  or  Isbash  equations.  Second,  an  integrated  solution  could 
be  developed  such  that  the  total  energy  of  the  direct  causal  mechanisms  (front  wave, 
gradients,  drawdown,  stem  wave,  return  current,  secondary  waves,  etc.)  is  related  to  the 
stability  of  the  stone.  Each  force  could  be  predicted  or  measured  and  plugged  into  a 
single  stability  equation.  The  disadvantages  of  this  approach  include  not  knowing  the 
relative  contributions  to  riprap  failure  of  each  force  component,  the  inability  to  account 
for  the  spatial  and  temporal  application  of  these  forces  on  the  slope,  and  the  variability  in 
the  equations  for  predicting  each  of  the  forces.  The  form  of  this  equation  is  not  known, 
but  it  is  probably  not  an  accepted  or  standard  solution.  The  third  approach  for  develop¬ 
ment  of  a  general  equation  for  riprap  stability  is  to  use  the  most  dominant  or  charac¬ 
teristic  force  or  forces  in  an  existing  stability  equation  with  empirical  coefficients 
designed  for  each  mechanism.  This  is  essentially  the  approach  taken  in  the  PIANC 
(1987)  guidance.  Evaluating  the  stone  size  from  either  a  return  current,  stem  wave,  or 
secondary  wave  equation  not  only  reduces  the  bias  of  looking  at  one  solution,  but  adds  a 
degree  of  safety  by  giving  the  designer  the  option  of  picking  the  most  conservative 
answer.  The  approach  is  relatively  simple.  Conversely,  the  results  or  empirical  coeffi¬ 
cients  developed  from  the  integrated  data  sets  are  likely  biased  because  the  parameters  are 
not  truly  separated.  Furthermore,  the  approach  may  oversimplify  the  problem.  In  a 
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related  approach,  one  or  more  variables  that  best  represent  an  integrated  solution  can  be 
selected.  Blauuw  et  al.  chose  a  “characteristic  stem  wave”  for  his  solution.  More  than 
one  characteristic  parameter  could  be  used  in  this  approach. 

The  third  method  of  stability  analysis  was  adopted.  Since  existing  methods  for  de¬ 
signing  riprap  are  based  on  either  a  shear  stress  approach,  velocity,  or  waves,  the  stability 
analysis  included  an  evaluation  of  these  approaches  with  several  dominant  mechanisms. 

Return  current 

An  attempt  was  made  to  use  the  shear  stress  approach  presented  by  Blaauw  et  al. 
(1984)  in  Equations  72-77  with  data  collected  from  this  research.  Observations  in  the 
model  did  not  comply  to  the  suggested  stability  criteria.  However,  this  may  be  due  to  the 
evaluation  of  the  drag  coefficient.  Based  on  information  from  Schlicting  (1979),  the 
equation  used  for  the  drag  coefficient  may  not  be  appropriate  for  gravel  and  riprap  size 
materials.  And  since  the  drag  coefficient  is  a  function  of  Reynolds  number,  it  may  not  be 
valid  for  the  conditions  tested  in  the  model. 

Both  Blaauw  et  al.  (1984)  and  PIANC  (1987)  ultimately  recommend  the  use  of  an 
Isbash  type  approach  for  evaluating  return  currents.  The  observed  data  are  compared  to 
Equation  78  in  Figure  38,  where  the  stone  size  is  related  to  the  maximum  return  currents. 
An  average  value  of  the  coefficient,  kj ,  of  1 .35  was  selected  for  the  analysis.  The  line 
represents  the  suggested  failure  threshold.  Points  above  the  line  should  be  stable  and 
below,  failed.  As  seen  by  the  figure,  all  observations  did  not  adhere  to  the  criteria.  That 
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Figure  38.  Analysis  of  return  current  using  Isbash  approach,  2H:1  V 

is,  many  failure  points  were  observed  above  the  Isbash  failure  threshold.  This  is 
probably  due  to  the  fact  that  the  return  current  does  not  cause  failure  of  the  stone  on  the 
upper  bank.  Other  mechanisms  are  involved  in  the  failure  process  that  are  not  captured  in 
this  approach. 

Waves 

Three  characteristic  waves  were  chosen  to  analyze  using  the  wave-based  stability 
approach  presented  by  Iribarren,  Hudson,  and  others.  They  were  the  maximum 


144 


drawdown  at  the  stem,  ?  the  maximum  secondary  wave,  H; ;  and  the  mayiTmim  stem 
wave,  Hnjax  •  Figures  39  and  40  show  how  the  final  adjusted  data  for  D50  and  2^^ , 
compare  to  Hudson's  equation  with  stability  coefficients  of  1.86  and  2.13  for  2H:1V  and 
3H:1V,  respectively.  Data  below  Hudson’s  threshold  would  indicate  potential  failure  of 
the  stone,  while  data  above  it  are  within  the  safe  threshold.  All  the  observed  failure 
points  fall  on  or  below  Hudson’s  safe  threshold. 

In  Figures  41  and  42  ,  the  peak  secondary  waves  associated  with  observed  data  points 
are  compared  to  Hudson’s  equation.  There  were  few  values  for  Hj ,  and  consequently, 
few  data  points  to  delineate  a  failure  threshold.  Recall,  secondary  wave  data  were  not 
available  for  the  TT  data.  Also,  if  any  wave  height  exceeded  the  minimum  required  to 
meet  the  Reynolds  criteria,  all  data  pertaining  to  that  data  set  were  retained  in  the  plots. 
Ignoring  values  of  Hi  less  than  approximately  1  ft,  the  failure  points  remaining  were 
still  below  the  threshold. 

Since  the  actual  height  of  the  stem  wave  can  differ  from  z^gx  =  H^ax  was  chosen  as  an 
alternative  for  a  characteristic  parameter  associated  with  the  stem  wave.  However,  both 
secondary  waves  and  the  maximum  stem  wave  showed  a  strong  relationship  to  Zn,ax  as 
was  shown  in  Figure  30.  Figures  43  and  44  show  the  relationship  of  the  observed  data  to 
Hudson’s  equation.  As  in  Figures  39  and  40,  the  separation  of  failed  and  stable  points 
adhered  to  the  Hudson  criteria.  Hudson’s  equation.  As  in  Figures  39  and  40,  the 
separation  of  failed  and  stable  points  adhered  to  the  Hudson  criteria. 
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Figure  39.  Analysis  of  using  Hudson  approach,  2H:1V 


Figure  40.  Analysis  of  using  Hudson  approach,  3H:  IV 
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Figure  41.  Analysis  of  H|  using  Hudson  approach,  2H:  IV 


Figure  42.  Analysis  of  H;  using  Hudson  approach,  3H:  IV 
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Figure  43.  Analysis  of  using  Hudson  approach,  2H:1  V. 


Figure  44.  Analysis  of  using  Hudson  approach,  3H:1V 
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In  an  attempt  to  couple  tow-induced  forces,  drawdown  and  return  current  were 
evaluated  as  a  total  head,  term,  Hj ,  such  that 


H 


T 


Z 


max 


(106) 


Replacing  wave  height  with  H-j-  and  comparing  the  values  for  from  Hudson  (1957), 
Figures  45  and  46  show  the  relationship  between  W50  and  Hj  for  each  slope,  respec¬ 
tively.  Since  the  velocity  head  due  to  return  current  is  small  relative  to  the  drawdown, 
little,  if  anything,  is  gained  in  using  this  method.  Separation  of  failed  and  stable  points 
does  not  appear  to  fit  the  cubic  relationship  of  the  Hudson  curve  any  better  than  using 
drawdown  alone. 


Selection  of  Ng 

The  next  step  in  the  analysis  is  to  determine  the  failure  threshold  for  the  data  given. 
Return  current  formulations  and  coupled  approaches  were  dismissed  from  the  analysis 
based  on  the  results  above.  Failure  thresholds  were  analyzed  as  a  function  of  ,  H,„ax , 
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Figure  45.  Combined  analysis  of  and  2^^  using  Hudson  approach,  2H:  1 V 


Figure  46.  Combined  analysis  of  Uj^  and  using  Hudson  approach,  3H:  1 V 


150 


and  Hj .  V alues  of  Ng  were  adjusted  in  each  of  the  plots  of  stone  weight  versus  wave 
height,  such  that  the  curve  passed  through  the  highest  failure  point.  That  is,  every  point 
at  and  below  the  failure  threshold  is  assumed  to  fail.  As  Ng  increases,  the  threshold 
shifts  downward,  indicating  that  the  required  stable  stone  size  is  smaller.  Figures  47-52 
show  failure  thresholds  through  these  data  for  ,  Hj ,  and  ,  respectively.  Data 
points  below  the  Reynolds  criteria  were  not  considered  in  the  selection  of  the  thresholds. 

These  values  of  Ng  do  not  represent  a  safe  design.  There  are  several  ways  to 
incorporate  safety  into  the  equation.  Ahren's  work  (Ahrens  1982)  resulted  in  such 
conservative  results  that  he  ultimately  recommended  a  method  using  the  damage  criteria 
where  zero  damage  is  defined  as  the  conditions  in  which  one  to  two  stones  are  moved. 
Broderick  and  Ahrens  (1982)  use  the  same  approach,  where  zero  damage  is  based  on  the 
movement  of  two  to  three  stones.  Hudson  (1957)  likewise  used  a  similar  approach  for 
zero  damage  plus  varied  the  significant  wave  height  to  determine  sensitivity.  If  there 
were  many  failure  points  along  the  threshold,  a  best  fit  line  can  be  used  to  determine  the 
failure  curve  and  safety  could  be  incorporated  by  .shifting  this  curve  through  the  outliers. 
This  method  was  used  by  Abt  et  al.  (1991)  in  the  design  of  riprap  grade  control  structures, 
The  recommended  safety  factor  resulting  from  that  methodology  and  for  a  subsequent 
model  study  at  the  WES  regarding  grade  control  structures  (Martin,  Knight,  and  Murphy 
1994)  was  1.25  times  Djg .  In  fact,  safety  factors  increasing  D50  have  been  used  by  a 
number  of  researchers.  Fuehrer,  Romische  and  Engelke  (1981)  increase  the  stable  stone 
size  from  their  tests  by  a  factor  of  3.  Maynord's  recommendations  in  Engineering 
Manual  1110-2-1601  (HQUSACE  1991)  was  a  multiple  of  1.1. 
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Figure  47.  Determination  of  N.,  for  2H:1V 
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Figure  49.  Determination  of  for  H,-,  2H:1V 


Figure  50.  Determination  ofNj  for  H;,  3H:1V 
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Figure  5 1 .  Determination  of  for  2H:  1 V 


An  approach  similar  to  the  Hudson  approach,  wherein  zero  damage  was  observed, 
was  used  in  this  analysis.  Values  of  the  zero  damage  stability  number,  ,  were 
selected  such  that  all  data  points  along  and  above  the  line  were  stable.  Observations  of 
"marginal"  had  to  fall  on  or  below  the  stable  threshold.  The  curves  shifted  through  the 
stable  points  are  also  shown  in  Figures  47-52.  The  safety  factor,  Sf ,  resulting  from  this 
analysis  can  be  calculated  by  the  ratio  of  N^/Ns^ .  Values  of  Nj  at  the  failure  threshold, 
at  zero  damage  and  Sf,  the  safety  factor  for  each  parameter  as  shown  in  Table  19. 


Table  19 

Final  Determination  of  Stability  Coefficients 

N. 

2H:1V 

N. 

2H:1V 

N. 

3H:1V 

3H:1V 

W  (/3 

< 

^ax 

1.82 

1.69 

1.08 

2.98 

1.71 

1.74 

2.72 

1.38 

1.97 

3.35 

1.54 

2.16 

1.89 

1.63 

1.16 

3.63 

1.94 

1.87 

Computed  safety  factors  ranged  from  1.08  to  2.16.  Safety  factors  used  in  open 
channel  flow  environments  can  be  typically  be  lower  than  in  breakwater  areas  or  near 
structures  since  human  safety  or  loss  of  property  would  be  minimal  in  these  areas.  The 
safety  factors  resulting  from  this  analysis  are  acceptable  for  this  application. 

For  the  2H:1  V  slopes  where  there  are  many  points  related  to  Hn,ax  and  ,  the 
values  of  and  Sf  are  similar.  This  is  not  a  surprise  since  both  parameters  are 

characteristic  terms  for  the  stem  wave.  These  results  also  follow  trends  in  the 
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hydrodynamic  analysis  suggesting  strong  correlations  between  the  two  parameters,  as 
was  shown  in  Figure  30.  In  fact,  it  follows  that  either  parameter  could  be  selected  to 
characterize  the  primary  wave  system  described  by  many  authors  as  the  dominant 
mechanism  for  failure.  The  3H:  1 V  curves  did  not  have  as  many  points  as  the  2H:  1 V 
curves  since  no  TT  data  were  taken  on  that  slope  angle.  For  and  ,  however, 
there  were  data  points  over  a  broad  range  of  wave  heights.  Stability  coefficients  in  the 
literature  (Table  1)  for  ranged  from  1.89  to  2.3.  The  zero  damage  stability 
coelBficients  from  this  analysis  for  z^^^  and  were  from  1.69  to  1.94.  Due  to  the 
similarities  between  z^^^  and  ,  values  of  were  averaged  for  each  slope 
yielding,  1.66  and  1.83  for  the  2H:1V  and  3H:1V  slopes,  respectively. 

The  values  of  related  to  the  secondary  waves  were  1.38  and  1.54  for  the  2H:1V 
and  3H:1  V  slopes,  respectively.  Values  in  the  literature  recommended  for  secondary 
waves  ranged  from  1.73  from  Sorensen  (1986)  to  2.4  from  PIANC  (1987).  Only  Blaauw 
etal.  (1984)  and  PIANC  (1987)  presented  formulations  for  both  Zn,ax  andHj.  Com¬ 
paring  those  values  of  stability  numbers  given  in  Table  1,  both  authors  suggest  that 
secondary  waves  have  higher,  not  lower,  stability  numbers  than  the  stem  wave.  The 
stability  numbers  related  to  Hj  were  dropped  from  the  analysis  based  on  the  following: 

a.  There  were  few  data  points  to  delineate  the  thresholds. 

b.  The  secondary  waves  were  implicitly  represented  by  the  stem  wave  parameters, 
since  secondary  waves  accompany  the  primary  wave  system  in  confined  reaches 
where  this  design  methodology  will  be  applied. 

c.  analysis  has  already  shown  the  strong  correlations  between  these  parameters. 
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d.  The  highest  measured  value  of  wave  height  was  less  than  2  ft  for  the  2H;  1 V  slope. 
/  Safe  values  of  the  stability  numbers  in  this  analysis  do  not  trend  in  the  same 
direction  as  the  literature  suggests. 

Slope  Angle  Relationship 


The  relationship  used  by  Hudson  (1957)  for  relating  stability  number  to  slope  angle, 
in  a  general  form  is 

=  a(cotcc)'/^  (107) 


Hudson's  stable  value  of  a  is  1.47.  He  states  in  his  paper  that  the  best-fit  line  was 
drawn  through  the  data  points  with  a  slope  of  one-third  to  simplify  the  formula  and  so 
that  it  could  be  plotted  as  a  straight  line  on  log-log  paper.  He  obtained  this  value 
experimentally  using  seven  different  slope  angles  between  cot  a  =  1.25  and  5.  PIANC 
(1987)  uses  this  formulation  and  the  coefficient,  a ,  is  1.5. 

Solving  this  equation  for  a  yields 


a 


(cota)'^^ 


(108) 


Then  calculating  a  based  on  values  of  associated  with  and  H^jg^  results  in  an 
average  a  of  1.30. 
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Figure  53  shows  the  relationship  between  the  stability  coefficient  and  the  slope  angle 


for  coefficients  of  a  from  Hudson,  PIANC  and  this  study. 


Figure  53.  Relationship  between  slope  angle  and  the  zero  damage  stability  number 

Since  tests  were  only  conducted  at  two  bank  slopes  in  this  study  and  since  the  data  points 
shown  in  Figure  53  follow  a  similar  slope  to  Hudson's,  the  zero  damage  stability  number 
for  this  analysis  is  given  by 


y 

=  1.30  (cot  a)  ^ 


(109) 
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Runup 


From  the  Shore  Protection  Manual  (1984),  the  runup  is  defined  by  Figure  54. 


Figure  54.  Definition  sketch  of  wave  runup 


The  manual  provides  a  method  for  calculating  runup  on  smooth  surfaces,  making 
corrections  for  scale  effects.  When  model  tests  for  riprapped  structures  are  not  available 
for  quantification  of  runup,  the  manual  also  provides  an  adjustment  for  riprap  from  the 
smooth  surface  estimate. 

In  subsequent  work  regarding  runup,  Ahrens  and  Heimbaugh  (1988)  present  an 
equation  for  runup  as  a  function  of  the  surf  parameter  and  slope  angle  based  on  two 
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laboratory  studies  on  riprap.  Their  equations  provide  a  simple  method  for  the  prediction 
of  runup  based  on  empirical  data.  Their  recommended  formulae  are: 

z  _  tan  a 

1/2 

(110) 

1.0  +  0.202 

where: 

is  the  surf  similarity  parameter  using  local  wavelength.  R^ax  the  mayimum  wave 
runup-  Hmo  is  the  energy-based,  zero-moment  wave  height.  Lp  is  the  Airy  wavelength 
calculated  using  Tp ,  the  period  of  peak  wave  energy  density  of  incident  wave  spectrum. 

Broderick  and  Ahrens  (1982)  definition  for  runup  is:  "Runup  was  visually  defined  as 
the  average  point  of  maximum  wave  uprush  on  the  riprap  surface  near  the  center  of  the 
wave  tank."  Small-scale  test  results  at  the  zero-damage  level  gave  more  conservative 
estimates  of  both  the  stability  number  and  the  runup.  In  fact,  they  conclude  that  runup  in 
the  small-scale  tests  is  approximately  20  percent  higher  than  prototype  measurements. 
From  their  tables,  it  appeared  that  the  magnitude  of  runup  was  on  the  order  of  three  times 
the  zero  damage  wave  height. 

Blauuw  et  al.  (1984)  define  the  upper  boundary  or  limit  of  riprap  protection  from  the 


imdisturbed  water  surface  with  the  wave  run-up,  R^ ,  formula 


R 


H 


2c 


-o.s 

tan  a 


(112) 


where  H  represents  secondary  waves,  L^j  is  the  wave  length,  and  c  is  a  coefficient. 
They  also  give  a  formula  to  determine  the  lowest  level  of  protection  from  secondary 
waves  below  the  undisturbed  water  level. 
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5  Summary  and  Conclusions 


Existing  methodologies  for  designing  stone  slope  protection  for  navigation-induced 
forces  produced  a  wide  variation  of  resulting  stone  sizes.  A  generalized  equation  for  the 
design  of  stone  protection  for  inland  navigable  waterways  in  the  U.S.  for  commercial  tow 
traffic  is  needed  to  optimize  safety  and  minimize  cost.  This  research  verifies  the  use  of 
the  Hudson  approach  to  riprap  design  for  waves  with  specific  recommendations  for  the 
zero  damage  stability  coefficient,  ,  for  tow-induced  forces  based  on  a  characteristic 
parameter  of  the  primary  wave  system,  maximum  drawdown. 

Furthermore,  lessons  learned  during  extensive  physical  model  testing  of  navigation 
effects  have  led  to  recommendations  for  future  model  studies  particularly  regarding  the 
scale  effects  at  low  Reynolds  values  and  related  to  the  stability  of  small  model  stones  due 
to  their  angle  of  repose. 

Recommendations  for  Physical  Model  Testing 

A  few  recommendations  are  suggested  regarding  physical  model  testing  of  riprap  in 
view  of  the  findings  of  this  research  regarding  scale  effects  at  low  Reynolds  numbers  and 
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for  small  stone  sizes.  Wave  Reynolds  values  should  be  approximately  greater  than  10"^ . 
This  conservatively  results  in  model  wave  heights  in  excess  of  approximately  0.6  in. 

(1.5  cm). 

Extremely  small  rock  sizes  should  be  avoided  unless  tests  can  be  conducted  to 
determine  not  only  the  relative  differences  in  angle  of  repose  between  gradations,  but  also 
the  importance  the  angle  of  repose  has  on  the  friction  coefficient  in  the  balance  of  forces. 
If  small  sizes  are  necessary,  adhere  to  the  method  prescribed  for  adjusting  model  sizes. 
Preferably,  angle  of  repose  should  be  measured  for  the  specific  model  gradation  tested. 
More  research  is  also  needed  to  verify  the  importance  of  the  friction  angle  in  the  balance 
of  forces. 

Finally,  regarding  navigation  effects  scale  testing,  adjustments  should  be  made  to  the 
tow  model  to  compensate  for  boundary  layer  problems.  Recommendations  regarding 
model  adjustments  can  be  found  in  Maynord  and  Martin  (1996)  based  on  the  most  recent 
comparisons  between  scale  testing  and  field  measurements. 

Design  Methodology  for  Sizing 
Riprap  for  Tow-induced  Forces 

The  completed  analysis  revealed  that  the  stem  wave  is  the  dominant  force  produced 
by  the  tow,  particularly  in  the  small  blockage  ratios  studied  here.  Secondary  waves  are 
small  or  negligible  for  slow-moving  vessels  and  are  of  the  same  order  magnitude  as  the 
stem  wave  for  laden  vessels.  There  were  no  tests  where  only  secondary  waves  were 
present.  In  the  case  where  tows  are  drafting  light  and  traveling  fast,  secondary  waves  can 
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be  of  greater  order  magnitude  than  the  primary  wave.  Since  this  analysis  integrates  these 
waves  with  the  primary  wave  system,  the  recommended  equation  conservatively  provides 
a  stable  rock  design. 

The  following  methodology  is  recommended  for  sizing  riprap  in  navigation  channels 
with  blockage  ratios  less  than  approximately  20  and  for  slopes  from  2H:1V  to  4H:1V: 

(1)  Determine  channel  dimensions  and  sailing  line. 

(2)  Select  the  critical  tow  dimensions  (the  largest  primary  wave  will  be  associated 
with  the  smallest  blockage  ratio).  The  largest  secondary  waves  will  be  associated 
with  the  highest  vessel  speed  and  perhaps  the  lightest  vessel  draft. 

(3)  Using  the  modified  Schijf  equations  presented  in  Maynord  (1996),  calculate  the 
vessels  limiting  speed  given  these  waterway  and  tow  characteristics. 

(4)  If  an  actual  maximum  vessel  speed  is  not  known  for  this  channel  reach,  select  a 
vessel  speed  of  approximately  95  percent  of  the  limiting  speed  to  calculate  the 
stem  wave  parameter,  » using  the  methods  in  Maynord  (1996)  or  PIANC 
(1987). 

(5)  Calculate  the  secondary  wave  height,  H; ,  using  the  methods  in  PIANC  (1987)  or 
in  Verhey  and  Bogaerts  (1989). 

(6)  Select  a  bank  slope.  Use  the  following  equation  to  determine  for  that  slope 

-  1.30  (cot  a)‘'^  (113) 
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(7)  Compare  to  H;  and  select  the  largest  value  for  design.  Using  the  maYimum 
wave,  H ,  calculate  the  stone  size  as  follows: 


D 


so 


N. 


H 


(114) 


(8)  Determine  stone  weight,  Wjq  ,  based  on  a  spherical  stone  diameter  such  that 


W 


so 


(115) 


Other  Design  Considerations 

There  are  other  considerations  in  the  design  of  riprap  in  addition  to  its  size.  Existing 
guidance  should  be  used  to  determine  quality,  gradation,  placement,  and  limits  of  bank 
protection. 

The  stability  of  the  riprap  is  dependent  upon  the  gradation,  placement,  and  thickness 
of  the  stone.  Stone  should  be  well  graded,  and  the  thickness  should  not  be  less  than 
I^ioo  •  Specific  guidance  for  determination  of  these  parameters  and  others  can  be  found  in 
Engineering  Manual  1110-2-1601  (HQUSACE  1991). 

In  a  riverine  situation,  pool  fluctuations  and/or  flood  stages  may  dictate  the  upper 
limit  of  bank  protection.  However,  if  the  waterway  is  being  designed  for  tow-induced 
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waves,  then  the  upper  extents  of  stone  protection  should  reflect  the  maximum  navigable 
pool  or  stage  plus  the  maximum  potential  wave  runup.  Generally,  protection  will  extend 
from  this  upper  limit  to  the  toe  of  the  slope.  Runup  should  be  calculated  using  guidance 
in  the  Shore  Protection  Manual  or  with  Ahren's  equations  until  further  research  verifies 
runup  values  appropriate  for  boat  waves  in  shallow  water  in  navigable  waterways. 

As  was  discussed  in  the  literature  review,  some  researchers  suggest  that  different  size 
stone  protection  may  be  used  for  upper  and  lower  banks  along  the  waterway.  In  that  case, 
the  upper  bank  protection  for  tow-induced  waves  should  extend  from  the  upper  limit  of 
protection  to  below  the  water  level  at  the  lowest  expected  navigable  conditions,  taking 
into  account  ^e  maximum  expected  drawdown.  If  stone  protection  is  required  for  the 
upper  slope,  the  lower  portion  of  slope  should  not  be  left  unprotected.  Studies  conducted 
at  the  WES  for  Himtington  District  suggest  that  unprotected  lower  slopes  could 
experience  failure  due  to  tow-induced  currents,  particularly  if  tows  travel  close  to  the 
bank.  More  research  is  needed  regarding  placement  of  different  rock  sizes  on  upper  and 
lower  slope. 

Location  of  the  protection  should  be  considered  when  selecting  the  hydrodynamic 
conditions.  In  lock  approaches,  it  is  highly  unlikely  that  tows  will  travel  near  the  limiting 
speed.  Reasonable  operating  conditions  should  dictate  design  methodologies  in  these 
areas.  In  the  case  of  long,  narrow  reaches  of  natural  rivers  or  man-made  canals  where 
vessels  underway  can  approach  the  limiting  speed,  critical  hydraulic  loadings  based  on 
approximately  95-percent  limiting  speed,  should  be  selected  for  bank  protection  design. 

If  wind  waves,  ambient  currents,  or  propeller  jets  could  potentially  impact  the  bank. 
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hydrodynamic  forces  should  be  estimated  for  each  of  the  conditions  using  available 
literature  (the  PIANC  guidelines  provide  an  acceptable  start)  and  stone  sizes  should  be 
estimated  using  riprap  equations  developed  for  those  conditions.  If  the  analysis  produces 
different  stone  sizes,  the  maximum  stone  protection  should  be  selected. 

In  design  of  the  waterway,  attention  should  be  given  to  expected  and  potential 
mooring  locations.  Observations  by  the  author  during  a  bank  erosion  field  survey  on  the 
Upper  Mississippi  and  Illinois  Rivers  suggested  that  the  integrity  of  the  bank  can  be 
jeopardized  at  locations  where  barge  mooring  occurs.  Once  a  nick  point  has  been  made 
in  the  bank  protection  (whether  rock  or  vegetation),  the  bank  becomes  subject  to  further 
unraveling  from  other  mechanisms.  This  can  result  in  dynamically  eroding  reaches  of 
bank  along  the  waterway.  Care  should  be  taken  to  either  prohibit  mooring  in  critical 
reaches,  design  specific  areas  for  mooring  of  barges,  and/or  be  prepared  to  repair  areas  of 
failure  through  routine  maintenance. 
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6  Current  and  Future  Studies  in 
Navigation  Effects 


Navigation  Effects  Studies 

The  WES  is  currently  funded  to  investigate  the  physical  effects  due  to  projected 
incremental  increases  in  navigation  on  the  Upper  Mississippi  River  System  -  Illinois 
Waterway  (UMRS)  as  a  part  of  the  environmental  portion  of  a  feasibility  navigation 
study.  The  environmental  concerns  for  this  project  regard  the  impacts  of  navigation  on 
adult  fish,  fish  larvae,  aquatic  vegetation,  and  benthic  organisms.  Of  primary  concern  is 
the  resuspension  and  transport  of  sediments  into  critical  biological  areas  such  as 
backwaters  or  side  charmels  adjacent  to  the  main  channel  waterways.  Also  important  are 
the  hydrodynamic  forces  in  the  main  channel  that  potentially  could  disturb  the  biota, 
including  the  mixing  due  to  tow  passage,  the  propeller  jet  forces,  the  wave  energy,  and 
the  current  reversals. 

A  new  funding  start  from  HQUSACE  in  fiscal  year  1995  allowed  the  research  and 
testing  to  continue  into  the  more  complex  phenomena  regarding  navigation  effects.  This 
work  has  focused  on  the  engineering  aspects  of  this  subject,  and  its  research  is  directed 
toward  design  of  chaimel  protection  in  navigation  charmels.  Particular  emphasis  for  the 
research  overlaps  the  objectives  of  the  UMRS  study  including  the  development  of 
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predictive  methods  for  return  currents  and  wave  heights  in  irregular  channels  with 
ambient  current  conditions.  Design  guidance  regarding  stable  riprap  will  also  be 
developed  for  propeller  jets  and  maneuvering  tows.  A  tow-induced  sediment 
resuspension  algorithm  will  be  investigated  for  both  studies  and  its  potential  inclusion  in 
a  2-D  sediment  transport  model. 

Developing  Technologies 

All  the  navigation  effects  studies  will  benefit  from  better  physical  model  facilities, 
relatively  new  technologies  in  instrumentation,  and  the  development  of  numerical  codes 
that  can  model  the  tow  movement.  Both  physical  and  numerical  model  output  will  be 
used  for  visualization  of  the  physics  of  the  problem  and  to  provide  a  better  understanding 
of  the  effects  of  navigation. 

A  new  flume  for  the  UMRS  Navigation  Effects  Study  provides  ample  length  and 
width  to  conduct  experiments  that  were  otherwise  not  possible  in  the  old  facility.  The 
flume  is  400  ft  (122  m)  long  by  70  ft  (21.3  m)  wide  by  4  ft  (1.2  m)  deep,  is  equipped  with 
a  towing  carriage,  and  has  a  flow  capacity  of  50  cfs  (1.4  cms).  Channels  with  irregular 
cross  sections  and  ambient  currents  can  be  modeled  in  the  facility. 

State-of-the-art  data  acquisition  and  instrumentation  are  available  on  this  model.  It 
is  equipped  with  ten  acoustic  doppler  velocity  probes  that  are  capable  of  taking  non- 
intrusive  3-D  velocities,  a  2-D  laser  doppler  velocimeter,  and  wave  capacitance  rods  for 
water  surface  measurements. 
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Recent  numerical  model  modifications  of  an  existing  hydrodynamic  code,  HIVEL2D 
(Stockstill  and  Berger  1992),  have  resulted  in  a  product  that  can  be  used  to  simulate  tow 
movement  through  a  channel.  This  code  has  been  verified  with  physical  model  data  so 
that  it  can  be  used  as  the  hydrodynamic  driver  for  the  sediment  transport  studies  and  as  a 
tool  for  predicting  velocities  and  drawdown  due  to  incremental  increases  in  navigation  in 
reaches  of  the  UMRS  (Stockstill,  Martin,  and  Berger  1995).  Because  it  is  a  2-D  model,  it 
is  only  valid  for  evaluating  far-field  effects,  areas  away  from  the  boimdaries  of  the  tow. 

The  products  of  these  studies  vrill  be  a  family  of  tools  that  will  predict  forces 
produced  by  navigation  and  provide  design  guidance  for  channel  protection  in  navigable 
waterways.  These  studies  will  help  in  the  avoidance  and  mitigation  of  environmental 
impacts  due  to  commercial  traffic  on  inland  waterways.  These  studies  have  already  and 
will  continue  to  advance  the  state  of  knowledge  regarding  navigation  effects. 
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Appendix  A 
Definitions 


Ambient  current:  The  river  current  undisturbed  by  the  presence  of  a  vessel. 

Blockage  ratio:  The  ratio  of  the  channel  cross-sectional  area  to  the  submerged  cross- 
sectional  area  of  the  vessel. 

Bottom  displacement  current:  The  current  beneath  the  vessel  acting  in  the  opposite 
direction  to  the  movement  of  the  vessel. 

Bow  current:  The  current  moving  ahead  of  the  vessel  in  general  direction  of  the  vessel. 

Bow  wave:  The  wave  generated  at  the  bow  of  the  vessel. 

Drawdown:  The  water  level  depression  that  forms  as  the  vessel  moves  forward  and  water 
is  displaced  from  bow  to  stem.  Drawdown  is  a  function  of  vessel  speed,  vessel 
size,  and  channel  geometry. 

Limiting  speed:  The  maximum  attainable  speed  of  a  vessel  in  a  confined  waterway  or  its 
critical  speed.  This  speed  is  limited  by  the  blockage  ratio  not  the  vessel 
horsepower. 

Primary  waves:  The  long  period  disturbance  caused  by  the  displacement  of  water  from  a 
vessel  moving  through  the  waterway.  This  disturbance  includes  the  front  wave, 
drawdown,  return  current,  and  transversal  stem  wave. 
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Propeller  jet:  The  highly  three-dimensional  currents  associated  with  propeller  jets  which 
cause  localized  disturbances  to  the  flow  field. 

Return  current:  The  current  adjacent  to  the  vessel  that  accompanies  the  drawdown  as  the 
water  moves  from  bow  to  stem  during  the  forward  motion  of  the  vessel.  This 
current  acts  in  the  direction  opposite  of  tow  movement  and  generally  parallel  to 
the  bank. 

Sailing  line:  The  path  of  the  vessel  in  the  navigation  channel.  Usually  stated  as  its 
distance  from  the  bankline. 

Secondary  waves:  The  waves  generated  by  a  moving  pressure  disturbance  in  deep  water. 
Transverse  and  diverging  waves  intersect  to  form  interference  peaks  or  ship 
waves. 

Slope  supply  flow:  The  current  produced  as  water  fills  in  behind  the  stem  to  replace  the 
water  displaced  as  the  vessel  moves  forward. 

Transversal  stern  wave:  The  transition  at  the  stem  of  the  vessel  from  the  drawdown  to 
the  normal  water  level  following  the  vessel. 
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